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Introduction 


Objective 

This  report  documents  the  results  of  salinity  intrusion  analyses  for  the 
Panama  Canal,  Panama.  The  objective  of  these  analyses  is  to  provide  an  early 
estimate  of  changes  in  salinity  intrusion  and  freshwater  usage  that  may  result 
from  an  expansion  of  canal  locks  capacity,  and  to  identify  measures  that  may 
reduce  salinity  intrusion  and  freshwater  consumption  if  unacceptable  increases 
are  forecast. 


Background 

The  Panama  Canal  is  an  80-kilometer  (50-mile)  long,  man-made  waterway 
connecting  the  Atlantic  Ocean  (Caribbean  Sea)  and  Pacific  Ocean  (Bay  of 
Panama)  for  purposes  of  navigation  (Figure  1-1).  It  became  operational  on 
August  15,  1914.  Ship  traffic  through  the  canal  has  steadily  increased  in  terms  of 
size  and  frequency  and  currently  about  thirty-eight  ships  transit  the  canal  each 
day.  The  canal  between  the  two  oceans  passes  through  two  natural  fresh  water 
lakes — Gatun  Lake  and  Miraflores  Lake  (see  Figure  1-2). 

Fresh  water  from  Gatun  and  Miraflores  lakes  is  used  for  filling  the  naviga¬ 
tion  locks  while  transiting  ships  in  and  out  of  the  lakes.  Salt  water  from  the 
ocean  gets  added  to  the  lakes  by  density  currents  and  fresh  water  from  the  lakes 
is  lost  to  the  sea  during  transit  of  ships.  The  net  loss  of  fresh  water  is  estimated 
by  the  Panama  Canal  Commission  to  be  about  fifty-two  million  gallons  per 
lockage  operation.  Fresh  water  is  also  supplied  to  the  local  population  in  two 
cities  for  drinking  purposes  at  the  rate  of  about  forty-seven  million  gallons  per 
day.  During  years  of  low  rainfall,  a  shortage  of  fresh  water  is  experienced  for 
canal  operations. 

The  Panama  Canal  Commission  is  considering  expanding  the  existing  ship 
transit  capacity  of  the  canal,  which  may  further  increase  freshwater  consumption 
and  increase  salt  intrusion  from  the  ocean.  The  Commission  needs  an  estimate  of 
the  potential  adverse  effects  of  capapacity  expansion  and  an  early  indication  if 
remedial  measures  can  be  expected  to  provide  protection  from  those  impacts. 
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Approach 

The  objectives  were  achieved  by  combining  desktop  reviews  and  analytic 
tools  with  simple,  unverified  numerical  models  that  provided  rapid  estimates  of 
salinity  intrusion  through  the  locks  and  freshwater  usage.  The  following  tasks 
were  performed:  data  on  ship  transits  and  observed  salinities  assembled  by 
Panama  Canal  Commission  staff  were  compiled  and  analyzed;  a  literature  review 
identified  potential  salinity  mitigation  methods;  spreadsheet  calculations 
characterized  salt  loading  and  water  usage;  a  two-dimensional  width-averaged 
(2DV)  numerical  model  was  run  to  define  density-driven  flow  out  of  the  upper 
locks;  and  a  two-dimensional  depth-averaged  (2DH)  numerical  model  was  run  to 
depict  the  dispersion  path  of  salt  water  if  it  were  released  into  Gatun  T  ake.  The 
numerical  models  were  applied  without  verifying  them  to  field  observations,  an 
approach  which  is  satisfactory  for  demonstrating  feasibility  of  plans.  If  a  canal 
capacity  expansion  plan  is  selected  for  design,  these  models  should  be  verified  to 
field  data  before  testing  the  design  plan. 

Three  conditions  -  existing  conditions,  use  of  holding  ponds  to  recirculate 
lock  lift  water  and  conserve  fresh  water,  and  use  of  a  Syncrolift  lock  (see  Part  2) 
to  increase  transit  capacity  -  were  examined,  and  a  number  of  salinity  intrusion 
mitigation  measures  were  examined  for  feasibility.  Model  results  for  two 
mitigation  methods  are  presented  here  -  a  salt  water  sump  and  a  gated  transition 
chamber. 
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2  The  Panama  Canal 


Design 

The  80  km  (50  mile)  long  man-made  Panama  Canal  connects  the  Atlantic 
Ocean  on  the  North  side  and  the  Pacific  Ocean  on  the  South  side  (Figure  1-2). 
On  its  way  it  passes  through  Gatun  Lake  and  Miraflores  Lake,  in  which 
controlled  water  levels  are  maintained  at  26  m  (85  feet)  and  16.5  m  (54  feet) 
above  mean  sea  level,  respectively.  During  high  rainfall  season  Gatun  Lake  is 
maintained  at  a  level  of  26.7  m  (87.5  feet). 

Both  the  Gatun  Lake  and  Miraflores  Lake  are  fresh  water  lakes.  Gatun  Lake 
has  an  estimated  storage  of  5.2  billion  cubic  m  (183  billion  cubic  feet)  at  a  lake 
level  of  26  m  (85  feet)  and  5.56  billion  cu  ft  (196  billion  cubic  feet)  at  an 
elevation  of  26.7  m  (87.5  feet).  It  has  a  surface  area  of 435  sq  km  (168  square 
miles)  at  an  elevation  of  26.7  m  (87.5  feet).  Gatun  Lake  was  created  by  placing 
an  earthen  dam  across  the  Chagres  River.  The  Gatun  Dam  is  nearly  2.4  km 
(1.5  mile)  long  and  0.8  km  (half  mile)  wide  at  the  base,  sloping  to  a  width  of 
30.48  m  (100  feet)  at  the  crest,  which  is  32  m  (105  feet)  above  the  sea  level,  or 
6  m  (20  feet)  above  the  normal  level  of  Gatun  Lake.  Another  series  of  dams  was 
constructed  under  the  Madden  Dam  Project,  which  stores  rain  water  from  the 
catchment  area  of  about  3367  sq  km  (1300  sq  miles).  The  main  features  include 
the  main  concrete  dam,  a  hydro  power  station  and  several  earth  and  gravel-fill 
dams.  Water  level  in  the  Gatun  Lake  is  regulated  by  releasing  excess  water  to 
the  sea  over  the  spillway  of  Gatun  Dam  and  by  letting  in  water  from  the  Madden 
Dam.  Water  level  in  Miraflores  Lake  is  regulated  through  the  Miraflores  Dam, 
which  is  a  gravity  concrete  dam.  The  dam  has  a  spillway  with  a  design  capacity 
of  2832  cu  m/s  (100,000  cfs)  at  the  highest  allowable  level.  Miraflores  Lake  is 
much  smaller  compared  to  the  Gatun  Lake.  It  has  a  surface  area  of  3.94  sq  km 
(1.52  sq  miles)  at  an  elevation  of  16.5  m  (54  feet).  Available  active  storage  of 
Miraflores  Lake  between  elevations  16  m  and  17  m  (53  and  55  feet)  is  approxi¬ 
mately  2.5  million  cu  m  (88  million  cu  feet)  and  the  total  storage  is  on  the  order 
of  1 9  million  cu  m  (670  million  cubic  feet). 

Fresh  water  from  these  two  lakes  is  used  for  filling  the  navigation  locks 
while  transiting  ships  through  the  canal.  Salt  water  from  the  ocean  gets  added  to 
the  lakes  during  transit  of  ships  from  the  ocean  to  the  lakes,  and  fresh  water  from 
the  lakes  is  lost  to  the  sea  while  transiting  ships  from  the  lakes  to  the  sea.  The 
net  loss  of  fresh  water  is  estimated  to  be  about  0.20  million  cu  m  (52  million 
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gallons)  per  lockage  operation.  Fresh  water  from  Gatun  Lake  is  also  supplied  to 
the  local  population  in  two  cities  for  drinking  purposes  at  the  rate  of  about 
0.18  million  cu  m  (47  million  gallons)  per  day.  During  years  of  low  rainfall,  a 
shortage  of  fresh  water  is  experienced  for  canal  operations. 

The  Panama  Canal  has  six  twin  locks,  grouped  in  three  sets,  namely 
Miraflores  Locks  (2  lifts)  and  Pedro  Miguel  Lock  (1  lift)  on  the  Pacific  side  and 
Gatun  Locks  (3  lifts)  on  the  Atlantic  side  (Figure  2-1).  The  three  navigation 
locks  on  each  end  lift  ships  from  the  ocean  to  the  Gatun  Lake  through  a  height  of 
.26  m  (85  feet)  above  sea  level  and  then  lower  them  through  three  locks,  bringing 
them  back  to  the  ocean  level.  The  channel  that  connects  Gatun  Lake  to  the 
Miraflores  Lake  is  called  Gaillard  Cut.  Ships  travel  for  about  37  km  (23  miles) 
from  Gatun  Locks  to  the  north  end  of  Gaillard  Cut.  The  twin  locks  serve  two 
purposes;  one  as  a  stand-by  in  case  of  closure  of  one  system  of  locks  and  the 
second  to  serve  as  an  alternate  system  of  locks  for  transiting  multiple  small  ships 
in  a  single  transit. 

Physical  characteristics  of  the  six  locks  are  given  in  Table  2-1  below. 


Table  2-1 

Lock  Characteristics 

Width  &  Length 

Lift  Height 

Lift  Volume 

Bed  Level 

Lock 

(ftV(m) 

mm) 

(cft)/(cu  m) 

(ft)/(m) 

Lower  Miraflores 

110  x  1051  ft 

26  ft 

3.005  x  106 

-52  ft 

Lock 

33.5  x  320.3  m 

7.9m 

0.851  x  10s 

-16  m 

Upper  Miraflores 

110  x  1071  ft 

28  ft 

3.298  x  106 

-20.3  ft 

Lock 

33.5  x  326.4  m 

8.5  m 

0.934  x  105 

-  6.2  m 

Pedro  Miguel 

110  x  1071  ft 

31  ft 

3.652  x  106 

+  11.0ft 

Lock 

33.5  x  326.4  m 

9.5  m 

1.034  x10s 

+  3.4  m 

Upper  Gatun  Lock 

110  x  1071  ft 

28  ft 

3.298  x  106 

+  13.0ft 

33.5  x  326.4  m 

8.5  m 

0.934  x  10s 

+  3.9  m 

Middle  Gatun 

110  x  1051  ft 

29  ft 

3.352  x  10s 

-15.3  ft 

Lock 

33.5  x  320.3  m 

8.8  m 

0.949  x  10s 

-  4.7  m 

Lower  Gatun  Lock 

110  x  1051  ft 

28  ft 

3.257~x  1U6 

-44.3  ft 

33.5  x  320.3  m 

8.5  m 

0.917  x  10s 

-13.5  m 

Balboa  Port  is  located  on  the  Pacific  Ocean  side  and  the  Port  of  Cristobal  is 
located  on  the  Atlantic  side.  Two  breakwaters  are  constructed  on  the  Atlantic 
side  to  form  Limon  Bay,  which  provides  shelter  for  ships  against  waves. 


Operation 

The  canal  operation  consists  of  two  distinct  parts,  namely  raising  a  ship  and 
lowering  a  ship.  Water  from  Gatun  Lake  is  released  into  the  higher  Gatun  Lock 
and  in  the  Pedro  Miiguel  Lock  for  lifting  ships  in  these  locks.  Water  from 
Miraflores  Lake  is  released  into  upper  Miraflores  lock  and  from  there  to  the 
lower  Miraflores  lock  for  lifting  ship  from  a  lower  elevation.  Large  culverts  are 
provided  for  filling  and  emptying  the  locks.  A  safe  rate  of  change  in  water 
elevation  in  a  lock  is  dependent  on  the  distribution  of  water.  For  the  Panama 
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Canal  locks  a  maximum  rate  of  2.3  m  (7.5  ft)  per  minute  with  a  good  water 
distribution  is  considered  optimum.  It  takes  about  10  minutes  to  fill  or  empty  a 
lock. 

A  ship  has  to  transit  through  six  locks  and  also  sail  through  an  80.5  km 
(50-mile)  long  navigation  channel.  It  requires  about  nine  hours  for  an  average 
ship  to  transit  the  canal.  Larger  ships  are  transited  during  daylight  hours  whereas 
smaller  vessels  are  transited  at  any  convenient  time.  Hydraulics  of  the  Panama 
Canal  Locks  is  described  in  detail  in  a  paper  by  Whitehead  (1915). 


Traffic 

Although  85  years  old,  the  strategically  located  canal  is  still  one  of  the  most 
highly  traveled  waterways  of  the  world.  It  operates  24  hours  a  day,  365  days  a 
year  transiting  vessels  of  all  nations.  Since  the  opening  of  the  waterway  to  world 
commerce,  the  canal  has  provided  transit  service  to  more  than  700,000  vessels. 
Toll  collection,  which  has  reached  about  500  million  dollars  a  year,  is  a  major 
source  of  income  for  the  Republic  of  Panama.  The  physical  dimensions  of  locks 
place  a  restriction  on  the  maximum  size  of  ship  that  can  transit  through  the  canal. 
This  size,  popularly  known  as  the  Panamax  ship,  is  shown  in  Figure  2-2.  The 
maximum  length,  width  and  draft  for  the  Panamax  ship  are  294  m  (965  ft), 

33.5  m  (1 10  ft)  and  12  m  (39.5  ft)  respectively. 

The  number  of  ships  that  transited  through  the  canal  during  1995  through 
1998  is  shown  in  Figures  2-3  through  2-6. 


Salinity 

Very  meager  data  are  available  on  the  salinity  of  water  in  the  two  lakes  and 
inside  and  outside  locks.  Historical  data  on  salinity  are  not  available  probably 
because  in  the  past  years  an  increase  in  salinity  of  lake  water  was  neither  noticed 
nor  threatened  under  normal  canal  operation.  Only  during  the  recent  years 
concern  over  possible  increase  in  lake  salinity  was  perceived  due  to  increase  in 
ship  transits  per  day  and  low  fresh  water  inflows.  It  is  learned  from  verbal 
communication  from  PCC  officials  that  the  salinity  of  Miraflores  Lake  increased 
quickly  immediately  after  commissioning  of  the  canal,  which  lead  to  changing 
the  location  of  intake  for  supply  of  fresh  water  to  the  city.  Salt  accumulation  was 
noticed  on  the  pipes  of  the  power  station.  A  steady-state  salinity  was  however 
reached  shortly  (maybe  within  the  first  year  or  so)  and  the  salinity  of  Miraflores 
Lake  has  been  less  than  one  part  per  thousand  ever  since.  Salinity  of  Gatun  Lake 
is  known  to  be  close  to  or  equal  to  zero.  It  is  reasonable  to  assume  that  the  salt 
intake  per  ship  transit  is  balanced  by  salt  output  from  both  the  lakes  resulting 
from  the  freshwater  released  per  transit  from  the  Gatun  Lake. 

Field  data  on  salinity  collected  by  PCC  are  summarized  in  Table  2-2  below. 
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Table  2-2 

Observed  Salinities 

Location  and  Time 

WSSSKKKM 

Pacific  Ocean  just  outside  LI 

10.0 

All  times. 

Inside  LI 

8.0 

Before  lift. 

Inside  L2 

1.5 

After  moving  ship  from  LI  to  L2 

Miraflores  Lake 

1.0 

All  times 

Proposed  Expansion 

The  present  system  has  reached  its  maximum  capacity  and  there  is  no  scope 
for  increasing  the  number  of  ships  transiting  the  canal  per  day.  In  order  to 
increase  transits,  the  Panama  Canal  Commission  is  considering  various 
alternatives.  Construction  of  an  additional  system  of  locks  is  a  very  expensive 
option.  Providing  a  Syncrolift  type  of  locks  at  both  ends  of  the  canal  is  being 
considered.  See  description  below. 


Syncrolift 

Raymond  Pearlson  of  Miami,  Florida,  USA,  invented  the  Syncrolift  shiplift 
system  in  1954.  In  1979,  the  original  company,  named  Pearlson  Engineering, 
became  Syncrolift  Inc,  which  is  a  subsidiary  of  Rolls  Royce,  UK.  The  system 
essentially  consists  of  a  large  platform,  which  can  be  lowered  into  water,  a  ship 
positioned  over  it  and  then  the  platform  raised  vertically  above  water  along  with 
the  ship.  It  has  three  principal  components:  (1)  a  structural  steel  platform, 

(2)  electrically  powered  wire  rope  hoists  for  raising  and  lowering  the  platform, 
and  (3)  an  electric  motor  system  to  operate  the  system.  Ships  are  dry  docked  on 
the  platform  using  conventional  procedure  for  setting  the  keel  blocks  and  placing 
the  bilge  blocks.  Ships  raised  above  water  are  transported  horizontally  over  rails 
to  the  desired  location  on  land.  The  first  Syncrolift  was  built  in  Miami,  Florida 
in  1957,  and  it  is  still  operational.  By  1998,  220  Syncrolifts  had  been  installed  in 
66  countries  with  lifting  capacities  from  101,600  kg  (100  tons)  to  60.9  million  kg 
(60,000  tons  ). 

The  advantages  of  Syncrolift  over  a  conventional  lock  system  are: 

Lower  initial  cost 
Shorter  construction  time 
Less  susceptible  to  biological  fouling 
-  No  under  water  work 
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No  use  of  fresh  water  for  raising  ship.  However,  freshwater  may  be 
lost  through  density  current  exchange  when  the  lock  doors  are 
opened. 


An  artist’s  impression  of  the  Syncrolift  lock  is  shown  in  Figure  2-7.  A 
schematic  diagram  of  Syncrolift  combined  with  flushing  of  salt  water  is  shown  in 
Figure  2-8.  Dimensions  of  Syncrolift  under  consideration  for  Panama  Canal  site 
are  shown  in  Figure  2-9.  This  information  and  additional  details  are  available  in 
an  article  published  in  the  Dock  and  Harbor  Authority,  January/February  1999. 


Evaluation  of  Mitigation  Options 

The  following  points  need  to  be  considered  while  evaluating  alternatives  for 
controlling  salinity  intrusion.  Combination  of  alternatives  may  have  to  be 
considered  to  achieve  the  intended  objectives. 

Construction  feasibility 

Low  initial  cost 

Low  operational  cost 

Long  life 

Ease  of  operation 

Low  vessel  throughput  time 

Low  increase  in  lake  water  salinity 

Low  fresh  water  consumption 
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Figure  2-1 .  Existing  locks  layout  and  dimesions. 
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Figure  2-2.  Panamax-size  vessel  and  required  lock  dimensions 


Figure  2-3.  Ship  transits  1995 


Figure  2-5.  Ship  transits  1997 


Figure  2-6.  Ship  transits  1998 


Figure  2-7.  Artist’s  impression  of  a  syncroTrft  lock.  (Source: 
Syncrolift,  Inc.  Reprinted  with  permission.) 
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Salt  Loading  and 
Freshwater  Usage 


This  section  describes  the  process  of  salinity  intrusion  and  some  straight¬ 
forward  analytical  methods  for  calculating  the  amount  of  salt  transported  into  the 
canal  through  the  locks  and  the  amount  of  fresh  water  lost  to  the  sea  through  the 
lockage  process.  The  analytical  methods  are  used  to  back-calculate  the 
maximum  salinity  that  can  occur  in  the  Syncrolift  lock  without  increasing  salinity 
intrusion  into  Gatlin  Lake  and  also  to  evaluate  one  option  for  conserving  fresh 
water  within  the  canal  -  recirculating  lift  water  by  collecting  it  in  holding  ponds 
and  pumping  it  back  for  reuse. 


Salinity  Intrusion 

Seawater  contains  about  35  grams  of  dissolved  salts  per  kg,  with  the  primary 
salt  constituents  being  chloride  and  sodium  ions.  This  concentration  is 
commonly  referred  to  as  “salinity,”  and  was  traditionally  defined  in  terms  of  the 
mass  concentration  of  halides  (principally  chlorine  forms)  in  the  water  and 
expressed  in  parts  per  thousand  (ppt).  The  standard  was  “Standard  Seawater”  or 
‘‘Normal  Water”  taken  from  the  North  Sea  and  prepared  for  scientific  use  and 
instrument  calibration  by  the  Hydrographic  Laboratories  in  Copenhagen, 
Denmark,  and  later  by  the  Institute  for  Oceanographic  Services,  Wormly, 
England.  The  definition  of  salinity  has  changed  slightly  several  times  in  response 
to  practical  demands  and  improved  instrumentation. 


In  1978  an  international  panel  sponsored  by  the  United  Nations  Educational, 
Scientific,  and  Cultural  Organization  (UNESCO)  recommended  that  the  old 
salinity  definitions  be  replaced  with  a  “Practical  Salinity  Scale”  based  on 
conductivity  and  employing  no  units,  replacing  the  ppt  units  with  the  sign  *10'3. 
The  new  scale  was  defined  to  yield  values  close  to  the  old  scales,  so  that  35  ppt 
was  about  35*10°.  That  notation  never  caught  on,  but  recently  the  units  of 
“practical  salinity  units”  (psu)  have  come  into  more  common  use.  Nevertheless, 
in  this  report  the  traditional  units  of  ppt  are  used  to  describe  salinities  defined 
under  the  new,  conductivity-ratio-based  definition. 

The  presence  of  dissolved  salts  increases  the  density  of  water,  such  that  at 
35  ppt,  water  is  about  2  percent  more  dense  than  at  0  ppt.  This  seemingly  small 
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difference  produces  density  currents.  Given  still  salt  water  on  one  side  of  a 
barrier  and  fresh  water  on  the  other  side,  hydrostatic  pressure  is  equal  at  the 
water  surface,  but  the  pressure  on  the  salt  water  side  increases  more  rapidly  with 
depth.  This  pressure  difference  produces  a  net  force  toward  the  fresh  water  side, 
which  will  elevate  the  fresh  water  side  if  the  barrier  is  removed.  The  combina¬ 
tion  of  a  raised  fresh  water  level  and  increased  salt  water  pressure  with  depth 
creates  a  flow  of  salt  water  under  the  fresh  and  of  fresh  water  over  the  salt  (see 
Part  5).  The  speed  of  the  density  underflow  will  be  proportional  to  the  square 
root  of  the  density  difference  (Keulegan,  1969). 

As  noted  previously,  density  flows  will  occur  from  lock  chambers  containing 
water  that  is  either  saltier  or  fresher  than  the  water  into  which  the  gates  are 
opened.  When  chamber  containing  fresh  water  is  opened  to  a  salt  water  basin, 
the  salt  water  will  flow  into  the  chamber  over  the  lower  portion  of  the  water 
column  and  fresh  water  will  flow  out  over  the  upper  portion  of  the  column.  If 
turbulence  is  present,  the  two  layers  will  mix,  at  least  partly. 


Calculating  Salinity  and  Exchanges  in  the  Locks 

The  amount  of  salt  intruding  into  Gatun  Lake  from  lockages  can  be 
calculated  by  analytical  or  numerical  means.  In  this  section  a  set  of  simple  mass 
balance  equations  are  used  to  estimate  salt  loading  and  freshwater  usage  by  the 
locks.  For  convenient  reference,  the  six  locks  in  the  system  are  given  the 
following  notation: 

LI :  Lowest  lock  on  the  Pacific  side 

L2:  Lock  between  Miraflores  Lake  and  LI . 

L3:  Lock  between  Miraflores  Lake  and  Gatun  Lake.  (Pedro  Miguel) 

L4:  Highest  lock  directly  connected  to  Gatun  Lake. 

L5:  Middle  lock  on  Atlantic  side. 

L6:  Lowest  lock  connected  to  Atlantic  Ocean. 

A  schematic  longitudinal  section  of  all  the  locks  is  given  in  Figure  A-l  of 
Appendix  A  and  in  Table  2-1,  which  provide  dimensions,  lifts,  and  relevant 
levels  of  all  the  locks. 

In  these  calculations,  the  following  assumptions  are  made: 


a.  Tidal  variation  in  sea  level  can  be  neglected.  All  computations  are  made 
for  mean  sea  level  (msl). 

b.  Mean  sea  level  can  be  assumed  to  be  roughly  equivalent  to  the  Canal’s 
Precise  Level  Datum  (PLD)  without  biasing  the  calculations. 

c.  Number  of  lockages  is  assumed  to  be  the  same  as  the  number  of  ship 
transits. 

d.  An  average  blockage  ratio  of  0.8  is  assumed  for  calculating  submerged 
volume  of  ship  from  the  length,  beam  and  draft  dimensions. 
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e.  Ocean  salinity  is  assumed  to  be  35  parts  per  thousand.  Observations  (see 
Part  2)  indicate  that  the  local  salinity  near  the  lowest  gate  is  much  less 
than  ocean  salinity  because  fresh  water  from  the  lake  is  added  to  the  sea 
for  each  ship  transit  and  a  transition  zone  exists  between  the  lowest  locks 
and  the  open  sea.  Based  on  the  limited  observations  it  is  assumed  that 
salinity  just  outside  the  lowest  lock  is  10  ppt  on  the  Pacific  side  and  8  ppt 
on  the  Atlantic  side. 

f.  Volume  of  water  in  Miraflores  Lake  is  about  1 9  million  cu  m 
(670  million  eft). 

g.  Net  transfer  of  fresh  water  from  the  lake  to  the  sea  is  about  0.2  million 
cu  m  (52  million  gallons)  per  lockage,  about  0.1  million  cu  m  (26  million 
gallons)  directly  to  fill  the  locks  and  about  0.1  million  cu  m  (26  million 
gallons)  pushed  into  the  highest  lock  by  density  flow. 

The  following  generic  lock  exchange  equations  can  be  used  to  characterize 
lock  water  and  salt  exchanges.  Given  the  schematic  lock  shown  in  Figure  3-1, 
the  left-hand-side  (LHS)  can  be  another  lock  or  the  sea  and  the  right-hand-side 
(RHS)  can  be  another  lock  or  Miraflores  Lake,  Gaillard  Cut,  or  Gatun  Lake.  The 
initial  conditions  (time  =  0)  for  lock  volume  and  salinity  are  given  by: 


5(0)  =  So  (r 

V(0)  =  Vt-Vl 

where  VT  =  total  lock  volume  when  filled  to  full  lift  height  and  VL  =  volume 
required  to  fill  lock  from  minimum  fill  condition  (min  in  Figure  3-1)  to  full  lift 
height  (max).  From  this  initial  condition,  a  one-way  upbound  ship  transit  from 
left  to  right  occurs  in  four  steps  as  given  below.  It  is  assumed  that  dissolved  salt 
in  the  lock  chamber  is  completely  mixed1  at  the  end  of  each  step. 

Step  1 :  Open  LHS  gates,  move  ship  into  lock,  close  gates. 

The  displacement  of  the  entering  ship  forces  an  equivalent  net  volume  of 
water  out  of  the  lock  to  the  LHS  and  density  currents  force  opposite  but  equal 
exchanges  of  saltier  water  from  the  LHS  for  fresher  water  in  the  lock.  The 
volume  of  water  exchanged  by  density  flows,  Vlhs,  can  be  expressed  as  a  fraction 
of  the  water  volume  in  the  lock: 

Vlhs  =  Elhs(Vt  -  Vs)  (2) 

where  ELhs  -  an  exchange  coefficient  between  0  and  1 ,  Vs  =  displacement 
volume  of  the  ship,  and  the  net  volume  (VT  -  Vs)  is  the  lesser  of  the  two 


The  complete  mixing  assumption  is  reasonable  overall,  since:  a)  lower  salinity  (thus  usually 
lower  density)  filling  water  enters  the  chamber  from  the  bottom  and  rises  through  higher  salinity 
water  with  considerable  turbulence;  and  b)  during  movement  from  one  chamber  to  the  next,  ship 
propwash  generates  substantial  mixing. 
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connected  chamber  volumes.  At  the  end  of  Step  1,  the  volume  and  salinity  of  the 
water  inside  the  lock  chamber  are  given  by: 

F(l)  =  V(0)  -  Vs 
=  Vr  =  Vl 

5( 0)F(0)  -  S(0)Vlhs  +  Slhs  Vuk  -  5(0)F.  (3) 

S\  '  ~~  T2V1  \ 


5(0)(F(0)  -  Vlhs  -  Vl)  +  SlhsVlhs 

_ 


where  Slhs  =  salinity  of  the  water  outside  the  lock  on  the  left-hand-side. 

Step  2:  Fill  lock  chamber 

Applying  conservation  of  mass  for  water  and  salt  once  more  yields  equations 
for  volume  and  salinity  within  the  chamber  at  the  end  of  Step  2: 


V(2)  =  Vt-Vs 


5(2)  = 


S{\)V(\)  +  SlVl 
V(2) 


(4) 


where  SL  =  salinity  of  the  lift  water,  which  can  be  the  salinity  of  the  water  outside 
the  right-hand-side  of  the  lock  (either  the  next  higher  lock  or  lake)  or,  if  a  holding 
pond  is  employed,  the  salinity  of  water  in  the  holding  pond. 

Step  3:  Open  right-hand-side  (RHS)  gates,  move  the  ship  out  of  the  lock,  and 
close  the  gates 

During  the  step  there  is  densimetric  exchange  with  the  RHS  and  water  from 
the  RHS  flows  into  the  lock  to  replace  the  ship  displacement.  The  densimetric 
exchange  volume  is  given  by: 

Vms  =  Erhs(Vt  —  Vs) 


At  the  end  of  the  step  the  volume  and  salinity  within  the  chamber  are: 


V(3)  =  Vt 


5(3) 


S(2)(V(2)  -  Vrhs)  +  Srhs(Vrhs  +  Vs) 

m 


(6) 


where  Srhs  =  salinity  of  the  water  outside  the  RHS  gates. 
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Step  4:  Drain  lift  volume  from  the  lock 


At  the  end  of  step  4  the  volume  and  salinity  are  given  by: 

V(4)  =  Vt-Vl 
5(4)  =  5(3) 


and  the  lock  is  in  the  same  fill  configuration  as  the  initial  condition  so  that 
another  cycle  can  begin. 

By  substitution  of  known  values  and  algebraic  manipulation,  these  equations 
can  be  used  to  derive  an  expression  for  the  salinity  in  the  lock  chamber  at  the  end 
of  N  one-way,  upbound  lockages: 


Sn  =  S(0)  +  N(S(4)-S(0)) 

=  So+n[  So  -  Vt —  Vl-Vs-Vrhs  )(  Vt-Vs-Vl-Vlhs  ) 

1  [  Vt  (  Vt-Vs  )  1 

+  &  Vl  (Vt -Vs -Vms)(  Vt-Vs-Vl-Vlhs)  " 

Vt  (Vt-Vs  )(  Vt-Vl) 

,  0  \  Vkhs  +  Vs  1  _  Vlhs  (  Vt-Vs-Vl-Vlhs  ) 

+  Srhs  -  +  Slhs  — - - - - i 

Vt  J  1_  Vr(Vt-Vs) 


A  similar  process  can  be  used  to  develop  expressions  for  the  salinity  with 
one-way  downbound  lockages  and  two-way  lockages. 

Equation  8  has  been  coded  in  spreadsheet  format  to  facilitate  evaluation  of 
some  salt  mitigation  alternatives  for  specific  locks,  and  results  for  existing 
conditions  are  shown  in  Tables  3-1  and  3-2  for  the  Miraflores  and  Gatun  Locks, 
respectively.  The  assumed  initial  condition  was  zero  salinity  in  the  locks  and 
Miraflores  Lake,  and  10  ppt  just  outside  the  lowest  lock  on  the  Pacific  side  and 
8  ppt  on  the  Atlantic  side  (see  Part  2).  The  exchange  coefficients,  Ews  and  Ems, 
can  be  determined  empirically  by  fitting  Equation  8  to  observed  salinities  in  the  ’ 
Panama  Canal  Locks.  Given  the  limited  observations  of  Part  2,  they  were 
estimated  to  be  about  0.5  (Keulegan,  1957),  and  Tables  3-1  and  3-2  show  them  to 
yield  results  that  are  reasonably  consistent  with  the  the  observed  salinities 
(Part  2).  After  about  20  lockages  the  lock  chamber  salinities  have  reached 
equilibrium  values  ranging  from  0.1  ppt  in  L4  to  3  ppt  in  LI.  Tuning  of  these 
coefficients  will  generate  lock  salinities  in  even  closer  agreement  with  the 
observations,  but  the  limited  number  of  observations  and  the  simplicity  of  the 
method  do  not  justify  it.  Appendix  A  employs  the  basic  equations  with  the 
exchange  coefficients  transformed  into  transfer  coefficients  in  order  to  precisely 
match  observed  lock  salinities  and  then  compute  salt  loadings  and  water  flows 
through  each  set  of  locks.  Those  calculations  show  the  salt  loading  to  be 
extremely  sensitive  to  the  calculated  in-lock  salinities.  For  example,  varying  the 
ambient  salinity  of  Miraflores  Lake  from  1.1  ppt  to  1.4  ppt  causes  the  net  salt 
input  to  the  lake  to  vary  from  +17,000  lbs  per  transit  to  -44,000  lbs  per  transit.  A 
near-equilibrium  salinity  occurs  if  the  lake  has  a  salinity  of  1.18  ppt. 
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Table  3-1 

Pacific  Locks  Salinity 


PARAMETER 


MIRAFLORES 


TOTAL  FILLED  VOL  (VT)  cf 


LIFT  VOL  (VL)cf 


SHIP  VOL  (VS)  cf 


EXCHANGABLE  VOL,  cf 


VT-VL 


vr-vs 


VT-VL-VS 


EXCHANGE  COEFS 
(ELHS  &  ERHS) _ _ 


EXCHANGE  VOL 
(VLHS.&VRHS)  _ 


TERM  1  |S0* 


TERM  2 


TERM  3  SLHS* 


TERM  4  SRHS* 


9.02E+06|  |  8.76E+06]  |  8.70E+08 


3.01  E+06  3.30E+06] _ |  3.65E+06 


1.14E+06  I  1.14E+06| _ |  1.14E+06 


6.01  E+06 1  5.46E+06]  5.46E+06|  5.46E+06|  8.66E+08|  5.06E+06| 


6.01  E+06  5.46E+06] _ |  8.66E+08 


7.88E+06)  |  7. 62 E+06 1 _ |  8.69E+08 


4.87E+06  4. 32 E+06] _ |  8.65E+08 


oTsl  05m  050  0.50  0.50  0.50 


8.71  E+06 


3.65E+06 


1.14E+06 


5.06E+06|  5.06E+06 


5.06E+06 


7.57E+06 


3.92  E+06 


0.50  0.50 


3.01  E+06 1  2.73E+06I  2.73E+06I  2.73E+06|  4.33E+08I  2.53E+06|  2.53E+06  2.53E+06 


5.63E-02 


9.07E'02 


2.18E-01 


4.29E-01 


3.79E-02 


7.86E-02 


2.00E-01 


4.42E-01 


4.93E-01 


4.17E-03 


4.96E-01 


4.22E-03 
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Table  3-2 

Atlantic  Locks  Salinity 


PARAMETER 


TOTAL  FILLED  VOL  (VT)  cf 


LIFT  VOLUME  (VL)  cf 


SHIP  VOLUME  (VS)  cf 


EXCHANGABLE  VOL,  cf 


VT-VS 


VT-VL-VS 


EXCHANGE  COEFS 
(ELHS&ERHS) 


EXCHANGE  VOL 
(VLHS&VRHS) 


TERM  1 


TERM  2  |SL* 


TERM  3  ISLHS* 


TERM  4  SRHS 


LOCKAGES 


Z50E+06]  2.50E+06 


3.08E-02 


LHS  SALINITY 


LOCK  6  SALINITY 

ppt 


7.58E-02 


1.96E-01 


LOCK  5  SALINITY 
ppt 


2.42E+06 


7.48E-02 


1.91  E-01 


4.19E-01 


LOCK  4  SALINITY 
ppt 


GATUN 
SALINITY  ppt 
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Tables  3-3  and  3-4  list  the  results  of  the  salt  loading  calculations. 


Table  3-3 

Estimates  of  Salinities  Resulting  from  Ships  Transiting  from  Pacific 
Entrance 

Location  and  stage 

Salinity  (ppt) 

Transfer  Coefficient  (%) 

Pacific  Ocean 

35.0 

100.0 

Pacific  Ocean 
just  outside  LI 

10.0 

28.6 

Inside  Lock  LI 
after  ship  enters 

8.0 

22.9 

Inside  Lock  LI 
after  full  lift 

5.5 

15.8 

Inside  Lock  L2 
after  ship  enters 

1.5 

4.3 

Inside  Lock  L2 
after  full  lift 

1.3 

3.7 

Inside  Lock  L2 
after  ship  leaves 

1.0 

2.9 

Miraflores  Lake 

1.0 

2.9 

Inside  Lock  L3 
after  ship  enters 

0.3 

0.7 

Inside  Lock  L3 
after  full  lift 

0.1 

0.4 

Table  3-4 

Estimate  of  salt  transfer  coefficients  from  Atlantic  Entrance 

Location  and  stage 

Salinity  (ppt) 

Transfer  coefficient  %) 

Atlantic  Ocean 

35.0 

100.0 

Atlantic  Ocean 
iust  outside  L6 

10.0 

28.6  ! 

Inside  Lock  L6 
after  ship  enters 

8.0 

22.8 

Inside  Lock  L6 
after  full  lift 

5.4 

15.6 

Inside  Lock  L5 
after  ship  enters 

1.5 

4.3 

Inside  Lock  L5 
after  full  lift 

0.8 

2.3 

Inside  Lock  L4 
after  ship  enters 

0.2 

0.4 

Inside  Lock  L4 
after  full  lift 

0.1 

0.2 
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Zero-Dimensional  Model  of  Gatun  Lake  Salinity 

Using  the  basic  approach  given  in  the  above  exchange  equations,  a  zero- 
order  numerical  model1  was  developed  to  describe  the  salinity  accumulation  in 
Gatun  Lake.  The  model  computes  the  salinity  of  Gatun  Lake  based  upon  the 
inflow  of  fresh  water  from  Madden  Lake,  local  tributaries,  exchange  with  the 
Pedro  Miguel  and  Gatun  locks,  and  the  influx  of  salt  from  the  proposed 
Syncrolift  Locks  on  both  ends  of  the  canal.  Equation  9  shows  the  mass  balance 
equation  used  for  the  model. 


VCL  *  Sot  +  (  VMadden  +  Vtrib  )  *  Sxiadden 

+  (F'z,3-Ls)*Sl  3  *EXl  3  *  N  +  (  V  'l  a  -Vs  )  *Sl  a*EXla  *  Nl 
v  +2  *  ((  Vsync  -  v££  )*Ssync*  EXsync  *  Nsync  ) 

VCL  +  VMadden  +  Vtrib 

+  ((Vl  3-Vs)*EXl  3  +  (V  'l4-Vs  )*EXl4  )*Nl 
,+2  *  ((Vsync  -  V  ;%)*  EXsync*  Nsync)  ■ 


where  Sgl+i  is  salinity  in  Gatun  Lake  at  the  end  of  the  time  step,  SgjJs  the 
salinity  in  Gatun  Lake  at  the  beginning  of  the  time  step  and  VGL  is  the  volume  of 
Gatun  Lake;  VMadden  and  V^b  are  volumetric  inflows  into  Gatun  Lake  from  Lake 
Madden  and  local  tributaries,  SMadden  is  the  salinity  of  the  inflows  from  Madden 
Lake  and  the  local  tributaries;  Vs  is  ship  volume,  V'L3}  Su,  and  EXU  are  the  lock 
volume,  salinity,  and  average  exchange  ratio  for  lock  3  (Pedro  Miguel  Locks); 
V'u,  Su,  and  EXU  are  the  lock  volume,  salinity,  and  exchange  ratio  for  L4 
(Uppermost  lock  at  Gatun  Locks);  NL  is  the  number  of  lockages  per  time 
interval  at  the  Gatun  or  Pedro  Miguel  Locks;  and  V^c>  vgg ,  S^c,  EX^,  and 

E'svnc  are  the  lock  volume,  ship  volume,  salinity,  exchange  ratio,  and  number  of 
lockages  per  time  interval  for  the  Syncrolift  facilities. 

The  observed  salinity  in  and  around  the  LI  (Miraflores  Locks)  provides  a 
means  of  estimating  the  in-lock  salinity  of  the  Pedro  Miguel  Locks.  The  salinity 
of  the  brackish  water  seaward  of  LI  was  measured  at  8-10  ppt,  a  reduction  by  a 
factor  of  nearly  4  compared  to  seawater.  It  seems  likely  that  the  Miraflores 
Lake’s  salinities  would  demonstrate  similar  tendencies  in  the  vicinity  of  Pedro 
Miguel  Locks.  Thus,  the  salinity  at  the  entrance  to  L3  would  be  less  than  0.7  ppt 
by  a  factor  of  3-4,  estimated  at  0.200  ppt.  The  dilution  that  occurs  during  filling 
of  L3  with  fresh  water  will  further  reduce  the  salinity  exposed  to  Gatun  Lake  to 
approximately  0. 1 00  ppt. 


Using  this  loading  rate  and  the  zero-order  model,  a  steady-state  salinity  can 
be  estimated  for  Gatun  Lake  at  about  0.030  ppt  (Figure  3-2),  which  is  well  below 
the  maximum  salinity  for  drinking  water  of  0.250  ppt  and  at  an  appropriate  level. 


EXCEL  spreadsheet  model  entitled  Lake  Gatun  Zero-Order  Salinitv  Intrusion  Model  vl  1  Sep 
1999.  '  V 
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since  there  are  no  reported  salinity  problems  within  Gatun  Lake.  The  zero-order 
model  also  shows  the  time  scale  of  salinity  equilibrium  into  Gatun  Lake  to  be  on 
the  order  of  6-10  years.  Table  3-5  gives  the  conditions  used  in  this  calculation 
and  for  other  examples  presented  in  subsequent  figures. 


Table  3-5 

Data  Used  in  Zero-Order  Model  Analysis 

Location/Process 

Variable 

Case  1  |  Case  2  |  Case  3  |  Case  4 

Gatun  Lake 

Time  interval 

1  yr 

Initial  salinity,  ppt 

0 

Volume,  cf 

1.82E+1 1 

Madden  Lake  and 
Tributary  Inflow 
Averaqe 

Annual  spill  volume,  cf 

5.73E+09 

Annual  power  volume,  cf 

6.66E+10 

Annua!  Tributary  Volumes,  cf 

4.53E+10  ! 

Salinity,  ppt 

0 

Exchange  during 
lift  operations 

Pedro  Miguel 

Locks 

Annua!  lock  volume,  cf 

5.12E+10 

Exchange  coefficient 

0.5  1 

Lock  salinity,  ppt 

0.1 

Number  of  lockages/d  ay 

38 

Exchange  during 
lift  operations 

Gatun  Locks 

Annual  lock  volume,  cf 

5.45E+10 

Exchange  coefficient 

0.5 

Lock  salinity,  ppt 

0.1 

Number  of  lockages/dav 

38 

Exchange  during 
lift  operations 
Miraflores  syncrolift 

Annual  lock  volume,  cf 

0 

8.74E+09 

Exchange  coefficient 

0 

0.5 

Lock  salinity,  ppt 

0 

"35  TT5  f03 

Number  of  lockages/day 

20 _ 

Exchange  during 
lift  operations 

Gatun  Syncrolift 

Annua!  lock  volume,  cf 

0 

8.74E+09 

Exchanqe  coefficient 

0.5 

Lock  salinity 

"o  ns  ro  las 

Number  of  lockages 

20 

The  Syncrolift  facilities  were  added  to  the  loading  equation  with  an  exchange 
coefficient  of  50  percent,  similar  to  the  conventional  locks.  Its  volume  was  based 
on  dimensions  of  243.8  m  (800  ft)  long  by  26.9  m  (88.5  ft)  wide  by  9.8  m 
(32  ft)  deep,  which  will  contain  a  ship  with  a  displacement  of  30,242  cu  m 
(1.068  million  cubic  feet).  With  an  in-lock  salinity  of  35.0  ppt  (seawater),  the 
steady  state  condition  within  Gatun  Lake  would  be  about  1 .7  ppt  (Figure  3-3)  for 
20  lockages  per  day.  Even  for  a  much  lower  salinity  loading  of  10  ppt  (brackish 
water  outside  of  Miraflores  Locks),  the  steady  state  concentration  is  estimated  at 
0.500  ppt,  which  is  double  the  drinking  water  maximum  of  0.250  ppt.  These 
results  clearly  indicate  that,  without  mitigating  technology,  salinity  in  Gatun 
Lake  could  become  unacceptably  high  once  additional  locks  are  put  into 
operation. 
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Varying  the  salinity  of  the  Syncrolift  lock  provides  insight  into  the  potential 
impacts  on  Gatun  Lake  salinity.  Without  exceeding  the  drinking  water  maximum 
in  Gatun  Lake,  the  allowable  salinity  in  the  Syncrolift  Lock  is  approximately 
4.5  ppt  (Figure  3-4).  However,  if  an  allowable  increase  in  Gatun  Lake  salinity  of 
0.010  ppt  is  adopted  (from  0.030  to  0.040  ppt),  then  the  salinity  within  the 
Syncrolift  Lock  could  not  exceed  0.300  ppt  (Figure  3-5). 


Recirculation  of  Lift  Water  in  Holding  Ponds 

In  this  scenario,  depicted  schematically  in  Figure  3-6,  the  lockage  process 
consists  of  the  steps  listed  above  for  a  one-way  transit  upbound,  and  when  the 
uppermost  lock  (L3  or  L4)  is  filled,  part  of  the  lift  volume  is  drawn  from  an 
upper  holding  pond  and  part  from  Gatun  Lake.  When  the  lowest  lock  (LI  or  L6) 
is  drained  to  accept  the  next  ship  at  sea  level,  the  lift  volume  is  drained  into  the 
lower  holding  pond.  Water  from  the  lower  holding  pond  is  pumped  to  the  upper 
pond  for  reuse. 


In  a  recirculation  operation,  the  salinity  of  the  water  used  to  fill  the 
uppermost  lock  will  be  the  salinity  of  the  upper  holding  pond,  which  at  the  end  of 
one  lockage  is  given  by: 


Sl  =  Suhp  (  1  ) 


Suhp  (  0  )  Vuhp  (  0  )  -  Suhp  (  0  )Vl  +  VpSlhp  (  0  ) 
Vuhp  (  1 ) 


(10) 


where  the  numbers  in  parentheses  indicate  the  lockage  number,  Suhp  is  salinity  of 
the  upper  pond,  SlhP  is  salinity  of  the  lower  pond,  VP  is  the  water  volume  pumped 
from  the  lower  pond  to  the  upper  pond  during  one  lockage,  and  VUHP  is  volume 
of  the  upper  pond,  given  by: 


Vuhp(  1)  =  Vuhp  (  0  )  -  Vl  +  Vp 


(11) 


The  salinity  of  the  lower  holding  pond  after  one  lockage  will  be: 


Slhp  (  1 ) 


Slhp  {0)Vlhp(Q)  +  Sll(0)Vl-Slhp  (0)Vp 
Vlhp  ( 1 ) 


(12) 


where  SLL  is  the  salinity  of  the  lowest  lock  chamber  at  the  end  of  a  lockage  and 
Vlhp  is  volume  of  the  lower  holding  pond,  given  by: 


Vlhp  (  1 )  =  Vlhp  (0)  +  Vl-Vp 


(13) 


During  each  one-way  lockage,  the  volume  of  fresh  water  consumed  by  the 
existing  procedure  over  N  lockages  can  be  approximated  by: 
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(14) 


N 

Vf(N)  =  [  Vl  ( 1  —  F  )  +  Vkhs  ,hl  +  Vs  ] 

n- 1 


where  Vrhs.ul  is  the  right-hand-side  exchange  volume  for  the  highest  lock.  With 
the  holding  pond  recirculation  procedure,  the  freshwater  consumption  will  be: 

N 

Vf  ,p  (N)  =  Vuhp  +  ^[Vl(\-F)  +  Vrhs  ,hl  +  Vs]  (15) 

n—\ 

where  F  is  the  fraction  of  the  lift  volume  drawn  from  the  upper  holding  pond. 

More  fresh  water  will  be  consumed  by  two-way  lockages  than  by  one-way, 
since  an  extra  lift  volume  is  used  to  fill  one  lock  chamber  that  cannot  be  filled 
from  the  next  higher  lock,  but  the  ratio  of  the  two  volumes  (Equation  15  divided 
by  14)  for  one-way  lockage  can  be  used  to  approximate  the  overall  savings  in 
fresh  water  consumption. 

Equation  8  and  Equations  10-15,  coded  in  a  spreadsheet1  and  combined  with 
the  following  assumptions,  provide  a  tool  for  examining  the  salinity  and  fresh 
water  usage  effects  of  this  operation.  The  additional  assumptions  are: 

a.  Each  pond  is  well-mixed  at  the  end  of  each  lockage. 

b.  The  salinity  just  seaward  of  the  lowest  locks  is  proportional  to  the  salinity 

of  the  water  in  the  lowest  lock  chamber,  but  cannot  rise  above  sea 
salinity  of  35  ppt. 

c.  The  small  differences  in  individual  lock  lift  volumes  is  neglected. 

Assumption  b  above  is  based  on  the  observation  that  salinities  just  seaward 
of  the  lowest  locks  were  just  8  to  10  ppt  instead  of  35  ppt,  apparently  because  of 
dilution  with  freshwater  from  lock  and  spillway  operations.  If  recirculating  lift 
water  reduces  freshwater  consumption,  those  seaward  salinities  will  rise.  A 
verified  numerical  model  of  the  waterways  seaward  of  the  locks  can  predict  the 
rise  in  salinity,  but  this  screening-level  evaluation  is  limited  to  the  reasonable 
approximation  that  the  ratio  of  salinity  just  outside  the  lock  to  salinity  just  inside 
the  lock  is  constant. 

In  the  spreadsheet  several  parameters  are  highlighted  in  yellow  to  indicate 
that  they  can  be  changed  to  evaluate  different  procedures  and  assumptions  for  the 
holding  pond  recirculation  operation.  In  the  following,  six  cases  are  explicitly 
considered,  as  listed  in  Table  3-6  with  the  equilibrium  salinity  results  for  Pedro 
Miguel  Lock  (L3)  and  highest  Gatun  Lock  (L4). 

Figures  3-7  and  3-8  show  the  results  for  Case  PI  for  the  Gatun  Locks  and 
Miraflores  Locks,  respectively.  The  Gatun  Locks  reach  an  equilibrium  condition 


1  EXCEL  spreadsheet  titled  “HOLDINGPOND.XLS.” 
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Tabl< 

Case 

3-6 

s  for  Holding  Pond  Alternative 

Case 

Parameter 

Results1  1 

Initial  Holding 
Pond  Salinity 
ppt 

Fraction  of 

Lift  Water 
from  Pond,  % 

Ratio  of 

Volume  of 
Upper  Holding 
Pond,  Vuhp/Vl 

Limiting 
Salinity  in 
Highest 
Lock,  ppt 

L3 

Equilibrium 

Salinity 

ppt 

L4 

Equilibrium 

Salinity 

ppt 

Fresh  Water 
Usage 
%  of  normal 

Base 

- 

0 

~ 

— 

0.09 

0.1 

100 

PI 

0 

100 

2 

0 

0.28 

1.3 

40 

P2 

0 

50 

2 

0 

0.16 

0.8  i 

70  i 

P3 

0 

100 

4 

0 

0.22 

1.3 

40 

P4 

0 

100 

2 

0.3 

0.16 

0.3 

60 

P5 

8 

100 

2 

0.3 

0.22 

0.3  i 

1  40 

P6 

4 

50 

2 

r0.3 

0.16 

0.3 

70 

I  Note  that  results  are  approximate  and  should  not  be  considered  precise  to  two  significant  digits.  Extra  digits  ar 
I  highlight  differences  in  results. 

e  shown  only  to 

after  about  80  lockages,  with  the  holding  ponds  reaching  a  salinity  of  about 
10  ppt  and  L4  reaching  1.3  ppt.  Despite  an  initial  increase  in  freshwater  con¬ 
sumption  due  to  filling  the  ponds,  Case  PI  reduces  long  term  consumption  sub¬ 
stantially,  to  about  40  percent  of  that  used  by  the  normal  method.  This  amounts 
to  one  lift  volume  of  0.1  million  cu  m  (3.6  million  eft)  per  lockage  once  the 
initial  two  lift  volumes  used  to  fill  the  pond  are  saved.  L3  salinity  does  not  rise 
as  high  as  L4  in  the  Gatun  Locks,  since  the  volume  of  Miraflores  Lake  buffers 
the  salinity  increase  in  L3  (at  the  cost  of  increasing  salinity  in  Miraflores  Lake), 
and  it  is  not  at  equilibrium  after  100  lockages.  Extension  of  the  calculations  to 
600  lockages  (the  number  shown  to  be  needed  for  equilibrium  by  the  zero-order 
model)  shows  L3  to  reach  a  maximum  of  about  0.28  ppt  for  case  PI. 

Cases  P4-P6  add  a  maximum  salinity  criterion  to  the  holding  pond  procedure 
to  reduce  salinity  intrusion  to  Gatun  Lake,  in  that  once  the  salinity  in  the  upper¬ 
most  lock  exceeds  a  preset  limit,  both  holding  ponds  are  flushed  and  refilled  with 
fresh  water.  The  criterion  is  expressed  as  a  limiting  salinity  as  shown  in  column 
5  of  Table  3-6.  Figure  3-9  illustrates  the  salinity  responses  of  L4  for  a  limiting 
salinity  of  0.3  ppt.  Once  the  salinity  limit  is  reached  at  about  12  lockages  it 
cycles  more  often  as  the  seaward  limit  salinity  rises,  requiring  refilling  the  ponds 
after  7  cycles  the  second  time  and  with  increasing  frequency  until  they  must  be 
flushed  and  refilled  every  3  lockages  at  100  lockages.  Despite  this  increased 
frequency,  the  method  is  seen  to  still  conserve  fresh  water,  reducing  consumption 
by  about  40  percent. 
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Figure  3-1 .  Schematic  lock  notation 
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Figure  3-2.  Case  1:  Gatun  Lake  Salinity,  38  locks/day  on  conventional  locks,  no 
lockages/day  on  Syncrolift  Locks,  lock  exchange  ratio  of  50  percent, 
salinity  of  0.100  ppt  in  conventional  locks  and  0.0  ppt  in  Syncrolift 
Lock 
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Figure  3-3.  Case  2:  Gatun  Lake  Salinity,  38  locks/day  on  conventional  locks, 
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Figure  3-4.  Case  3:  Gatun  Lake  Salinity,  38  locks/day  on  conventional  locks, 
20  lockages/day  on  Syncrolift  Locks,  lock  exchange  ratio  of 
50  percent,  salinity  of  0.100  ppt  in  conventional  locks  and  4.5  ppt 
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Figure  3-5.  Case  4:  Gatun  Lake  Salinity,  38  locks/day  on  conventional  locks, 
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50  percent,  salinity  of  0.100  ppt  in  conventional  locks  and  0.300  ppt 
in  Syncrolift  Lock 
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Figure  3-6.  Schematic  of  recirculating  lift  water  through  holding  ponds  and  pumping 


Figure  3-7.  Gatun  Locks  salinities  and  freshwater  consumption  with  holding 
pond  recirculation  for  Case  1:  ponds  holding  two  lift  volumes  with 
initial  salinity  of  0  ppt;  100  percent  of  lift  volume  from  holding  pond 
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Figure  3-9.  Gatun  Locks  salinities  and  freshwater  consumption  with  holding 

pond  recirculation  for  Case  P4:  ponds  holding  two  lift  volumes  with 
initial  salinity  of  0  ppt;  100  percent  of  lift  volume  from  holding  pond, 
and  limiting  salinity  of  0.3  ppt 
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Methods  of  Reducing 
Saltwater  Intrusion  in 
Navigation  Locks 


A  review  of  the  literature  was  performed  to  identify  methods  which  have 
been  shown  to  successfully  reduce  salt  water  intrusion  through  navigation  locks. 
Methods  and  designs  from  the  United  States,  the  Netherlands,  France,  and  China 
are  briefly  described  here.  This  section  draws  on  the  published  literature, 
Appendices  B  and  C,  and  an  unpublished  review  by  the  Panama  Canal 
Commission  (Kiamco  1999). 


Survey  of  Methods 

Flushing  canal  with  fresh  water 

A  repressive  method  of  combating  salinity  intrusion  through  locks  is  to  flush 
the  “collector”  canal  with  fresh  water  in  the  direction  of  the  sea.  However,  this 
method  is  not  very  effective  when  used  alone.  The  difference  in  density  allows  a 
comparatively  large  volume  of  salt  water  to  enter  the  canal  through  the  lock  and  a 
comparatively  large  volume  of  fresh  water  would  be  required  to  mix  with  it  and 
carry  it  away.  Tests  have  shown  that  the  dimensions  of  canals  very  greatly  affect 
the  degree  of  salt  intrusion.  A  small  increase  in  depth  causes  a  comparatively 
large  rise  in  the  salinity  of  the  water  in  a  canal. 


Pneumatic  barriers 

Pneumatic  barriers  are  curtains  of  air  bubbles  rising  from  perforated  pipes 
placed  on  the  bottom  of  a  channel  or  lock  chamber  perpendicular  to  the  axis  of 
the  channel  (See  Figure  4-1).  The  rising  bubbles  create  a  vertical  current  of 
water,  which  reduces  the  magnitude  of  any  exchange  currents  passing  over  it, 
provided  the  volume  of  air  is  adequate.  The  mechanics  of  flow  and  mixing  near  a 
pneumatic  barrier  are  described  by  McAnally  (1973). 

Pneumatic  barriers  reduce  salt  intrusion  that  is  caused  by  density  exchange 
currents.  This  is  important,  as  salt  intrusion  is  largely  due  to  density  exchange 
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currents.  The  efficiency  of  pneumatic  barriers  in  reducing  exchange  currents 
during  locking  operations  has  been  studied  extensively  in  the  laboratory  and  the 
field.  Field  studies  were  conducted  in  the  Ijmuiden  locks,  the  Netherlands. 
Salinity  intrusion  can  be  reduced  by  50  to  80  percent  depending  upon  the  amount 
of  air  used.  The  efficiency  decreases  with  increasing  water  depth  in  the  locks. 
Experience  shows  that  the  pneumatic  barriers  can  be  expensive  to  operate  and 
need  careful  operation  and  maintenance.  Further,  since  they  induce  intense 
mixing,  their  use  precludes  capture  techniques  such  as  selective  withdrawal  once 
mixing  has  begun. 


Water  barriers 

Theoretically,  a  vertical  current  like  that  produced  by  a  pneumatic  barrier  can 
also  be  induced  by  pumping  water,  thus  creating  a  water  barrier,  which  would 
also  reduce  the  magnitude  of  the  exchange  currents.  Initial  experiments  with  this 
method  have  been  carried  out  on  models. 


Selective  withdrawal  of  salt  water  during  intrusion 

The  salt  water  tending  to  enter  the  canal  when  the  inner  gates  of  a  lock  are 
opened  can  be  withdrawn  and  discharged  back  into  the  sea  during  the  exchange 
process  through  a  slot  in  the  bottom  of  the  lock  located  at  the  point  where  the 
lock  meets  the  adjoining  fresh  water  basin.  The  volume  of  fresh  water  wasted 
corresponds  to  the  quantity  of  water  withdrawn  through  the  slot.  If  withdrawal 
cannot  be  entirely  selective,  the  amount  of  fresh  water  wasted  during  each 
complete  locking  cycle  will  exceed  the  volume  of  the  lock  chamber.  This  system 
has  been  applied  at  the  Temeuzen-locks,  the  Netherlands.  Salt-water  intrusion  is 
prevented  at  the  cost  of  a  relatively  large  amount  of  fresh  water  due  to  limitations 
in  the  selective  withdrawal  of  salt  water. 


Selective  withdrawal  of  salt  water  after  intrusion 

The  intrusive  salt  water  is  collected  in  a  sump  of  adequate  volume  made  by 
deepening  part  of  the  canal,  selectively  withdrawn  from  it  and  discharged  in  a 
steady  stream  (see  Figure  4-2).  The  volume  of  fresh  water  wasted  corresponds  to 
the  quantity  of  water  withdrawn  from  the  sump.  Since  the  salt  water  is  diluted  as 
it  enters  and  while  it  remains  in  the  sump,  more  fresh  water  will  be  wasted  than 
when  using  the  method  described  in  item  4  above.  A  drawback  of  this  method  is 
that  not  all  the  intrusive  salt  water  is  trapped  in  the  sump;  some  of  it  passes  on 
into  the  canal  and  has  to  be  removed  by  flushing.  Large  volumes  of  fresh  water 
are  required  to  expel  even  small  quantities  of  salt  water.  The  basin  has  to  be 
large  enough  to  avoid  overtopping  by  salt  water  and  to  allow  selective  with¬ 
drawal  of  the  salt  water.  This  method  is  employed  at  the  Chittenden  Locks  in 
Seattle  Washington  as  described  in  a  subsequent  section. 
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Equalization  by  pumping  instead  of  gravity 

If  there  is  a  considerable  difference  between  average  canal  level  and  high 
water  on  the  seaward  side  of  the  lock,  salt  intrusion  can  be  greatly  reduced  by 
pumping  the  water  in  the  lock  chamber  back  into  the  sea  instead  of  allowing  it  to 
reach  the  canal.  Large-capacity  pumps  would  be  required  to  remove  the  large 
quantities  of  water  in  lock  chambers  in  a  short  time.  The  pumping  capacity 
required  could  be  reduced  by  storing  the  water  in  a  basin  between  the  lock  and 
the  sea  and  artificially  keeping  the  level  in  the  basin  lower  than  that  in  the  lock 
chamber.  This  advantage  of  the  method  is  that  there  is  more  time  for  pumping 
water  out  of  the  basin. 


Exchanging  while  gates  are  closed 

The  salt  water  can  be  withdrawn  from  the  lock  chamber  through  orifices  near 
the  bottom  and  replaced  by  letting  fresh  water  into  the  chamber  through  orifices 
near  the  surface  of  the  water.  Fresh  water  can  be  replaced  by  salt  water  in  like 
manner,  the  salt  water  being  admitted  near  the  bottom  and  the  fresh  water  being 
withdrawn  near  the  surface.  If  this  exchange  process  is  completed  before  the 
gates  are  opened,  waters  of  different  density  are  prevented  substantially  from 
mixing  thus  obviating  exchange  currents.  There  would  be  a  certain  amount  of 
mixing  while  the  exchange  is  being  effected;  the  vessels  in  the  lock  chamber 
would  also  cause  some  fresh  water  to  be  wasted  due  to  mixing.  The  fresh  water 
can  be  withdrawn  and  admitted  through  the  wall  of  the  chamber  or  through  the 
gates  at  the  fresh-water  end. 


Movable  sill 

One  of  the  measures  used  to  reduce  the  exchange  flows  is  to  have  a  movable 
sill,  which  is  adjusted  to  the  draft  of  the  ships  to  be  locked  (see  Figure  4-3).  This 
option  involves  risk  of  mechanical  damage,  which  may  make  the  option 
ineffective  or  may  block  navigation.  It  is  employed  at  the  Chittenden  Locks  in 
Seattle  Washington  as  described  in  a  subsequent  section. 


Simultaneous  exchange  of  salt  and  fresh  water 

Locks  of  this  type  have  perforated  false  floors  to  evacuate  or  to  introduce  salt 
water  and  have  openings  in  the  side  walls  for  the  exchange  of  fresh  water  while 
the  lock  gates  are  close.  The  system  has  been  applied  in  the  Mardyck  sluice  at 
Dunkirk  with  dimensions  of  145  x  12  x  4  m  (Ribes  and  Blanchet,  1965, 
Monadier,  1981)  and  the  Kreekrak  sluices,  the  Netherlands,  with  dimensions  of 
330  x  24  x  6  m  (Kolkman  and  Slagter,  1976).  The  saltwater  intrusion  can  be 
reduced  to  about  5  percent  of  the  lock  chamber  volume  with  a  freshwater  loss  of 
about  60  percent  of  the  lock  volume. 
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Canal  gate 

Another  option  for  the  lakeside  canal  is  to  have  two  or  more  deeper  basins 
and  a  gate  at  the  freshwater  end,  so  that  the  length  of  the  canal  is  available  for 
temporarily  trapping  salt  water.  By  itself,  such  a  design  will  offer  few  significant 
advantages,  but  in  combination  with  other  measures  it  can  be  used  to  retain 
brackish  water  until  such  time  as  flushing  can  be  accomplished.  For  example,  if 
a  movable  sill  and  selective  withdrawal  were  used,  but  sometimes  failed  to 
capture  all  the  salt  water,  a  long  canal  with  sills  at  intervals  would  trap  the  excess 
salt  water  for  an  extended  period,  then  it  could  be  flushed  (item  1)  when  an 
ample  supple  of  freshwater  became  available 


New  developments 

Various  new  lock  types  have  been  designed  and  tested  to  reduce  freshwater 
losses  (Van  der  Kuur,  1985,  Kolkman,  1986,  and  NN,  1986).  An  example  is  the 
lift  lock,  which  has  an  inner  lock  chamber  which  moves  vertically.  The  outer 
lock  chamber  is  twice  as  deep  as  normal.  The  inner  lock  chamber  is  always  filled 
with  salt  water.  Other  lock  types  have  been  studied,  such  as  locks  with  artificial 
reed  screens  to  decrease  exchange  flows  and  locks  with  horizontal  flexible 
membranes  to  prevent  mixing  between  salt  and  fresh  water  (van  der  Kuur,  1985). 
Freshwater  losses  with  this  type  of  lock  are  reduced  to  about  5  percent  of  the  lock 
chamber  volume,  whereas  saltwater  intrusion  is  only  about  5  percent  of  the  lock 
volume. 


Selecting  The  Method 

For  existing  locks 

A  method  that  might  be  considered  for  canals  linked  with  existing  locks  is  to 
flush  them  daily  with  copious  quantities  of  fresh  water.  The  quantity  of  water 
required  to  flush  wide,  deep  canals  is  considerable;  it  would  not  always  be 
possible  to  satisfy  the  maximum  permissible  salinity  requirements  in  certain 
places.  If  flushing  does  not  produce  the  desired  result,  a  stretch  of  canal  close  to 
the  lock  could  be  deepened  and  used  as  a  sump  in  which  the  salt  water  would 
collect.  A  pumping  station  operating  with  a  deep  suction  nozzle  would  withdraw 
the  salt  water  and  pump  it  into  the  sea.  If  necessary,  a  wall  could  be  built 
between  the  pumping  station  and  the  sump  with  orifices  near  the  bottom  of  the 
latter.  Alternatively,  the  salt  water  could  be  removed  by  gravity  if  the  water  level 
in  the  canal  were  either  temporarily  or  permanently  higher  than  that  of  the  sea. 

The  quantity  of  salt  water  intruding  into  an  existing  canal  owing  to 
exchange  flow  caused  by  differences  in  density  can  be  reduced  by  means  of  an 
air-bubble  curtain  placed  near  the  lock  gates.  This  option  could  be  combined 
with  the  sump  method  described  above  or  with  an  adjustable  barrier.  Whether  a 
single  method  or  a  combination  of  methods  is  adopted  depends  on  their  local 
feasibility. 
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For  projected  locks 


Installing  water  curtain  equipment  in  existing  locks  is  difficult  because  of  the 
size  of  the  pipes  required.  It  is  even  questionable  whether  this  method  would  be 
so  much  better  than  the  other  options  as  to  warrant  its  adoption  for  new  locks. 
The  qualitative  control  of  fresh  water  constantly  has  to  be  stepped  up  and  it  is 
essential  to  be  economical  with  fresh  water  for  many  reasons. 

There  are  four  methods  is  addition  to  the  three  already  described  for  existing 
locks  that  may  be  considered  for  new  locks,  viz. 

a.  Construction  of  a  sump  next  to  the  lock  on  the  landward  side  from  which 
the  intrusive  salt  water  can  be  selectively  withdrawn 

b.  Replace  salt  water  in  the  lock  chamber  by  fresh  water  while  the  lock 
gates  are  closed 

c.  Pump  the  water  required  for  equalization  into  the  sea. 

d.  Adjustable  barrier 


Navigation  Locks  in  the  Netherlands 

The  most  important  1 5  navigation  locks  in  the  Netherlands  and  the  methods 
employed  in  them  for  salinity  intrusion  control  are  described  below.  The  main 
features  of  options  are: 

a.  The  flushing  method  is  used  for  all  locks;  water  from  the  canal  is  flushed 
through  openings  placed  as  low  as  possible  in  the  lock  or  sluice 
structures. 

b.  The  velocity  of  the  current  passing  through  the  openings  may  exceed  a 
certain  magnitude  to  ensure  that  the  water  discharged  will  contain  the 
maximum  amount  of  salt. 

c.  Since  the  salt  water  pumped  from  the  polders  in  various  places  as  well  as 
that  coming  in  through  the  locks  reaches  the  canals,  the  fresh  water 
flowing  towards  the  locks  also  helps  to  overcome  the  salinization  of  the 
country’s  water  and  soil. 

d.  It  should  be  noted  that  it  would  always  be  necessary  to  flush  the  systems 
with  a  certain  quantity  of  fresh  water,  no  matter  how  effective  any  steps 
taken  to  reduce  the  waste  of  fresh  water  through  locks  may  be. 

e.  Air  curtain  equipment  has  already  been  installed  in  a  number  of  locks 
and  others  are  to  be  equipped  with  the  system  in  the  near  future. 

f  Supplementary  systems  are  being  considered  for  certain  locks. 
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g.  One  or  more  locks  at  seven  sites  in  the  Netherlands  are  equipped  with 
pneumatic  barriers,  viz.  at  Harlingen,  Den  Helder,  Ijmuiden,  Hartel 
canal/Brielse  Meer,  Haringvliet,  Volkerak  and  Temeuzen. 

Abraham  et  al.  (1973)  have  provided  a  brief  description  of  salinity  intrusion 
control  measures  at  several  docks  in  the  Netherlands.  A  section  of  that  report  is 
reproduced  verbatim  below. 

1  Delfzijl 

Salinization  of  the  canals  running  from  the  locks  is  countered  by  replacing 
the  water  in  the  canals  at  regular  intervals  by  rainwater  and  by  water  drawn 
from  Lake  IJssel.  The  locks  are  provided  with  apertures  located  well  below  the 
water  line  so  as  to  make  the  replacement  process  as  effective  as  possible. 

2  Lauwerszee 

Damming  off  the  Lauwerszee  has  very  greatly  reduced  the  salinity  of  the 
adjoining  parts  of  the  provinces  of  Friesland  and  Groningen.  Sluices  are 
incorporated  in  the  dam  through  which  excess  water  may  be  discharged. 
Provision  has  been  made  for  the  comparatively  small  lock  to  befitted  with  air- 
curtain  equipment ,  if  required. 

3  Harlingen 

Salt  intrusion  through  the  locks  is  countered  by  replacing  the  water  in  the 
canals  at  regular  intervals  by  rainwater  and  by  water  drawn  from  Lake  IJssel. 
Air-curtain  equipment  is  also  used  to  further  reduce  salinization. 

4  Zuyder  Zee  Barrier  Dam 

There  are  locks  in  the  dam  beside  the  two  sets  of  sluices  with  which  the  water 
in  Lake  IJssel  is  controlled.  Arrangements  are  being  made  to  equip  all  the  locks 
with  pneumatic  barriers. 

5  Den  Helder 

Salt  intrusion  through  the  three  locks  and  the  sea  dock  lock  is  countered  by 
replacing  the  water  in  the  canals  at  regular  intervals  by  rain  water  and  by  water 
drawn  from  Lake  IJssel.  Arrangements  are  being  made  to  replace  gravity 
discharge  by  pumping.  A  sump  is  to  be  built  near  the  pumping  station  in  which 
to  store  the  salt  water  from  the  locks  prior  to  its  being  pumped  into  the  sea.  The 
most  important  lock  is  equipped  with  a  pneumatic  barrier. 

6  Ijmuiden 

Salt  intrusion  through  the  North  Sea  Canal  is  countered  by  replacing  the 
water  at  regular  intervals  by  water  from  the  adjoining  regions  and  from  Lake 
IJssel.  The  four  locks  are  equipped  with  a  pneumatic  barrier.  Arrangements  are 
being  made  to  install  a  pumping  station  to  supplement  gravity  discharge  through 
a  sluice.  The  construction  of  a  deep  sump  to  collect  the  salt  water  is  also  being 
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contemplated  together  with  a  salt  baffle  to  maximize  the  salinity  of  the  water 
discharge  into  the  sea. 

7  New  Waterway-  New  Maas 

There  are  a  number  of  locks  along  the  New  Waterways  and  the  New  Maas, 
both  of  which  are  open  to  the  sea.  Salt  intrusion  in  the  canals  running  from  the 
locks  is  countered  by  replacing  the  water  by  rain  water  and  water  drawn  from 
the  River  Rhine.  The  quantity  of  river  water  required  could  be  somewhat 
reduced  by  installing  pneumatic  barriers  in  the  locks  but  there  are  as  yet  no 
plans  to  make  any  such  provisions. 

8  H artel  Canal  and  Brielse  Meer 

The  Hartel  canal  is  the  landward  shipping  link  between  Europoort  and  the 
Old  Maas.  There  are  pneumatic  barriers  in  the  two  locks  that  give  access  to  the 
Old  Maas  to  reduce  salinization  of  the  latter.  The  canal  is  also  flushed  with  river 
water ,  which  runs  out  through  the  locks,  which  are  also  equipped  with  pneumatic 
barriers. 

The  Brielse  Meer  is  a  fresh-water  lake  south  of  the  Hartel  canal;  a  lock  links 
it  with  the  Old  Maas.  The  lake  must  be  protected  as  well  as  possible  against 
salinization  if  the  salt  intrudes  beyond  the  lock  when  Rhine  is  low.  Consequently, 
a  great  deal  has  been  done  to  counter  salt  intrusion  at  this  point.  There  are  two 
pumps  in  two  culverts  by  means  of  which  water  can  be  pumped  out  of  the  lock 
chamber.  There  is  air-curtain  equipment  and  a  deep  sump  in  the  canal  next  to 
the  lock  to  which  the  pumps  can  be  connected  to  remove  any  salt  water  that  may 
get  through  the  lock. 

9  Haringvliet 

The  Haringvliet  lock  is  fitted  with  pneumatic  barriers.  There  are  deep 
culverts  in  the  sluice  next  to  the  lock  connected  with  a  sump  from  which  any  salt 
water  can  be  withdrawn. 

10  Volkerak 

There  are  two  locks  linking  the  rivers  with  Zeeuwse  Meer.  The  water  south 
of  the  dam  will  be  strongly  saline  until  the  Eastern  Schldt  dam  is  finally  closed  in 
1978;  meanwhile,  salt  intrusion  through  the  locks  is  being  reduced  by  means  of 
pneumatic  barriers. 

11  Canal  Across  Walcheren 

Salt  intrusion  through  the  locks  at  flushing  was  investigated  by  conducting 
model  experiments  to  discover  whether  a  canal  flushing  system  would  be 
effective  or  whether  the  locks  will  have  to  be  equipped  with  pneumatic  barriers  at 
both  ends. 
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12  Canal  Across  South  Beveland 

A  modem  system  of  countering  salt  intrusion  is  being  considered;  the  salt 
water  in  the  lock  chambers  would  be  replaced  by  fresh  before  opening  the  gates 
giving  access  to  the  fresh-water  lake.  The  fresh  water  used  would  serve  to  clean 
the  adjoining  canal. 

13  Scheldt-Rhine  Canal 

Two  locks  are  being  incorporated  in  the  canal  linking  Antwerp  with  the 
Rhine.  Salt  intrusion  will  be  countered  by  replacing  the  salt  water  in  the  lock 
chambers  by  fresh  water  flowing  in  through  apertures  in  the  walls  while  the 
gates  are  closed  As  in  the  new  system  described  in  13,  the  fresh  water  used 
would  serve  to  clean  the  adjoining  stretch  of  the  fresh-water  lake. 

14  Terneuzen 

There  is  one  old  and  two  new  locks  in  the  canal  between  Terneuzen  and 
Gent.  Pneumatic  barriers  reduce  salt  intrusion  through  the  new  inland- 
waterways  lock  A  sump  has  been  constructed  next  to  the  new  sea  lock  with  deep 
culverts  the  orifices  of  which  are  built  into  the  lock  to  carry  the  salt  water 
trapped  in  the  basin  towards  the  sea.  The  level  of  the  water  in  the  canal  is 
almost  always  above  that  of  the  high  tides,  so  discharge  is  by  gravity.  Moreover, 
salt  intrusion  is  further  reduced  by  means  of  pneumatic  barriers. 

15  Kreekrak 

The  Kreekark  lock  in  the  Netherlands  has  dimensions  of 320  x  24  x  5  m.  The 
lock  is  used  in  achieving  a  water  level  change  of  1.3  to  2.3  m.  It  is  navigable  for 
4-barge  push-tows  of 11,000  metric  tons  (10,000  tons),  with  the  design  annual 
traffic  capacity  of  66  million  metric  tons  (60  million  tons).  Fresh  water 
reservoirs  are  provided  on  both  sides  of  the  chamber.  A  salt-water  reservoir 
with  a  size  of  40  hectares  is  also  provided.  The  area  of  the  salt-water  reservoir 
is  26  times  that  of  the  chamber.  The  lock  is  operational  since  1986  and  is 
functioning  well.  The  Kreekrak  lock  is  an  improved  version  of  Mardyck  lock  in 
France.  It  also  uses  a  complete  exchange  system;  however,  the  fresh  water  loss 
is  reduced  to  about  15  to  30  percent.  The  reasons  for  improved  performance 
are:  a)  open  fresh  water  channel  used  instead  of  long  culverts,  b)  lift  gate  is 
provided  at  the  immediate  outlet  of  culverts,  which  can  be  adjusted  according  to 
tidal  level  changes,  thus  preventing  salt  water  intrusion,  and  c)  the  bottom  slab  is 
channel  shaped  and  water  being  discharged  through  outlets  on  the  sides  of 
channel  eliminate  vertical  currents,  which  produce  salt  and  fresh  water  mixing. 

16  Philip 

The  Philip  lock  has  dimensions  of 280  x  24  x  5  m  and  is  navigable  for  4- 
barge  pushtow  of 11,000  metric  tons  (10,000  tons).  The  lock  is  subjected  to  two- 
way  water  head  and  water  level  difference  of  1.7  m.  Philip  lock  made  the 
following  improvements  over  the  Kreekrak  lock:  1)  Since  outlets  in  channeled 
bottom  are  located  at  the  middle  of  the  channel  sides,  water  below  outlets  cannot 
drain  well,  so  bottom  slabs  of  T  and  inverted  T  shapes  are  used.  2)  Weir  gates  at 
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the  freshwater  outlets  are  not  closed  to  the  top;  thus  the  freshwater  between  the 
operating  valves  and  the  weir  gate  will  get  salted  with  a  loss  of  5%.  3)  When  salt 
water  surface  travels  towards  the  bottom  during  exchange ;  there  will  be  vertical 
currents  in  culverts  at  the  bottom  of  salt-water  body  causing  greater  flushing.  4) 
When  the  chamber  is  full  of  fresh  water  and  a  ship  enters ,  the  water  levels  in  the 
front  and  back  of  the  ship  are  different .  This  provides  opportunity  for  salt  water 
to  enter  the  chamber  and  fresh  water  to  enter  culverts. 


Hiram  Chittenden  Locks,  Seattle  WA,  USA 

Salt  water  intrusion  problems  in  USA  are  predominantly  related  to  estuaries 
and  inland  navigation.  A  summary  note  on  these  problems  is  given  in  Appendix 
B.  Other  problems  have  been  related  to  salt  water  intrusion  at  dams.  The 
primary  U.  S.  location  where  salt  intrusion  through  locks  has  been  addressed  is 
the  Hiram  Chittenden  Locks  on  the  Lake  Washington  Ship  Canal,  Seattle,  WA. 

Construction  of  these  locks  was  accomplished  during  191 1-1917  to  provide 
access  to  and  from  the  salt  water  regime  of  Puget  Sound  and  the  fresh  water  lakes 
of  Lake  Union  and  Lake  Washington  in  the  state  of  Washington,  USA.  The 
intrusion  of  salt  water  into  the  fresh  water  regime  is  controlled  through  a  number 
of  operational  and  structural  methods.  Foremost  are  the  salt  water  scavenging 
drain,  maintaining  salt  water  barrier  in  the  upright  position,  and  use  of  a  small 
lock  to  prevent  salt  water  intrusion  into  Lake  Washington.  During  critical  years, 
the  potential  for  intrusion  becomes  an  intensified  coordination  process  of  water  / 
lock  usage  with  municipal  water  supply  authorities,  environmental  agencies  and 
the  general  public  on  operational  uses  of  fresh  water  at  the  lock  structures. 
Appendix  C  describes  the  design  and  operations  employed  at  the  Chittenden 
Locks. 


Navigation  Locks  in  China 

The  Haihe  lock  in  China  has  made  the  following  improvements. 

a.  In  view  of  limited  land  area  available,  the  design  rejected  the  fresh  water 
reservoir  option.  Instead,  a  fresh  water  culvert  plan  is  adopted. 

b.  The  lock  design  adopts  grade  crossing  of  four-way  culverts  with  a  control 
valve  on  each  of  the  culverts  at  the  cross. 

c.  The  available  land  allows  to  build  a  reservoir  only  5  times  that  of  the  lock 
chamber  area.  Hence  a  pump  station  with  a  capacity  of  20  cubic  meters 
per  second  has  been  provided. 

d.  The  lock  has  less  flow  dividing  piers. 

e.  Froude  number  of  about  1 .0  is  achieved  at  the  point  of  water  exchange  at 

the  outlets.  ~ 
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Navigation  Locks  in  France 

The  Mardyck  lock  at  Dunkirk  Harbor  in  France  uses  a  complete  exchange 
system.  The  lock  dimensions  are  148.6  x  12.0  x  4.0  m  and  water  level  difference 
is  about  2.0  m.  It  is  navigable  for  push  tows  of  3,000-ton  barges  and  its  complete 
water  exchange  requires  12  minutes.  Since  the  lock’s  salt  prevention  measure  is 
a  combination  of  water  exchange  and  salt-water  trough,  fresh  water  losses  may 
be  as  high  as  55  percent.  The  lock  has  been  running  well  over  its  30  years  of 
operation. 


Salinity  Intrusion  Problems  in  Norway 

The  problems  in  Norway  are  related  to  strong  currents  and  ice  formation  at 
surface.  Use  of  air  bubbles  is  found  to  be  effective  in  breaking  down  density 
stratification. 

The  Norwegian  coast  is  basically  relatively  steep,  with  high  mountains 
surrounding  long  and  deep  fjords,  with  a  brackish  layer  on  the  top.  Therefore, 
the  saltwater  intrusion  problem  as  experienced  in  other  countries  like  the 
Netherlands,  or  England,  is  not  so  common  in  Norway.  In  the  southeastern  part 
of  the  country,  mainly  in  the  Glomma  estuary  area,  we  find  saltwater  intrusion 
problems  which  are  somewhat  similar  to  what  is  typically  found  abroad.  In  the 
rest  of  the  country,  saltwater  intrusion  problems  are  of  a  different  kind.  Two 
types  of  problems  have  been  experienced  in  Norway  which  can  be  related  to 
navigation  problems  as  well.  Both  problems  are  connected  to  the  presence  of 
saltwater  intrusion  in  estuary  or  fjord  areas. 

The  first  problem  is  related  to  the  outflow  of  fresh  or  brackish  water  through 
a  narrow  entrance.  The  presence  of  the  saltier  water  masses  below  creates  a 
strong  density  stratification  which,  in  turn,  creates  a  strong  concentration  of 
surface  currents  in  the  upper  water  masses.  Since  the  area  under  consideration  is 
used  as  a  fishing  harbour  as  well,  these  strong  surface  currents  have  created  some 
navigation  problems  for  the  fishing  vessels  moving  in  and  out  of  the  harbour 
area.  In  particular,  the  maneuvering  conditions  can  be  difficult  when  incoming 
waves  from  the  open  sea  are  entering  the  harbour  area  and  then  meet  the 
opposing  brackish  layer  flow,  causing  the  waves  to  steepen  considerably. 

The  other  case  is  of  a  somewhat  different  kind.  In  the  winter  time  and  under 
calm  weather  conditions,  the  fresh  water  from  the  river  mouth  will  spread  out  in 
the  surface  layers  in  the  fjord  area.  Due  to  the  stable  density  stratification  caused 
by  the  saline  water  below,  the  surface  layer  will  be  exposed  to  cooling  and 
subsequent  ice  formation.  The  saline  water  stratifications  will  in  this  case  act  as 
a  preventor  for  transport  of  excess  heat  from  below.  This  ice  might,  in  turn, 
cause  difficulties  for  navigation  to  and  from  harbours  in  the  area. 

In  the  latter  case,  the  problem  was  solved  by  means  of  an  air-bubble  plant 
installed  in  the  river  mouth  area.  In  order  to  break  down  the  stable  density 
stratification,  the  intrusive  salt  water  in  the  river  mouth  area  was  forced  by  the  air 
bubble  plant  to  mix  with  the  fresh  water. 
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Syncrolift  Lock 


The  Syncrolift  Lock  described  in  Part  2  as  being  considered  for  installation 
in  the  Panama  Canal  presents  a  somewhat  different  challenge  for  preventing  salt 
water  intrusion.  Since  the  lock  is  filled  to  the  same  level  at  all  times,  water 
within  it  is  not  diluted  by  the  addition  of  fresh  lift  water. 

During  discussions  of  the  Syncrolift  lock  numerical  experiments  presented  in 
Part  5  of  this  report,  the  PCC-COE  design  team  formulated  a  conceptual  design 
for  preventing  salinity  intrusion.  Depicted  in  Figure  4-4,  it  consists  of  a  small 
gated  antechamber  at  the  upper  lock  position  in  the  dam,  separated  from  the  lake 
by  a  movable  gate.  After  the  lock  is  lifted  into  position  at  the  top  and  seated 
against  the  dam,  the  lock  gates  are  opened  so  that  saline  water  in  the  lock  and 
fresh  water  in  the  antechamber  begin  to  mix.  Salty  water  is  then  withdrawn  from 
the  bottom  of  the  antechamber  as  freshwater  is  released  at  the  surface  at  the  same 
flow  rate.  These  flows  are  continued  until  the  lock  water  is  freshened  to  the 
required  degree,  and  then  the  antechamber  gates  are  opened  to  permit  the  ship  to 
move  into  the  lake.  This  concept  is  modeled  in  Part  5  of  this  report. 


Summary  and  Conclusions 

a.  All  over  the  world  salt  intrusion  through  navigation  locks  connected  to 
ocean  is  being  increasingly  regarded  as  a  serious  threat  to  the  quality  of 
water  in  the  canals  that  link  ocean  with  the  locks. 

b.  There  are  several  ways  of  preventing  salt-water  intrusion  in  the 
navigation  locks.  The  selection  of  method  is  based  on  factors  such  as 
site  conditions,  cost,  volume  and  frequency  of  ship  traffic  etc. 

c.  The  most  effective  and  at  the  same  time  most  expensive  method  is 
exchanging  fresh  water  for  salt  water  within  the  lock  while  gates  are 
closed;  however,  it  can  also  lead  to  adverse  motion  of  vessels  within  the 
lock. 

d.  The  option  of  equalizing  levels  by  pumping  instead  of  gravity,  which 
only  concerns  the  salt  water  required  for  equalization,  is  not  very 
effective. 

e.  The  selective  withdrawal  of  salt  water  usually  calls  for  substantial 
alterations  or  additions  to  locks  and  may  require  fairly  large  quantities  of 
fresh  water. 

/  The  water  barrier  method  has  not  been  tested  thoroughly  but  should 
perform  about  the  same  as  the  pneumatic  barrier. 

g.  Pneumatic  barriers  constitute  a  simple  method  of  countering  salt 
intrusion.  This  method  will  not  keep  all  the  salt  out  but  may  reduce 
intrusion  quite  considerably.  Pneumatic  barriers  have  been  successfully 
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installed  at  several  locks  and  can  be  combined  with  other  systems.  The 
pneumatic  barrier  method  of  preventing  salt  intrusion  is  found  to  be 
suitable  for  practically  all  the  locks  adjoining  the  sea  in  the  Netherlands. 
Pneumatic  barriers  are  already  operating  in  13  locks. 

h.  Establishing  a  canal  on  the  upstream  side  of  the  lock  and  flushing  the 
canal  with  fresh  water  is  less  effective  than  the  other  methods  and 
involve  wasting  large  quantities  of  fresh  water.  It  is  nevertheless  the 
oldest  way  of  getting  rid  of  intrusive  salt.  It  can  be  retained  as  a  useful 
supplement  to  one  or  more  of  the  other  methods. 

i.  For  the  proposed  new  locks,  the  interests  served  by  any  measures  taken 
to  combat  salt  intrusion  are  considered  for  each  case  individually. 
Methods  other  than  the  pneumatic  barrier  system  might  be  chosen  based 
on  certain  considerations  and  the  locks  will  be  designed  accordingly. 

NOTE:  See  Appendix  E  for  non-cited  references  used  in  Part  4. 
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Figure  4-1 .  Bubble  Screen 


Figure  4-4.  Schematical  synerolift  lock  with  gated  salt  water  extraction  chamber 
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Modeling  Density  Flows 


The  two-dimensional  width-averaged  numerical  model,  CE-QUAL-W2,  was 
used  to  investigate  density  flow  from  the  uppermost  Gatun  Lake  lock  into  a 
navigation  channel.  A  channel  sump  for  collecting  and  retaining  salt  water  from 
the  locks  and  a  salinity  barrier  with  bottom  extraction  were  simulated.  The 
results  of  these  simulations  are  summarized  here. 


Model  and  Assumptions 

CE-QUAL-W2  is  a  time- varying,  mechanistic,  two-dimensional  (vertical- 
longitudinal)  hydrodynamic  and  water  quality  reservoir  model.  The  model  is 
capable  of  simulating  an  extensive  list  of  chemical  and  biological  and 
hydrodynamic  processes  that  affect  water  quality  of  a  reservoir.  This  includes 
the  impacts  of  density  stratification  due  to  temperature  or  salinity'  The  model 
assumes  lateral  homogeneity  and  is  best  suited  for  water-bodies  exhibiting 
longitudinal  and  vertical  water  quality  gradients.  The  model  has  been  applied  to 
rivers,  lakes,  reservoirs,  and  estuaries  and  simulates  the  effects  on  reservoir 
water.  Details  of  the  model  are  provided  by  Environmental  and  Hydraulics 
Laboratory  (1986)  and  Cole  and  Buchak  (1995). 

Two  computational  grids  were  developed  for  modeling  Gatun  Lake  salinity 
intrusion.  The  segment  distribution  (longitudinal  cell  lengths)  used  for  both 
computational  grids  is  shown  in  Figures  5-1  and  5-2.  The  ship  lift  was  250  m 
long  by  28  m  wide  by  10  m  deep,  the  sump  was  500  m  long  by  100  m  wide  by 
10  m  deep.  The  upper  entrance  canal  was  750  m  long  by  100  m  wide  by  10  m 
deep.  The  numerical  model  does  not  allow  moving  boundaries,  thus  a  ship  was 
not  simulated  in  the  grid.  The  first  grid  simulated  a  Syncrolift  lock  chamber 
opening  into  a  100-m-wide  canal  with  a  10-m-deep  by  500-m-long  salinity 
collection  sump  (Figure  5-3).  During  lock  operation,  the  sump  was  to  collect  salt 
water  and  then  drain  the  salt  water  back  to  the  sea.  The  second  grid  simulated  a 
barrier  gate  (Figure  5-4)  to  retain  escaping  salt  water  with  simultaneous 
extraction  through  a  bottom  withdrawal.  The  barrier  gate  was  located  35  m 
upstream  of  the  entrance  to  the  ship  lift  with  a  crest  set  at  3  m  deep.  The 
extraction  sump  was  10-m  deep  between  the  barrier  gate  and  the  ship  lift.  For  the 
extraction  simula-tions,  a  large  reservoir  width  (1000  m  at  the  upstream  end)  was 
simulated  outside  of  the  immediate  vicinity  of  the  ship  lift  to  provide  a  large 
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reservoir  surface  area.  This  allowed  the  withdrawal  of  water  without 
significantly  lowering  the  pool. 


Experimental  Conditions 

Several  data  files  were  assembled  to  describe  the  initial  salinity  conditions 
across  the  computational  grid.  In  the  simulated  part  of  Gatun  Lake,  the  salinity 
was  set  at  0,  while  salinity  within  the  ship  lift  was  homogeneous  and  initially  set 
at  5  ppt  and^lOppt  for  different  simulations.  Water  temperature  was  set  at  a 
constant  30  C  for  the  entire  grid  and  did  not  vary  during  model  operation. 

T"he  .™tial  condltlons  were  defined  to  simulate  the  instantaneous  opening  of 

the  ship  lift  gates.  The  model  runs  simulated  up  to  2  hours  of  time. 


Plans  Tested 


Tested  plans  included  several  variations  of  a  sump  to  capture  salt  water 
carried  out  of  the  lock  by  density  flows,  following  the  design  of  the  Lake 
Washington  Ship  Canal  saltwater  sump.  Two  plans  are  presented  here  a  Ion- 
sump  in  a  lock  approach  channel  with  a  salt  water  drain  in  the  end  nearest  the 
lock,  as  depicted  m  Figure  4-2,  and  a  gated  chamber  (near-field  barrier)  built  into 
the  lock  structure  to  capture  salt  water  and  replace  it  with  fresh  water  before  the 

forfhip  movement-  The  tetter  design  concept  was  developed  by 
the  LUL-PCC  study  team  during  the  testing  procedure.  For  both  cases  the 
bottom  salt  water  extraction  rate  was  set  at  60  m3/s. 


lock. 


Either  of  these  plans  can  be  used  with  the  existing  locks  or  the  Syncrolift 


Results 


Salt  Water  Sump.  The  dense  water  within  the  ship  lift  slumped  out  toward 

the  upstream  channel,  cascading  down  into  the  salinity  sump  in  the  approach 
channel  OjdeoJ).  The  salinity  filled  the  sump  with  concentrations  greater  than 
U.D  ppt.  Even  though  the  numerical  model  was  more  diffusive  than  the  actual 
conditions  would  be,  these  results  showed  that  more  fresh  water  would  likely  be 
used  to  drain  the  salinity  sump  than  would  be  used  for  a  conventional  lock 

operation.  Thus,  this  alternative  was  modified  to  simulate  the  near-field  barrier 
with  bottom  extraction. 


,  Barrier.  The  results  of  this  simulation  are  presented  in  Video ; 

and  Video  s.  The  salt  water  slumps  into  the  trap  when  the  lock  gate  is  opened 
and  is  withdrawn  through  the  bottom.  Fresh  water  from  the  approach  channel 
s  lms  in  over  the  barrier  gate  to  replace  the  salt  water  that  has  been  withdrawn 
The  velocity  of  the  freshwater  over  the  barrier  gate  prevents  salt  water  from 
escaping  into  the  approach  channel.  Video  4  shows  the  particle  tracks  from  this 
simulation.  This  alternative  captures  the  salt  water  and  removes  it  before  it  can 
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escape  into  the  lake,  but  introduces  currents  over  the  barrier  that  may  cause  ship 
motion  within  the  lock. 


Recommendations 

This  numerical  model  offers  the  capability  to  quickly  simulate  and  assess 
stratified  flow  during  lock  operation.  Once  it  has  been  verified,  it  can  be  used  in  a 
design-level  effort. 

Of  the  two  alternatives  simulated  in  the  model,  the  near-field  salinity  barrier 
with  bottom  extraction  appears  most  feasible.  The  canal  sump  for  trapping  salt 
water  would  likely  result  in  excessive  use  of  freshwater  to  flush  the  large  volume 
of  saltwater/freshwater  mixture  from  the  sump.  The  near-field  barrier  with 
bottom  extraction  showed  better  promise  of  salinity  exclusion  provided  by  the 
velocity  pattern  over  the  barrier.  Further  research  should  be  conducted  into  the 
extraction  alternative  to  determine  the  maximum  time  allowable  for  the  saltwater 
exchange  and  the  degree  to  which  this  alternative  can  meet  this  requirement. 
Effects  on  vessels  within  the  lock  and  variations  in  barrier  depth  and  extraction 
rates  and,  in  concert  with  other  alternatives,  such  as  in-lock  flushing,  should  be 
included  in  more  detailed  studies  for  design  purposes. 
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flTable  5.1 

BLake  Gatun  Zero-Order  Salinity  Intrusion  Model  vl.1  Sep 

I  Panama  Canal 

|  Flows  and  time  interval  are  based  on  the  13870 

I  number  of  lockages 

i  I  Average  Madden  Lake  and  Tributary 

B  Inflow  for  time  interval 


Time 

Interval 


0 


1 


2 


3 


Volume  - 
LG,  cf 


1.82E+1 1 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+1 1 


1.82E+1 1 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


Spill  -  ML, 
cf 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


Qpower  - 
ML,  cf 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


Trib  -  LG, 
cf 


Salinity  - 
ML,  cf 


Exchange  during  lift  operations 
Pedro  Miguel  Lock 


Phi 

portion  of  Sal 
Qlock,  cf  lock  loc 


Exchange  during  lift 
operations  Gatun  Lock 


Phi 

portion  of 

Qlock,  cf  lock  Salin 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+1 0 


4.53E+10 


mil 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


IKIHglill 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


IBH3XE1I 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


Chapter  5  Modeling  Density  Flows 


ange  during  lift  operations 
Pedro  Miguel  Lock  _ 


Phi 

portion  of  Sal 
fc,  cf  lock  loc 


■I 

SS93 


Exchange  during  lift 
operations  Gatun  Lock 


Phi 

portion  of 

Qlock,  cf  lock  Salin 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


Number  of  Syncro  Lift 
Lockages  per  time  interval 


Mass  of  Exchange  during  lift  operations 
salt,  lbs  Miraflores  Syncrolift _ 


Phi 

portion  of  Sal 
Qlock,  cf  lock  loc 


3.30E+08  |  8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74  E+09 


8.74E+09 


8.74  E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8. 74  E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8. 74  E+09 


8.74E+09 


8.74  E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8. 74  E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74  E+09 


8. 74  E+09 


8.74E+09 


8. 74  E+09 


8.74  E+09 


Exchange  during  lift  operations 
Gatun  Syncrolift 


Phi 

portion  of 

Qlock,  cf  lock  Salini 


8.74  E+09 


8.74E+09 


8.74  E+09 


8.74E+09 


8.74  E+09 


irwzrereiMiKTiiiMi 


8.74E+09 


8.74E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74  E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8. 74  E+09 


8.74E+09 


8. 74  E+09 


8.74E+09 


8. 74  E+09 


8. 74  E+09 


1.64E+08 


m 


of  Syncro  Lift 
per  time  interval 

7300 

Computation  of  Trib  Inflow  to  | 
Gatun  1 

:ring  lift  operations 
res  Syncrolift 

Exchange 

Gi 

during  lift  operations 
itun  Syncrolift 

Mass  of  salt,  lbs 

Phi 

>rtion  of 
lock 

Salinity  - 
lock,  ppt 

Qlock,  cf 

Phi 

portion  of 
lock 

Salinity 

lock,  ppt 

Total 
Inflow,  cf 

Less  Madden 
Releases,  cf 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

1.64E+08 

1.18E+11 

4.53E+10 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

\ 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

).500 

0.300 

8.74E+09 

0.500 

0.300 

).500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

3.500 

0.300 

8.74E+09 

0.500 

0.300 

).500 

0.300 

8.74E+09 

0.500 

0.300 

).500 

0.300 

8.74E+09 

0.500 

0.300 

).500 

0.300 

8.74E+09 

0.500 

0.300 

).500 

0.300 

8.74E+09 

0.500 

0.300 

).500 

0.300 

8.74E+09 

0.500 

0.300 

).500 

0.300 

8.74E+09 

0.500 

0.300 

>.500 

0.300 

8.74E+09 

0.500 

0.300 

).500 

0.300 

8.74E+09 

0.500 

0.300 

).500 

0.300 

8.74E+09 

0.500 

0.300 

).500 

0.300  | 

8.74E+09 

0.500 

0.300 

\ 

).500 

0.300 

8.74E+09 

0.500 

0.300 

'.500 

0.300 

8.74E+09 

0.500 

0.300 

1.500 

0.300 

8.74E+09 

0.500 

0.300 

^.500 

0.300 

8.74E+09 

0.500 

0.300 

.500 

0.300 

8.74E+09 

0.500 

0.300 

.500 

0.300 

8.74E+09 

0.500 

0.300 

.500 

0.300 

8.74E+09 

0.500 

0.300 

.500 

0.300 

8.74E+09 

0.500 

0.300 

.500 

0.300 

8.74E+09 

0.500 

0.300 

.500 

0.300 

8.74E+09 

0.500 

0.300 

.500 

0.300 

8.74E+09 

0.500 

0.300 

(Continued)^ 
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Table  5.1  (Concluded)  


Panama  Canal 


Flows  and  time  interval  are  based  on  the 

number  of  lockages _ _ 


Average  Madden  Lake  and  Tributary 
«  Inflow  for  time  interval 


Time 

Interval 


52 


53 


54 


55 


56 


57 


58 


59 


60 


61 


62 


63 


64 


65 


66 


67 


68 


69 


70 


71 


72 


73 


74 


75 


76 


77 


78 


79 


80 


81 


82 


83 


84 


85 


86 


87 


88 


89 


90 


91 


92 


93 


94 


95 


96 


97 


98 


99 


100 


Volume  - 
LG,  cf 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


lEEsaai 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


Spill  -  ML, 
cf 


5.73E+09 


5.73E+09 


5.73E+09 


Qpower  - 
ML,  cf 


Trib  -  LG, 
cf 


■aus&JI 


IKK^SSa 


1.82E+11 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


imssmi 


5.73E+09 


5.73E+09 


5.73E+09 


EE3EI 


5.73E+09 


5.73E+09 


5.73E+09 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


IKEESm 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


1.82E+11 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


5.73E+09 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


6.66E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


EMSEII 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


4.53E+10 


Salinity  - 
ML,  cf 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


Exchange  during  lift  operations 
Pedro  Miguel  Lock 


Phi 

portion  of  Sal 
Qlock,  cf  lock  loc 


0.500  0.100  5.45E+10 


5.12E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


0.500  0.100  5.45E+10 


5.12E+10  I  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10  I 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  0.100  5.45E+10 


5.12E+10  0.500  I  0.100  I  5.45E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


5.12E+10 


Exchange  during  lift 
operations  Gatun  Lock 


Phi 

portion  of 

lock  Salin 


0.500 


0.500 


0.500 


0.500 


0.500 


0.500 


0.500 


0.500 


0.500 


5.12E+10 

0.500 

5.12E+10 

0.500 

5.12E+10 


5.12E+10 


5.45E+10 


5.45  E+ 10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


BE3E1I 


5.45E+10 


5.45E+10 
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unng  lift  operations 
>  Miguel  Lock 


Phi 

>ortion  of  Sal 
lock  loc 


m 

HRS! 


Exchange  during  lift 
operations  Gatun  Lock 


Phi 

portion  of 

Qlock,  cf  lock  Salin 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


5.45E+10 


Number  of  Syncro  Lift 
Lockages  per  time  interval 


Mass  of  Exchange  during  lift  operations 
salt,  lbs  Miraflores  Syncrolift  _ 


Phi 

portion  of 
Qlock,  cf  lock 


Computal 


r.TH  [«■  •]  • 


Exchange  during  lift  operations 
_ Gatun  Syncrolift 


Phi 

portion  of 

Qlock,  cf  lock  Salinity  -  lock, 


Mass  of  salt,  lbs 


HR 


Total 
Inflow,  cf 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74  E+09 


8.74E+09 


msm 


8. 74  E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8. 74  E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8. 74  E+09 


8.74E+09 


8.74E+09 


8.74  E+09 


8. 74  E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8. 74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8. 74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8. 74  E+09 


8.74  E+09 


8.74E+09 


8. 74  E+09 


8.74E+09 


8. 74  E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8.74E+09 


8.74  E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8. 74  E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8. 74  E+09 


8.74E+09 


8.74E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74  E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


8.74E+09 


m 


jmber  of  Syncro  Lift 
cages  per  time  interval 

7300 

Compi 

itation  of  Trib  Inflow  to 
Gatun 

jge  during  lift  operations 
iraflores  Syncrolift 

Exchange 

Gi 

during  lift  operations 
atun  Syncrolift 

Mass  of  salt,  lbs 

cf 

Phi 

portion  of 
lock 

Salinity  - 
lock,  ppt 

Qlock,  cf 

Phi 

portion  of 
lock 

Salinity 

lock,  ppt 

Total 
Inflow,  cf 

Less  Madden 
Releases,  cf 

59 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

39 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

99 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

)9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

)9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

)9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

)9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

)9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

>9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

)9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

>9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

>9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

>9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

r 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

\ 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

i  8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

9 

0.500  ! 

0.300 

8.74E+09 

0.500 

0.300 

) 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

7 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

7 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

) 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

) 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

T 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

T 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

? 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

i 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 

0.500 

0.300 

8.74E+09 

0.500 

0.300 
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Table  5.2 

Zero-Order  Miraflores  Lake  Salinity 


Modelvl.0  Aug  1999 


Flows  are  based  on  the 
number  of  lockages 


Inflow  from  Pedro 
Miguel  Lock 
Operation 


Exchange  during  lift  Operations 


Qlock,  cf 


9.55E+07 


9.55E+07 


9.55E+07 


9.55E+07 


9.55E+07 


9.55E+07 


9.55E+07 


9.55E+07 


9.55E+07 


9.55E+07 


9.55E+07 


9.55E+07 


9.55E+0 


8.23E+08  9.55E+07 


8.23E+08  9.55E+07 


8.23E+08  9.55E+07 


8.23E+08  9.55E+07 


8.23E+08  9.55E+07 


8.23E+08  9.55E+07 


8.23E+08  |  9.55E+07 


8.23E+08  9.55E+07 


8.23E+08  9.55E+07 


8.23E+08  9.55E+07 


8.23E+08  9.55E+07 


8.23E+08  9.55E+07 


Salinity, 

ppt 


Qlock,  cf 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1 .40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


1.40E+08 


Phi, 
portion 
of  lock 


Salinity  - 
lock,  ppt 


Exchange  during  lift 
operations 

Phi, 

portion  Salinity  -  Mass  of 
Qlock,  cf  of  lock  lock,  ppt  salt,  lbs 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


1.49E+08  0.500 


6.98E+06 


25 

0.49 

8.23E+08 

9.55E+07 

0.000 

1.40E+08 

0.500 

0.10 

1.49E+08 

0.500 

1.500 

26 

0.49 

8.23E+08 

9.55E+07 

0.000 

1.40E+08 

0.500 

0.10 

1.49E+08 

0.500 

1.500 

27 

0.49 

8.23E+08 

9.55E+07 

0.000 

1.40E+08 

0.500 

0.10 

1.49E+08 

0.500 

1.500 

28 

0.49 

8.23E+08 

9.55E+07 

0.000 

1.40E+08 

0.500 

0.10 

1.49E+08 

0.500 

1.500 

29 

0.49 

8.23E+08 

9.55E+07 

0.000 

1.40E+08 

0.500 

0.10 

1.49E+08 

0.500 

1.500 

30 

0.49 

8.23E+08 

9.55E+07 

0.000 

1.40E+08 

0.500 

0.10 

1.49E+08 

0.500 

1.500 

31 

0.49 

8.23E+08 

9.55E+07 

0.000 

1.40E+08 

0.500 

0.10 

1.49E+08 

0.500 

1.500 

32 

0.49 

8.23E+08 

9.55E+07 

0.000 

1.40E+08 

0.500 

0.10 

1.49E+08 

0.500 

1.500 
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Table  5.2  (Continued) 


Flows  are  based  on  the 
number  of  lockages 


Inflow  from  Pedro 
Miguel  Lock 
Operation 


Volume  Salnity 

cf  Qlock,  cf  ppt 


Salinity  - 
ML,  ppt 


Exchange 


Qlock,  cf 


during  lift 


Phi, 
portion 
of  lock 


operations 


Salinity  - 
lock,  ppt 


8.23E+08  |  9.55E+07  |  0.000  |  1.40E+08 


8.23E+08  |  9.55E+07  |  0.000  |  1.40E+08 


8.23E+08  |  9.55E+07  |  0.000  |  1.40E+08 


8.23E+08  I  9.55E+07  I  0.000  I  1.40E+08 


8.23E+08  |  9.55E+07  |  0.000  |  1.40E+08 


8.23E+08  |  9.55E+07  I  0.000  I  1.40E+08 


8.23E+08  |  9.55E+07  |  0.000  |  1.40E+08 


8.23E+08  |  9.55E+07  |  0.000  |  1.40E+08 


8.23E+08  I  9.55E+07  I  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  |  9.55E+07  j  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  I  9.55E+07  I  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  |  9.55E+07  |  0.000 


8.23E+08  I  9.55E+07  I  0.000 


8.23E+08  |  9.55E+0 


8.23E+08  |  9.55E+0 


8.23E+08  |  9.55E+0 


8.23E+08  |  9.55E+0 


8.23E+08  |  9.55E+0 


8.23E+08  I  9.55E+0 


Exchange  during  lift 
operations 
“PhiT 

portion  Salinity  -  Mass  of 
Qlock,  fc  of  lock  lock,  ppt  salt,  lbs 


1.49E+08  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1 .49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  I  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  |  0.500 


1.49E+08  I  0.500 
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Figure  5-1.  Details:  Segment  distribution  for  Lake  Gatun  computational  grids 
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Figure  5-4.  Longitudinal  Profile  of  width  used  in  salinity  trap  and  extraction  calculations 


6  Modeling  Salt  Movement 
Across  Gatun  Lake 


Objective 

The  objective  of  this  portion  of  the  work  is  to  examine  how  salinity  moves 
into  Gatun  Lake  from  the  locks  if  the  salt  water  is  not  captured  immediately  in 
front  of  the  lock.  The  results  will  help  Canal  managers  compare  capacity 
expansion  and  fresh  water  conservation  plans  and  to  assess  the  need  for  remedial 
measures  to  maintain  the  desired  salinity  level  across  the  lake. 


Models  and  Assumptions 

The  hydrodynamics  of  Gatun  Lake  were  modeled  in  two  dimensions  by 
using  the  TABS-MD  numerical  modeling  system  of  the  U.S.  Army  Corps  of 
Engineers.  TABS-MD  consists  of  a  number  of  computer-based  numerical  models 
for  open-channel  flow  and  transport.  For  the  effort  described  here,  the  depth- 
averaged  two-dimensional  hydrodynamic  model  RMA-2  and  companion 
transport  model  RMA-4  were  used.  RMA-2  provides  a  finite  element  solution  of 
the  Reynolds  form  of  the  two-dimensional,  constant  density  Navier-Stokes 
equations  for  turbulent  open  channel  flows.  RMA-4  solves  the  two-dimensional 
advection-diffusion  equation  for  multiple  dissolved  constituents.  (See 
Appendix  D  for  discussion  of  the  governing  equations  and  assumptions.) 

Salinity  flow  out  of  the  locks  will  have  a  significant  three-dimensional 
density-driven  structure  that  will  be  fully  captured  only  by  a  three-dimensional 
model  that  has  been  verified  to  field  observations.  However,  the  objectives  of 
this  effort  can  be  met  by  application  of  two-dimensional  models  as  described  in 
this  section  and  the  preceding  section  of  this  report.  The  two-dimensional  width- 
averaged  numerical  model  of  Part  5  accurately  depicts  the  density-driven  vertical 
flow  structure  near  the  lock,  and  the  depth-averaged  results  presented  in  this 
section  will  be  more  accurate  farther  from  the  locks,  where  diffusive  transport 
dominates,  and  where  the  waterway  is  much  wider.  The  models  used  here  have 
not  been  verified  to  field  observations,  but  the  calculations  are  physics-based  and 
have  been  applied  to  many  sites  by  the  Corps  of  Engineers  and  others  (Heath 
et  al.,  1999)  so  the  results  are  considered  qualitatively  reliable  and  quantitatively 
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useful  when  one  set  of  model  far-field  results  is  compared  with  another,  e.g., 
existing  conditions  versus  the  Syncrolift  plan  and  wet  versus  dry  season.. 

For  the  purposes  of  this  study,  a  mesh  representing  Gatun  Lake  and  the 
Panama  Canal  with  its  Atlantic  and  Pacific  entrances  was  developed  (See 
Figure  6-1).  This  mesh,  constructed  using  maps  provided  by  the  Panama  Canal 
Commission,  was  employed  by  RMA-2  and  RMA-4  to  solve  the  governing 
equations  across  the  area  of  interest.  Elevations  of  the  lake  bottom  are  shown  by 
shading  in  Figure  6-1. 


Experimental  Conditions 

The  hydrodynamic  calculations  assumed  that  38  lockage  operations  occur  per 
day  at  both  Atlantic  and  Pacific  entrances;  thus,  boundary  rating  curves  were 
applied  at  both  entrances  to  provide  outflows  equal  to  the  lift  volume  required  to 
execute  that  many  lockages.  These  curves  remove  water  from  the  lake  by  taking 
into  account  the  volume  of  the  lock  and  the  water  surface  elevation  in  the  vicinity 
of  the  lock  at  the  time.  Outflows  over  the  spillways  were  specified  as  zero  during 
the  dry  season  and  were  calculated  from  water  surface  elevations  during  the  wet 
season.  Fresh  water  inflows  to  the  system  were  applied  at  four  locations: 

Madden  Lake,  Gatun  River  at  Ciento,  Ciri  Grande  River  at  Los  Canones,  and 
Trinidad  River  at  El  Chorro.  Evaporation  was  considered,  as  well  as  possible 
outflows  from  the  lake  for  municipal  uses  and  power  generation.  Model  runs 
were  performed  for  a  two  month  period,  which  was  sufficient  to  reach  a  quasi¬ 
steady-state  condition. 

The  existing  conditions  were  evaluated  for  two  seasons,  wet  and  dry.  During 
dry  season,  which  lasts  approximately  5  to  6  months  out  of  the  year,  the  inflows  * 
at  the  locations  mentioned  above  are  significantly  lower  than  those  present 
during  wet  season.  As  a  result,  the  Gatun  Dam  spillway  does  not  operate  during 
this  time.  On  the  other  hand,  the  wet  season  inflows  are  up  to  4  times  larger  than 
the  diy  season,  thereby  increasing  the  water  surface  elevation  of  the  lake  and 
opening  the  spillway  for  operation.  For  the  modeling  process,  Gatun  Lake  was 
kept  at  a  water  surface  elevation  of  26  m  (85  feet)  during  dry  season  and  at  27  m 
(87.5  feet)  for  wet  season. 

Based  on  field  observations  and  the  calculations  shown  in  Part  3  and 
Appendix  A,  a  salinity  of  0.083  ppt  was  applied  as  a  boundary  condition  at  the 
Gatun  Dam  locks.  Likewise,  a  salinity  value  of  0.13  was  applied  as  the  boundary 
condition  at  the  Pedro  Miguel  Locks.  For  the  Syncrolift  lock  plan,  a  unit 
concentration  was  specified  at  the  lock  location.  Since  the  depth-integrated 
model  does  not  calculate  water  density  and  density  currents,  this  unit 
concentration  can  be  approximately  any  specified  salinity  in  the  lock  by  simply 
multiplying  the  model  results  by  the  specified  lock  concentration.  The  existing 
lock  results  may  be  similarly  extended  to  other  in-lock  concentrations  by  a  ratio 
of  the  specified  concentration  to  the  values  listed  above.  Extending  these  results 
to  other  lock  salinities  and  comparing  the  syncrolift  tracer  results  to  existing  lock 
results  is  only  approximate,  since  diffusion  rate  in  the  model  is  proportional  to 
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the  concentration  gradient.  Thus,  more  rapid  initial  diffusion  will  be  computed 
for  the  1 .0  ppt  lock  condition  than  for  the  0.08  ppt  lock  condition. 


Plans  Tested 

Two  lock  configurations  -  the  existing  locks  and  an  additional  lock 
(Syncrolift)  adjacent  to  each  set  of  locks  -  and  two  hydrologic  conditions  -  wet 
and  dry  seasons  -  were  examined.  In  summary,  four  model  runs  were  made: 

a.  Existing  conditions  during  dry  season 

b.  Existing  conditions  during  wet  season. 

c.  Syncrolift  lock  (only)  at  both  entrances  during  dry  season. 

d.  Syncrolift  lock  (only)  at  both  entrances  during  wet  season. 


Results 

Figure  6-2  shows  a  snap-shot  comparison  of  wet  and  dry  seasons  for  existing 
conditions  in  the  vicinity  of  the  Pedro  Miguel  Locks  after  two  months.  The  dry 
season  permits  salinity  to  intrude  much  further  than  the  wet  season,  and  salinity 
contours  are  seen  to  intrude  nearly  to  Madden  Lake.  However,  during  the  wet 
season  the  salt  is  swept  back  to  the  locks  and  out  via  the  spillway. 

Figure  6-3  shows  a  comparison  between  wet  and  dry  seasons  at  the  Atlantic 
side  of  Gatun  Lake.  As  seen  in  the  figure,  salinity  concentrations  remain  higher 
on  this  side  more  readily  during  wet  season,  since  the  salt  is  not  confined  in  a 
narrow  channel  as  on  the  Pacific  side.  Nevertheless,  salinity  values  do  not 
exceed  0.027  ppt  inside  the  lake  even  in  the  worst  case,  which  corresponds  to  the 
lack  of  observed  salinity  buildup  in  Gatun  Lake.  The  Gatun  spillway  is  seen  to 
play  a  key  role  during  wet  season  in  removing  higher  salinity  water  from  the 
lake,  for  the  low  concentration  plume  is  drawn  to  the  spillway. 

Similar  scenarios  play  out  when  adding  the  Syncrolift  locks  to  the  system. 
Figure  6-4  portrays  how  tracer  in  the  Syncrolift  locks  diffuses  from  the  Pedro 
Miguel  Locks  for  wet  and  dry  seasons.  The  contours  on  Figures  6-4  and  6-5 
show  the  fraction  of  the  initial  lock  concentration  present  in  different  areas  of  the 
mesh  after  the  model  has  run  for  two  months. 

Figure  6-5  shows  a  closer  look  at  the  Gatun  Locks  area,  where  we  again  see 
that  the  Gatun  spillway  is  important  in  alleviating  salinity  concentrations  during 
the  wet  season.  We  should  note,  though,  that  the  increased  tracer  concentrations 
observed  in  the  vicinity  of  the  Gatun  Locks  during  dry  season  are  mostly  a 
product  of  tracer  diffusing  up  through  the  lake  from  the  Pedro  Miguel  Locks. 
Film  loops  accompanying  this  report  show  this  situation  clearly. 
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Recommendations  for  Design  Phase 

This  modeling  approach  depicts  what  happens  to  salt  that  is  not  captured 
close  to  the  locks.  If  the  final  design  for  expanding  canal  capacity  ensures  that 
no  more  salt  water  will  be  released  from  the  locks  than  under  existing  conditions, 
Gatun  Lake  need  not  be  modeled  during  the  design  phase.  However,  this 
modeling  approach,  either  in  two-dimensions  or  three-dimensions,  can  be  used  to 
predict  salinities  in  the  waterways  at  the  sea  ends  of  both  locks  to  improve  the 
boundary  conditions  assumed  in  Part  3  of  this  report. 

If  Gatun  Lake  transport  is  modeled  in  the  design  phase,  the  three- 
dimensional  version  of  these  models  (TABS-10)  should  be  used  to  ensure  that 
densimetric  transport  is  reproduced.  The  model  should  be  verified  to  field 
measurements  of  current  speed  and  direction  and  salinity,  with  low-threshold 
instruments  employed  to  successfully  measure  the  small  flow  speeds  and 
salinities  of  the  lake. 
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Figure  6-2.  Salinity  Profiles  at  Pedro  Miguel  Locks  and  Surrounding  Areas  for  dry  and 
wet  seasons 
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Figure  6-4.  Concentration  patterns  for  the  Syncrolift  lock  plan  in  dry  and  wet  seasons 
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Figure  6-5.  Concentration  patterns  at  Gatun  Locks  for  dry  and  wet  seasons  with  both 
existing  locks  and  Syncrolift  in  operation 
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A  review  of  salinity  control  measures  at  locks  worldwide  shows  a  variety  of 
methods  being  employed,  including  salt  water  barriers,  sumps  to  collect  and 
drain  salty  water,  and  mixing  devices  to  limit  intrusion.  A  detailed  evaluation  of 
the  Hiram  Chittenden  Locks  in  Seattle,  Washington,  USA,  shows  that  a  salt 
water  sump  with  a  hinged  barrier  is  an  effective  approach  to  salt  water  intrusion 
control. 

Salt  water  is  slightly  more  dense  than  freshwater,  causing  it  to  flow  under 
fresh  water  in  a  density  current  unless  mixed  by  strong  currents,  ship  propwash, 
or  artificial  mixing  devices  such  as  bubble  screens.  Density  currents  may  lead  to 
greater  saltwater  intrusion  than  would  occur  in  well-mixed  waters,  but  offers  an 
opportunity  to  capture  some  of  the  salt  water  before  it  intrudes  into  the 
freshwater  zone  upstream  of  the  locks. 

A  set  of  tools  has  been  developed  that  supports  evaluation  of  salinity 
intrusion  and  freshwater  consumption  through  the  Panama  Canal  locks  for  Gatun 
Lake.  They  include: 

a.  A  simple  mixing  analysis  coded  in  spreadsheet  form  that  provides  zero¬ 

dimensional  models  for  the  salinity  of  individual  locks  and  for  Miraflores 
Lake  and  Gatun  Lake  plus  a  freshwater  consumption  model  for  Gatun 
Lake. 

b.  A  depth-integrated  two-dimensional  numerical  model  for  dispersive 

transport  of  salt  water  into  Gatun  Lake. 

c.  A  width-integrated  two-dimensional  numerical  model  for  densimetric 
advection  and  dispersion  of  salty  water  from  the  locks  into  a  channel. 

The  numerical  models  employ  the  graphical  user  environment  of  the  U.  S. 
Army  Corps  of  Engineers  Surface  Water  Management  System  (SMS)  to  manage 
the  models  and  display  their  results.  The  zero-dimensional  models  have  been 
adjusted  to  yield  results  roughly  comparable  to  limited  field  observations  of 
salinity  in  the  Panama  Canal.  The  two-dimensional  models  are  physics-based 
and  have  been  successfully  employed  in  many  applications;  however,  they  have 
not  been  verified  to  field  observations  in  Gatun  Lake. 

These  tools  have  been  applied  in  a  feasibility-level  evaluation  of  the 
following  alternatives  for  expansion  of  Panama  Canal  capacity. 
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a.  Use  of  holding  ponds  to  conserve  fresh  water  by  recirculating  lift  water. 

b.  Installation  of  a  Syncrolift  lock  chamber  adjacent  to  the  Miraflores  and 
Gatun  Locks. 

c.  A  sump  in  the  channel  upstream  of  the  locks  to  capture  more  dense  salty 
water  and  drain  it  before  it  reaches  the  lake. 

d.  A  gated  chamber  upstream  of  the  locks  to  hold  salty  water  while  it  is 
drained  to  the  sea  and  replaced  with  freshwater. 

Analyses  using  the  tools  listed  above  show  that: 

a.  Freshwater  consumption  can  be  reduced  by  up  to  about  60%  by 
recirculating  lift  water  through  holding  ponds,  albeit  at  the  cost  of 
increasing  salinity  in  the  ocean  approach  channels  and  ultimately  in  the 
uppermost  locks.  If  this  option  is  employed,  the  potential  for  increased 
salinity  intrusion  and  necessary  remedial  measures  should  be  addressed. 

b.  A  sump  in  the  channel  upstream  of  the  locks  can  capture  salt  water 
through  selective  withdrawal,  and  is  most  effective  for  higher  salinity 
(and  thus  more  dense)  water  in  the  locks;  however,  higher  salinity  in  the 
locks  (now  about  0.1  ppt  in  the  uppermost  Miraflores  and  Gatun  Locks) 
risks  greater  salinity  contamination  of  the  lake  if  any  lock  water  escapes 
the  sump.  A  sump  requires  some  sort  of  confined  channel  upstream  of 
the  lock  to  train  the  densimetric  flow  out  of  the  lock. 

c.  A  salinity  of  greater  than  0.3  ppt  in  the  uppermost  locks,  including  a 
Syncrolift  lock,  risks  unsatisfactory  buildup  of  salt  water  in  Gatun  Lake 
unless  some  means  of  capturing  most  of  the  salt  water  is  incorporated  in 
the  facility. 

d.  A  gated  chamber  on  the  lake  side  of  the  uppermost  lock  offers  the  option 
of  draining  salty  water  from  the  lock  and  replacing  it  with  fresh  water 
before  opening  for  ship  passage.  This  technique  provides  positive 
salinity  control  at  the  locks,  either  standard  or  Syncrolift,  at  the  cost  of  a 
delay  in  the  vessels’  exit  of  the  lock. 


Recommendations 

We  recommend: 

a.  Further  consideration  of  recirculating  lift  water  through  holding  ponds  as 
a  means  of  freshwater  conservation. 

b.  Further  consideration  of  a  gated  chamber  at  the  uppermost  lock  entrance 
on  each  side  to  reduce  salinity  intrusion  into  Gatun  Lake  if  either  the 
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Syncrolift  lock  or  the  holding  pond/recirculation  alternatives  are 
considered  for  design. 

c.  Use  of  the  tools  developed  here  in  the  design  phase  of  any  canal 

improvement  plans.  They  should  be  verified  to  a  comprehensive  field 
data  set  of  water  levels,  velocities,  salinities,  and  temperatures  before  use 
in  a  design  effort. 

d.  Addition  of  a  model  of  the  sea  approaches  to  calculate  the  salinity  there 

under  any  alternative  that  might  alter  the  salinity,  such  as  the 
recirculation/holding  pond  option. 

e.  Field  data  collection  in  support  of  further  modeling  should  be  synoptic  to 

the  maximum  degree  possible,  should  include  acoustic  Doppler  current 
profiling,  and  should  include  salinity  meters  capable  of  measurements  to 
hundredths  of  a  ppt.  If  field  instruments  are  not  available  to  meet  this 
requirement,  water  samples  may  be  collected  and  analyzed  in  a 
laboratory  to  the  required  accuracy. 
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Appendix  A 

Salt  Loading  Example 

Calculations 


Introduction 

The  objective  of  this  portion  of  the  work  is  to  compute  the  amount  of  salt 
transferred  from  the  ocean  to  the  two  fresh  water  lakes.  The  computations  are 
based  on  the  simple  mixing  model  of  Chapter  3  of  this  report. 

The  local  datum  used  at  Panama  Canal  is  the  Precise  Level  Datum  (PLD). 
Mean  sea  level  (msl)  is  1 .0  foot  (0.30  m)  above  the  PLD  on  the  Pacific  side  and 
0.204  feet  (0.06  m)  on  the  Atlantic  side.  For  purposes  of  this  study,  the  differ¬ 
ence  between  PLD  and  msl  is  ignored  and  all  levels  are  referred  with  respect  to 
the  msl.  This  assumption  does  not  make  any  significant  change  in  the  results. 

The  two  locks  on  the  Pacific  Ocean  side  are  the  Miraflores  Locks.  The  lock 
between  the  Miraflores  Lake  and  the  Gatun  Lake  is  called  Pedro  Miguel  lock. 
The  three  locks  on  the  Atlantic  Ocen  side  are  called  the  Gatun  Locks.  For 
convenience  of  reference,  the  six  locks  in  the  system  are  given  the  following 
notation. 

LI :  Lowest  Miraflores  lock  on  the  Pacific  side 

L2:  Lock  between  Miraflores  Lake  and  LI. 

L3:  Lock  between  Miraflores  Lake  and  Gatun  Lake.  (Pedro  Miguel) 

L4:  Highest  Gatun  lock  directly  connected  to  Gatun  Lake. 

L5:  Middle  Gatun  lock  on  Atlantic  Side 

L6:  Lowest  Gatun  lock  connected  to  Atlantic  Ocean. 

A  schematic  longitudinal  section  of  all  the  locks  is  given  in  Figure  Al.  This 
figure  provides  dimensions,  lifts  and  relevant  levels  of  all  the  locks. 

Physical  characteristics  of  the  six  locks  at  the  Panama  Canal  are  given  in 
Table  Al.  Salinity  values  obtained  from  limited  field  data  are  given  in  Table  A2. 
Assumptions  made  for  computations  are  listed  in  Table  A3.  Notation  used  in 
salinity  computations  is  given  in  Table  A4. 

Calculations  presented  in  Appendix  A  were  performed  in  non-SI  units.  To  convert  values  to 
SI  units,  use  the  conversion  factors  given  on  page  ix  of  this  report. 
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Table  A1 

Physical  Character^ 

>tics  of 

Panama  Canal  Locks 

Lock  # 

Lock  Width  & 
Length  (ft) 

Lift 

Jft} _ 

Lift  Volume 
(eft) 

Bed  Level 
(ft) 

Sill  Level 
(ft) 

Height  of  Sill 
above  Bed 
(ft) 

Volume  of  water 
below  Sill  (eft) 

L  1 

110  x  1051 

26 

3.005  x  106 

-52.0 

-18.3 

33.7 

3.896  x  106 

L  2 

110  x  1071 

28 

IHB1 

-  20.33 

+  11.3 

31.63 

3.487  x  106 

L  3 

110  x  1071 

31 

B99KE91 

+  11.0 

+  37.3 

26.3 

3.334  x  106 

L  4 

110  x  1071 

28 

HBI 

+  13.0 

+  37.3 

24.3 

1.566  x  106 

L  5 

110  x  1051 

29 

- 15.33 

+  15.0 

30.0 

L  6 

110  x  1051 

28 

3.237  x  106 

-  44.3 

-13.6 

30.7 

Note:  All  levels  are  given  with  reference  to  the  mean  sea  level(msl). 

The  difference  between  the  mean  sea  level  and  the  Panama  Canal  Precise  Datum 
Level  is  ignored. 


Table  A2 

Salinity  values  derived  from  field  observations 

Location 

Salinity  (ppt) 

On  the  sea  side  of  lowest  lock  on  Pacific  side  (lock  LI) 

10.0 

Inside  lock  LI 

8.0 

Inside  lock  L2 

1.  5 

Miraflores  Lake 

1.0 

Gatun  Lake 

0.0 

Limitations  of  field  data  on  salinity: 

a.  The  time  of  measurement  is  not  given.  Hence  range  and  stage  of  tide  are 
not  known. 

b.  The  status  of  water  level  in  the  lock  is  not  given;  whether  the  level  was 
after  full  lift  or  after  the  ship  entered  at  a  lower  water  elevation. 

c.  The  size  of  ship  is  not  reported. 

d.  Salinity  readings  are  given  to  the  first  decimal  value. 


A2 


Appendix  A  Salt  Loading  Example  Calculations 


Table  A3 

Assumptions  made  in  salt  loading  computations 


Assumption  _ 


Tidal  variation  in  sea  is  not  accounted  for.  All  computations  are  made  for 
mean  sea  level  (msl).  _ 


The  mean  sea  level  is  assumed  to  be  the  same  as  Precise  Level  Datum 

(PLD). _ _ _ _ _ 


Number  of  lockages  is  assumed  to  be  the  same  as  the  number  of  ship 
transits.  _ 


Blockage  ration  of  0.8  is  assumed  for  calculating  submerged  volume  of  ship 
from  the  length,  beam  and  draft  dimensions.  The  average  submerged  volume 
of  a  ship  is  assumed  to  be  1,140,000  eft. _ 


Fifty  percent  mixing  is  assumed  when  lock  gates  are  opened. 


Ocean  salinity  is  assumed  to  be  35  parts  per  thousand.  The  local  salinity  near 
the  lowest  gate  is  expected  to  be  much  less  because  fresh  water  from  the 
lake  is  added  to  the  sea  for  each  ship  transit  and  a  transition  zone  exists 
between  the  lowest  locks  and  the  open  sea.  Hence  it  is  assumed  that  salinity 
just  outside  the  lowest  lock  is  10  ppt  and  salinity  inside  the  lowest  lock  is  8 
PPt. _ - _ 


Volume  of  water  in  Miraflores  Lake  is  670  million  eft. 


Net  transfer  of  fresh  water  from  the  lakes  to  the  sea  is  equal  to  52  million 
gallons  per  transit.  _ 


Table  A4 
Notations  used 


in  salinity  computations 


Notation 


Volumes 


Parameter 


Magnitude 


Volume  of  LI  below  msl  _ 


Volume  of  LI  above  sill  after  full  lift  _ 


Average  volume  of  submerged  ship  _ 


Volume  of  water  in  LI  below  msl  with  ship  inside 


Lift  Volume  of  LI  _ 


The  volume  of  water  from  LI  that  will  mix  with  water  in  L2.  (i.e.  volume  of  water  above 
sill  on  the  L2  side  minus  the  submerged  ship  volume  with  water  in  L2  (i.  e.  the  volume 
above  sill  on  the  L2  side  minus  the  submerged  ship  volume 


Volume  of  water  in  L2  above  bed  at  +26  ft  _ 


Volume  of  water  in  L2  at  +26  feet  after  ship  enters  L2  _ 


Lift  volume  in  L2  _ 


Volume  of  water  in  L2  at  +54  ft  without  ship  inside 


Volume  of  water  in  L2  at  +  54  ft  with  ship  inside  _ 


Volume  of  water  in  L3  at  +54  ft  without  ship  inside 


Volume  of  water  in  L3  at  +  54  ft  with  ship  inside  _ 


Assumed  salinity  outside  LI  _ 


Measured  salinity  inside  LI  _ _ 


Measured  salinity  inside  L2  _ _ 


Computed  salinity  in  LI  after  operation  LI /I 


Measured  salinity  in  Miraflores  Lake 


Computed  salinity  in  L2  after  operation  L2/1  _ 


Computed  salinity  in  L2  after  operation  L2/2  _ 


Salinity  in  L2  after  ship  leaves  L2  and  enters  Miraflores  Lake _ _ 


Salinity  in  L3  after  ship  enters  L3 


Salinity  in  L3  after  raising  ship  to  +85  ft  _ 


6,011,720  eft 


5,121,523  eft 


1,140,000  eft 


4,871,720  eft. 


3,005,860  eft 


3,981,523  eft 


5,458,137  eft 


4,318,137  eft 


3,298,680  eft 


8,756,817  eft 


7,616,817  eft. 


5,056,830  eft 


3,916,830  eft 


10.0  ppt. 


18.0  ppt _ 


[1-5  ppt 


[5.52  ppt 


1  PPt 


3.2  ppt 


1.28  ppt 


1.003  ppt 


0.25 
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Computations  for  ship  transiting  from  Pacific 
Ocean 

Computations  for  Lock  -  LI 

Data 

Size  of  LI 
Lift 

Bed  Elevation  below  msl 
Volume  of  LI  below  msl  (VI) 

Lock-2  side  sill  level  of  LI 
Water  depth  above  sill  at  full  lift 
Volume  of  LI  above  sill  (V2) 

Average  volume  of 
submerged  ship  (V3) 

Assumed  salinity  outside  LI  (SI) 

Assumed  salinity  inside  LI  (S2) 

Assumed  salinity  of  water  in  L2  (S3) 

Initial  Status: 

Lock  gates  to  the  sea  are  open. 

Ship  is  waiting  out  in  the  sea  to  get  in. 

Sea  level  =  msl. 

Water  Elevation  in  LI  =  msl 

Operation  Ll/1 

Bring  the  ship  in  LI  from  sea  and  close  gates. 


1051  feet  x  110  feet 
26  feet 
-  52  feet 

=  llOx  1051  x52 
=  6,01 1,720  cubic  feet. 

- 18.3  feet 

=  (26  +  18.3)  =  44.3  feet 
=  110x1051  x  44.3 
=  5,121,523  cubic  feet 

=  1,140,000  cubic  feet 
=  lO.Oppt. 

=  8.0  ppt 

=  1 .5  ppt 


Result 

The  volume  (Vj)  will  be  reduced  by  the  submerged  volume  of  ship  (V3) 
Volume  of  water  in  LI  below  msl  with  ship  inside  (V4)  =  (Vj)  -  (V3) 

=  6,011,720-  1,140,000 
=  4,871,720  cubic  feet. 


Operation  Ll/2 

Raise  water  level  in  LI  by  26  feet  by  releasing  water  from  L2  to  LI 

Result 

Volume  of  water  added  (V5)  =  Lift  volume  of  LI 

=  110  x  1051  x  26 
=  3,005,860  cubic  feet 
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Entire  volume  of  water  in  LI  participates  in  mixing  with  water  from  L2 
because  the  inflow  is  from  the  bottom  of  lock. 

Total  volume  in  LI  =  (V4)  +  (V5) 

=  4,871,720 
+  3.005.860 
7,877,580 

Salinity  of  mixture  in  L1(S4)  =  [  (V4)  x  (S2)  ]  +  [  (V5)  x  (S3)  ]  / 

[(V4)  +  (V5)] 

=  14.87 1.720Y81  +  G.005.860¥1.5) 
(7,877,580) 

-  38-973.760  +  4.508.790 
(7,877,580) 

=  5.52  ppt 

Exchange  coefficient  =  (5.52/10.0)  =  0.552 

Computations  for  Lock  -  L2 

Data 

Size  of  L2 
Lift 

Bed  Elevation  below  MSL 
Water  depth  at  +26  ft 
Volume  of  water  in  L2 

at  +26  feet  (V 7) 


Lift  Volume  of  L2  (V9) 

Salinity  (S3) 

Initial  Status: 

After  releasing  water  from  L2  to  LI,  water  level  in  L2  drops  from  +54  to 
+26  feet  Water  elevation  in  LI  =  +26’ 

Salinity  in  LI  =  5.52  ppt 
Salinity  in  L2  =  1 .5  ppt 
Ship  is  in  LI 

Operation  L2/1 

Open  gates  between  L2  and  LI 
Move  ship  from  LI  to  L2 
Assume  100%  mixing 

Comments 

The  volume  of  water  from  LI  that  will  mix  with  water  in  L2  is  the 
volume  above  sill  on  the  L2  side  minus  the  submerged  ship  volume 


1071  feet  x  1 10  feet 
28  feet 
-  20.33  feet 
46.33  feet 

=  110  x  1071  x  46.33 
=  5,458,137  cubic  feet 

=  110x1071  x  28 
=  3,298,680  cubic  feet 
=  1.5  ppt  (assumed) 
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=  (V2)-(V3) 
=  5,121,523 
-  1,140,000 


(V6) 


3,981,523  cubic  feet 

This  volume  will  mix  with  the  entire  volume  of  water  in  L2  above  bed 

(V7). 

Result 


Salinity  of  mixture  in  L2:  (S6)  = 

[  (V6)  x  (S4)  ]  +  [(V7)  x  (S3)  ]  /  [(V6)  +  (V7)  ] 
=  (3,981.523  x  5.52)  +  (5.458.137  x  1 .51 
3,981,523  +  5,458,137 


=  21.978.007  +  8.187.205 
9,439,660 

=  30.165.212 
9,439,660 

=  3.2  ppt 


The  measured  salinity  (S3)  =  1.5  ppt 

Hence  mixing  =  47%,  i.e.  coefficient  of  mixing  =  0.47 

Status  at  the  end  of  operation  L2/1 
Ship  is  in  L2 

Volume  of  water  in  L2  at  +26  ft  without  ship  (V7)  =  5,458,137 

Subtract  ship  volume  (V3)  1,140.000 

Volume  of  water  in  L2  at  +26  feet 

with  ship  inside  (V8)  4,3 18,137 

Salinity  of  L2  (S3)  =1.5  ppt 

Operation  L2/2 

Close  gates  between  L2  and  LI 

Raise  water  level  in  L2  by  28  feet  by  releasing  water  from 
Miraflores  Lake 

Result 

Salinity  of  mixture:  (S7) 

=  [(V8)  x  (S3)  ]  +  [(V9)  x  (S5)  ]  /  [(V8)  +  (V9)  ] 

=  (4,318.137  x  1.51  +  (3.298.680  x  1  01 
(7,616,817) 
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=  6.477.205  +  3.298.680 
7,616,817 

=  1.28  ppt 

Status  after  completion  of  operation  L2/2 

Water  elevation  in  L2  and  Miraflores  Lake  =  +  54  feet 

Volume  of  water  in  L2  at  +54  ft  without  ship  (Vjo) 

=  110  x  1071  x  74.33 
=  8,756,817 

Subtract  ship  volume  (V3)  =  1.140.000 

Volume  of  water  in  L2  with  ship  (Vj  i)  7,616,817  eft. 

Operation  L2/3 

Open  gates  between  L2  and  Miraflores  Lake 
Move  ship  from  L2  into  Miraflores  Lake 


Result 

Volume  of  water  in  Miraflores  Lake  between  +53  ft  and  +55  ft  =  87  x 

106  eft. 

Assume  total  volume  =  870  x  106  cubic  feet. 

Assume  that  100%  mixing  with  lock  water  takes  place  with  this  volume 

of  lake. 

(Instead  of  assuming  that  only  the  volume  of  water  above  sill  mixes  with 
lake  water.) 

Salinity  of  mixture  in  L2:  (S8) 

=  f870x  106x  1.01 +  (7.616.817  x  1.281 
(877.6  x  106) 

=  879.76  x  106 
877.6  x  106 

=  1.003  ppt 

Also,  note  that  26  x  106  gallons  of  fresh  water  is  added  per  ship  transit,  from 
Pacific  side,  which  is  equal  to  3.47  x  106  eft.  This  quantity  with  approximately 
zero  salinity  is  about  50%  of  the  magnitude  of  the  lock  water  volume  with  a 
salinity  of  1.28  ppt.  This  helps  in  reducing  the  salinity  of  Miraflores  Lake  water. 


Computations  for  Lock  -  L3 

Status 

Ship  is  in  Miraflores  Lake 

Water  elevation  in  L3  and  Miraflores  =  +54  ft 
Salinity  in  Miraflores  =1.0  ppt 
Salinity  in  L3  =  0 
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Operation  L3/1 

Open  gates  between  L3  and  Miraflores. 
Move  ship  from  Miraflores  to  L3. 


Result 

water  elevation:  +  54  ft  msl 
lock  bottom  +1 1 

water  depth  =  (54-11)  =  43  ft 

Volume  of  water  in  L3  at  +54  ft  without  ship:  (V12) 

=  110  x  1071  x  43 
=  5,056,830  eft 

subtract  ship  volume  (V3)  =  1.140.000 

Volume  of  water  in  L3  at  +54  ft  with  ship  inside  (V13)  3,916,830  eft 

If  we  assume  that  half  of  this  volume  mixes  fully  with  Miraflores  Lake  water 
with  salinity  of  1 .0  ppt,  then  salinity  in  L3  will  be  s9  =  0.25  ppt. 

Operation  L3/2 

Close  gate  between  L3  and  Miraflores  Lake 
Raise  water  level  in  L3  by  3 1  ft: 


Result 

Fresh  water  from  Gatun  Lake  added  to  L3  for  raising  water  level  by  3 1  ft 

=  110x1071x31 
=  3,652,1 10  eft 

Salinity  of  mixture  in  L3  after  raising  ship  to  +85  ft  (S10) 

=  1Y3.652.100)  x  (Oil  +  3.916-830  0.25] 

7,577,940 

=  0.13  ppt 

Summary.  Estimates  of  salinities  in  three  locks  (LI,  L2  and  L3)  resulting  from 
ship  transiting  from  Pacific  Ocean  are  given  in  Table  A5.  Estimated  transfer 
coefficients  are  given  in  Table  A6. 
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Table  A5 

Estimates  of  salinities  resulting  from  ships  transil 
Pacific  Entrance 

iing  from 

Ship  Transit 

Location  and  stage  of  transit 

Computed 

Salinity  (ppt) 

From  Pacific  Ocean 

To  Lock  1 

Inside  LI 

After  full  lift 

5.52 

From  Lock  1 

To  Lock  2 

inside  L2 

After  moving  ship  from  LI  to  L2 

1.50 

Inside  L2 

After  full  lift 

1.28 

From  Lock  2 

To  Miraflores  Lake 

Inside  L2 

After  moving  ship  from  L2  to  Miraflores  Lake 

1.003 

From  Miraflores  Lake 

To  Lock  3 

Inside  L3 

After  moving  ship  from  Miraflores  Lake  to  L3 

0.25 

Inside  L3 

After  full  lift 

0.13 

Table  A6 

Estimate  of  transfer  coefficients  for  salt  loading  from  Pacific 
Entrance 

Location  and  stage 

Salinity  (ppt) 

Transfer  Coefficient  (%) 

Pacific  Ocean 

35.0 

100 

Pacific  Ocean 
just  outside  LI 

10.0 

28.571 

Inside  Lock  LI 

After  ship  enters 

8.0 

22.857 

Inside  Lock  LI 

After  full  lift 

5.52 

15.771 

inside  Lock  L2 

After  ship  enters 

1.5 

4.286 

Inside  Lock  L2 

After  full  lift 

1.28 

3.657 

Inside  Lock  L2 

After  ship  leaves 

1.003 

2.866 

Miraflores  Lake 

1.00 

2.857 

Inside  Lock  L3 

After  ship  enters 

0.25 

0.71 

Inside  Lock  L3 

After  full  lift 

0.13 

0.37  : 

Miraflores  Lake  Salt  Loading  Estimate 


Salt  loading  per  transit  (set  1) 

Inflow 

:  Volume  of  water  in  L2  at  +54  ft  with  ship  inside  (VI 1) 

=  7,616,817  (with  ship  inside) 
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Salinity  inside  L2  after  full  lift  (S7)  =  1 .28  ppt. 


Assume  that  50%  of  this  mixes  with  lake  water. 


Weight  of  salt 


=  (7,616,817)  (62.6  Ib/cft)  x  0.28/2)  lb 

1000  lb 


=  304,185  lbs. 


Outflow 

When  the  water  level  in  L2  is  lowered  by  28  ft,  salt  moves  out. 
Decrease  in  volume  of  lake  water 
=  110  x  1071  x  28 
=  3,298,680  eft 


Lake  salinity  =  1 .0  ppt 


Weight  of  salt  out 


=  (3,298,680)  (62.4)  xLO 

1000 


=  205,838  lb. 


Net  salt  added  per  transit  =  304, 1 85  lb 

-  205,838  lb 


98,  347  lb  per  transit 


Salt  loading  per  transit  (set  2) 

Assume  the  same  inflow  as  in  setl . 

Outflow 

Assume  that  the  local  salinity  of  Miraflores  Lake  just  outside  L2  was  1.15  ppt 
(instead  of  1.0  assumed  in  setl) 

Salt  out  =  (3,298,680)  (62.4)  (1.15/1000) 

=  236,713  lb  /  transit 

Salt  input  =  304,1 85  Ib/transit 

Net  salt  added  per  transit  =  304, 1 85  lb 

236,713  lb 

67,472  lb  per  transit 
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Salt  loading  per  transit  (set  3) 


It  is  seen  that  a  small  change  in  the  magnitude  of  salinity  of  water  in 
Miraflores  form  1.0  to  1.15  changes  the  magnitude  of  net  salt  addition  per  transit 
from  98,347  lb  per  transit  to  67,472  lb  per  transit.  With  an  estimated  14,000 
transits  per  year,  this  makes  a  difference  of  432  million  pounds  per  year. 

Instead  of  assuming  50  percent  of  lock  water  salt  mixing  with  lake  water,  we 
assume  40  percent,  the  salt  input  changes  from  304,185  to  243,348  lbs  per  transit. 
The  following  results  of  sensitivity  analysis  computations  illustrate  the 
importance  of  the  accuracy  of  salinity  magnitudes. 


Table  A7 

Sensitivity  analysis  illustrating  impact  of  Miraflores  Lake  salinity  magni 

itude 

Miraflores  Lake  Salinity  Outside  L2 

Salt  Inflow  (lb) 

Per  transit 

Salt  Outflow  (lb) 

Per  transit 

Net  Salt  Input  (lb) 

Per  transit 

1.00 

243,348 

205,838 

37,510 

1.10 

243,348 

226,421 

16,927 

1.15 

243,348 

236,713 

6,635 

1.16 

243,348 

238,772 

4,576 

1.17 

243,348 

240,830 

2,518 

1.18 

243,348 

242,888 

459 

1.182 

243,348 

243,300 

48 

Computations  for  ship  transiting  from  Atlantic  Ocean 


Notation  is  given  in  Table  A8. 


Computations  for  Lock  L6 

Status 

Lock  gates  to  the  sea  are  open. 

Ship  is  waiting  outside  to  get  in. 

Sea  level  =  msl  =  P.L.D. 

Data 

Bottom  level  of  L6  =  -44.3  ’ 

Volume  of  water  in  L6  at  msl  (VI 4) 

=  110x  1051  x44.3 
=  5,121,523  cubic  feet 
Assume  salinity  inside  L6  (SI  1)  =  8  ppt 
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All 


Operation  L6/1 

Bring  the  ship  in  and  close  gates. 


Result 

Volume  of  water  in  lock  will  be  reduced  by  the  submerged  volume  of  ship 

V14  5,121,523 

V3  -  1.140.000 

Volume  of  water  in  L6  with  ship  (VI 5)  4,98 1 ,523  eft. 

Operation  L6/2 

Raise  water  level  by  28  ft  above  msl  by  using  water  from  L5. 


Result 

Volume  of  water  added  (VI 6) 

Assume  salinity  of  L5  (SI 2) 
Salinity  of  mixture  in  L6  (SI  3) 


=  110x1051  x  28 
=  3,237,080  eft 

=  1 .5  ppt 

=  (4,981.523  x  8V+  f3.237.080  x  1 J) 
8,218,603 

=  39.852.184 +  4.855.620 
8,218,603 

=  5.44  ppt. 


Computations  for  Lock  -  L5 

Status 

Water  elevation  in  L6  =  +28  ft 
Water  elevation  in  L5  =  +28  ft 
Salinity  in  L6  after  full  lift=  5.44  ppt 
Salinity  in  L5  =  1.5  ppt 

Operation  L5/1 

Open  gates  between  L6  and  L5 
Move  ship  from  L6  to  L5 


Result 

Total  water  depth  =  (28  +  15.33)  =  43.33  ft 
Volume  of  water  in  L5  without  ship  at  elevation  +28 
V17  =  110x  1051  x43.33 

=  5,009,381  eft 


For  conservative  estimate,  assume  100%  mixing,  with  entire  volume  of  water  in 
L6 
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Salinity  of  mixture 
In  L5  (SI 4) 


=  t4.981.523  x  5.441  +  (5.009.381  x  1.5) 
9,990,904 

=  27.099.485  +  7.514.071 
9,990,904 

=  3.46  ppt 


Salinity  assumed  in  L5  was  1.5 

Mixing  =  43.35%  say  40% 

Status  after  moving  ship  from  L6  to  L5 

Volume  of  water  in  L5  (V17)  =  5,009,381 

-  submerged  ship  volume  (V3)  -  1,140,000 

Volume  of  water  in  L5  with  ship  (V18)  =  3,869,381 

Operation  L5/2 

Raise  water  level  by  29  ft  from  elevation  +28  ft  to  57  ft  by  adding  water 

from  L4 
Result 

Volume  of  water  added  (VI 9)  =  1 10  x  1051  x  29 

=  3,352,690  eft 

In  the  absence  of  field  data,  salinity  of  L4  is  assumed  to  be  zero. 

Salinity  of  mixture  in  L5  (SI  5)  =  (3.869.381  x  1.5)  +  (3.352.690  x  0) 

7,222,071 


=  0.8  ppt 


Computations  for  Lock  -  L4 

Status 

Water  elevation  in  L5  =  water  elevation  in  L4  =  +57  ft 
Ship  is  in  L5 

Operation  L4/1 

Open  gates  between  L5  and  L4  and  move  the  ship  to  L4 
Bed  elevation  =  +13 
Water  depth  =  44  ft 


Result 

Volume  of  water  in  L4  above  sill  at  water  at  elevation  +57  ft 

Volume  (V20)  =  1 10  x  1071  x  44 
=  5,183,640  eft. 
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Salinity  of  water  in  L5  =  0.8  ppt. 
Assume  100%  mixing  after  gates  open. 


Volume  of  water  in  L5  with  ship  inside  (V21)  =  3,869,381  eft. 
Salinity  of  water  in  L4  =  0.15  ppt  (assumed) 

Salinity  of  mixture  in  L4:  (SI  6)  = 

(3.869,381  x  0.8)  +  (5.1 83.640  x  0.1 5) 
(3,869,381  +5,183,640) 

=  3.095.505  +  777.546 
9,053,021 

=  0.43  ppt 

It  is  assumed  that  salinity  of  water  in  L4  is  0.15  ppt. 

This  gives  [0.15  /  0.43]  =  35%  mixing 

Operation  L4/2 

Raise  water  level  in  L4  from  elevation  +57  to  +85’  (by  28  feet),  using 
water  from  Gatun  Lake. 

Result 

Volume  of  water  in  L4  without  ship  (V21)  =  5,1 83,640 

Subtract  ship  volume  .  1.140.QQQ 

Volume  of  water  in  L4  with  ship  inside  (V22)  4,043,640 

Salinity  in  L4  =  0.15  ppt. 

Volume  of  fresh  water  added 
=  110  x  1071  x  28 
=  3,298,680  eft. 


Salinity  of  mixture  in  L4:  (SI  8)  =  4,043.640  x  015 

7,342,320 


=  0.083  ppt 


Computations  for  Gatun  Lake 

Water  level  = 

Sill  level 

Difference  = 

Volume  of  salt  water  above  sill  = 

Subtract  ship  volume  = 

Weight  of  water  =  279.5  x  106  lb 
Salinity  =  0.083  ppt 


+  85  ft 
+  37.3  ft 
47.7  ft 

1 10  x  1071  x  47.7  eft 

5.619.537 
1.140,000 

4.479.537  eft 
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Weight  of  salt  =  279.5  x  1 06  x  0,083 

1000 


=  23,198  lbs. 

Assuming  50%  mixing,  weight  of  salt  added  per  transit  =  1 1 ,600  lbs. 

Summary:  Estimates  of  salinities  in  three  locks  (L4,  L5  and  L6)  resulting  from 
ship  transiting  from  Ocean  are  given  in  Table  A9.  Estimated  transfer  coefficients 
are  given  in  Table  A 10. 


Table  A8 

Notations  used  in  salinity  computations 

Notation 

Parameter 

Magnitude 

Volumes 

V14 

Volume  of  L6  below  msl 

5,121,523  eft 

V15 

Volume  of  water  in  L6  above  bed  with  ship  inside 

4,981,523  eft  | 

V16 

Lift  Volume  of  L6 

3,237,080  eft 

V17 

Volume  of  water  in  L5  at  +28  ft  without  ship  inside 

5,009,381  eft. 

V18 

Volume  of  water  in  L5  at  +28  ft  with  ship  inside 

3,869,381  eft  j 

V19 

Volume  of  water  added  to  L5  from  L4 

3,352,690  eft 

V20 

Volume  of  water  in  L4  at  elevation  +57  ft 

5,183,640  eft 

V21 

Volume  of  water  in  L4  at  +57  ft  without  ship  inside 

5,183,640  eft 

V22 

Volume  of  water  in  L4  at  +57  ft  with  ship  inside 

4,043,640  eft 

Salinity 

S11 

Assumed  salinity  inside  L6 

8.0  ppt. 

S12 

Salinity  inside  L5 

1.5  ppt 

S13 

Computed  salinity  in  L6  after  full  lift 

5.44  ppt 

S14 

Computed  salinity  in  L5  after  operation  L5/1 

3.46  ppt 

S15 

Salinity  in  L5  after  full  lift 

0.8  ppt 

S16 

Computed  salinity  in  L4  after  operation  L4/1 

0.43  ppt 

S17 

Computed  salinity  in  L4  after  operation  L4/2 

0.083  ppt 
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Table  A9 

Estimate  of  salinities  resulting  from  ships  transiting  from  Alantic 
Entrance 

Source 

Volume  of  water 
(cubic  feet) 

Location  and  Time 

Salinity 

<ppt) 

Atlantic  Ocean 

Atlantic  Ocean  just  outside  LI 

All  times. 

10.0 

Lock  L6 

Inside  L6 

All  times 

8.0 

Lock  L6 

Lock  L6 

After  full  lift 

5.44 

From  Lock  L6 

To  Lock  L5 

4,981,523 

Lock  L5 

After  exchange  with  L6 

1.50 

Lock  L5 

Lock  L5 

After  full  lift 

0.80 

From  Lock  L5 

To  Lock  4 

3,869,381 

Lock  L4 

After  moving  ship  from  L5  to  L4 

0.43 

Lock  L4 

Lock  L4 

After  full  lift 

0.083 

Gatun  Lake 

All  times 

0.0 

Table  A10 

Estimate  of  salt  transfer  coefficients  from  Atlantic  Entrance 

Location  and  stage 

Salinity  (ppt) 

Transfer  Coefficient  (%) 

Atlantic  Ocean 

35.0 

100 

Atlantic  Ocean 
just  outside  L6 

10.0 

28.6 

Inside  Lock  L6 

After  ship  enters 

8.0 

22.8 

Inside  Lock  L6 

After  full  lift 

5.44 

15.54 

Inside  Lock  L5 

After  ship  enters 

1.5 

4.29 

inside  Lock  L5 

After  full  lift 

0.8 

2.29 

Inside  Lock  L4 

After  ship  enters 

0.15 

0.43 

Inside  Lock  L4 

After  full  lift 

0.083 

0.24 
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Figure  A1 .  Schematic  plan  and  elevations  for  Panama  Canal  Locks 
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Appendix  B 
Summary  of  Saltwater 
Intrusion  Problems  at  Inland 
Navigation  Channels  in  the 
United  States 


A  literature  search  for  information  on  the  environmental  effects  of  saltwater 
intrusion  in  the  United  States  (USA)  yielded  no  bona  fide  before  and  after  studies  of 
saltwater  intrusion  due  to  inland  navigation  channels.  Likewise,  Brogdon  (1986) 
surveyed  coastal  Districts  of  the  U.S.  Army  Corps  of  Engineers  to  obtain  pre-  and 
post-channel  deepening  salinity  data,  but  the  survey  yielded  very  little  useful  data. 
Intuitively,  any  inland  navigation  channel  which  involved  significant  increases  in  the 
natural  depths  of  estuaries  or  which  opened  a  new  inlet  to  the  ocean  likely  caused 
saltwater  intrusion.  Thus,  many  inland  navigation  channels  have  probably  caused 
saltwater  intrusion,  but  neither  the  magnitude  nor  the  extent  has  been  documented. 
Navigation  channels  discussed  herein  are  limited  to  those  which  are  definitively 
known  to  have  caused  saltwater  intrusion,  based  on  knowledgeable  sources,  and 
those  which  limited  data  indicate  that  saltwater  intrusion  is  likely  to  have  occurred. 

The  Committee  on  Tidal  Hydraulics  reviewed  44  U.S.  Army  Corps  of  Engineers 
(US  ACE)  projects  in  Technical  Bulletin  No.  17.  Of  the  44,  the  following 
navigation  channels  are  definitively  stated  by  the  Committee  to  have  caused 
saltwater  intrusion  after  construction:  Moriches  Inlet,  New  York,  Mississippi  River- 
Gulf  Outlet,  Louisiana,  and  Calcasieu  River  and  Pass,  Louisiana  (USACE,  1971). 
The  Gulf  Intracoastal  Waterway  (GIWW)  has  probably  caused  saltwater  intrusion 
in  Lake  Borgne,  Louisiana  (Gagliano,  1973).  Both  the  Houma  navigation  Canal 
from  Cat  Island  pass  to  Houma,  Terrebonne  Parish,  Louisiana,  and  the  Barataria 
Bay  Waterway  in  Jefferson  Parish,  Louisiana  have  caused  saltwater  intrusion 
(Gagliano,  1973).  Although  not  specifically  stated  by  the  Committee  on  Tidal 
Hydraulics  as  the  cause  of  saltwater  intrusion,  the  Lake  Washington  Ship  Canal, 
Washington,  includes  a  saltwater  basin  upstream  of  the  main  lock  to  retard  saltwater 
intrusion  into  Lake  Washington  (USACE,  1971).  Gosselink  et  al.,  (1979)  indicate 
that  the  Sabine-Neches  Ship 
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Channel  has  caused  saltwater  intrusion  in  the  Sabine  and  Neches  Rivers,  the  ship 
channel,  and  the  north  end  of  Sabine  Lake,  Texas,  USA. 

The  case  studies  of  saltwater  intrusion  examined  by  Cook  and  McGaw  (1981) 
in  the  National  Waterways  Study  included  the  Chesapeake  and  Delaware  Bays,  the 
Gulf  Intracoastal  Waterway  (GIWW)  West  I  (New  Orleans  to  the  Calcasieu  River), 
and  San  Francisco  Bay.  Cook  and  McGaw  (1981)  concluded  that  there  have  been  ' 
no  documented  cases  of  saltwater  intrusion  due  to  existing  Corps  of  Engineers 
navigation  projects  in  the  Chesapeake  and  Delaware  Bays.  The  Calcasieu  River 
Salt  Water  Barrier  and  the  Catfish  Point  Control  Structure,  Schooner  Bayou  Control 
Structure,  Calcasieu  Lock  and  Vermillion  Lock  in  the  Mermentau  Basin  were 
constructed  to  prevent  saltwater  intrusion  from  navigation  channels.  Saltwater 
intrusion  due  to  the  Houma  Canal,  which  is  a  15-foot  (4.6m)-deep  by  150-foot 
(45.7-m)-wide  channel  providing  a  direct  route  between  the  GIWW  and  deep  water 
in  the  Gulf,  has  adversely  impacted  the  city  of  Houma=s  water  supply  since 
construction  in  1962.  The  authors  indicate  that  navigation  channels  are  one  of 
several  factors  which  influence  salinity  in  San  Francisco  Bay/Delta  (Cook  and 
McGaw,  1981);  however,  they  did  not  distinguish  the  relative  contribution  of 
navigation  channels  to  saltwater  intrusion  from  other  causes. 

Browder  and  Moore  (1981)  summarized  recent  research  on  freshwater  inflow  in 
Florida,  Louisiana  and  Texas  estuaries.  Although  saltwater  intrusion  occurs  in 
many  of  the  estuaries  mentioned  in  the  summary,  the  authors  did  not  specifically  tie 
it  to  navigation  channels.  In  the  case  of  Apalachicola  Bay,  however,  they  state  that 
tongues  of  high-salinity  water  often  extend  into  the  bay  along  the  bottom  through 
passes  (Livingston  et  al.,  1974).  (Browder  and  Moore,  1981).  The  passes  referred 
to  are  natural  deep  water  (Indian  Pass  and  West  Pass)  and  deep  water  due  to  a 
navigation  channel  (Sike=s  Cut)  (Livingston  et  al.,  1974). 

Hackney  and  Yelverton  (In  Press)  showed  that  sea  level  rise  and  deepening  and 
widening  the  Wilmington  Harbor  Federal  Navigation  Channel  in  the  Cape  Fear 
River  Estuary,  North  Carolina,  have  caused  saltwater  intrusion. 

Saltwater  intrusion  has  occurred  in  the  East  Bay  Basin,  Texas,  due  to  opening 
of  Rollover  Pass  (Gosselink  et  al.,  1979). 

In  a  study  of  salinity  in  the  Columbia  River  Estuary,  Oregon,  USA,  McConnell 
et  al.  (1981)  showed  that  during  monthly  neap  tides  (tides  with  decreased  range)  in 
September  and  October  of  1977  and  1978,  higher  salanities  extended  further 
upstream  and  lasted  longer  than  in  previous  studies.  Freshwater  inflows  during  the 
1977  study  period  were  considered  low  and  during  the  1978  study  period  were 
considered  normal  (McConnell  et  al.,  1981).  McConnell  et  al.  (1981)  indicate  that 
deepening  and  widening  the  Columbia  River  Bar  navigation  channel  in  the  spring 
and  summer  of  1977  could  have  contributed  to  the  increased  saltwater  intrusion; 
however,  the  scope  of  the  study  did  not  involve  a  determination  of  the  cause  of 
saltwater  intrusion. 
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Appendix  C 

Salinity  Intrusion  at  Lake 
Washington  Locks  and  Ship 
Canal  Project  Seattle,  WA,  USA 


David  A.  Schuldt,  P.E. 


Introduction 

The  U.  S.  Corps  of  Engineers,  Lake  Washington  Locks  and  Ship  Canal 
(LWSC)  project  is  located  in  the  metropolitan  area  of  Seattle,  Washington.  Formal 
name  for  the  lock  structures  is  Hiram  M.  Chittenden  Locks,  named  after  the  Corps= 
district  engineer  during  design  of  the  project.  Construction  was  accomplished 
between  1 9 1 1  - 1 7  to  provide  access  to  and  from  the  salt  water  regime  of  Puget 
Sound  and  the  fresh  water  lakes  of  Lake  Union  and  Lake  Washington.  At  the  time, 
purpose  of  the  project  was  to  provide  for  deep-draft  and  fishing  vessels,  log  rafts, 
and  other  tug  and  barge  commodities  access  to  and  from  commercial  facilities  in 
Lake  Union  and  Lake  Washington.  Present  day  lock  useage  is  primarily  recreational 
and  commercial  fishing  vessels,  with  tug  and  barge  transport  still  common,  along 
with  a  few  log  raft  make-ups  and  medium-draft  commercial  vessels.  Present  day 
lockages  are  highest  during  the  late  spring-early  fall  time  period  because  of  high 
recreation  and  commercial  fishing  transits.  These  high  lockages  also,  unfortunately, 
coincide  with  low  precipatation  in  the  Seattle  area,  and  the  periods  of  high 
household  water  consumption  (home,  lawn,  air  conditioning  consumption,  etc.)  and 
the  need  for  discharge  of  fresh  water  attraction  waters  for  andromenous  salmon 
migrating  to  fresh  water  spawning  areas,  and  desires  to  maintain  high  water 
elevations  along  the  lake’s  shoreline  for  recreational  vessel  docks.  Lake  levels  are 
maintained  at  their  highest  levels  during  summer  to  provide  for  lockage  water  usage, 
andromenous  fish  attraction,  and  reduction  of  saline  water  intrusion  into  the  lakes. 

In  early  fall  the  lakes  are  generally  lowered  about  two  feet  to  prevent  wind-driven 
wave  damage  along  the  lake’s  shoreline  and  storage  for  fresh  water  flood 
discharges. 
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The  1900’s  project  consisted  of  the  locks  structure,  an  13  km  (  8  mile)  long, 

9  m  (30-foot)  deep  navigation  channel,  lowering  of  Lake  Washington  about  eight 
feet,  and  relocation  of  Cedar  River  discharge  into  Lake  Washington.  Photograph 
and  project  layout  are  shown  in  Figure  1  and  general  project  features  are  shown  in 
Figure  2.  Location  map  of  the  LWSC  project,  project  features  of  the  Hiram  M. 
Chittenden  locks,  and  typical  yearly  Lake  Washington  hydrographic  elevations  are 
shown  in  Plates  1  through  5. 


Hydraulic/Physical  Components 

The  locks  and  spillway  dam  regulate  the  water  surface  elevations  of  the  lakes 
tributary  to  the  locks  structure.  The  lakes  are  normally  regulated  at  elevation  +  6. 1 0  m 
(+20  feet)  to  +  6.71  m  (+22  feet)  mean  lower  low  water  (MLLW)  Puget  Sound  datum. 
Lake  surface  areas  are  about  90  sq  km  (35  square  miles)  and  shorelines  total  about 
144  km  (90  miles).  Fresh  water  inflow  during  summer  months  is  usually  less  than  56 
cu  m/s  (2,000  cfs),  and  as  low  as  28  cu  m/s  (1 ,000  cfs)  during  drought  years.  Winter 
flow  highs  are  in  the  227  to  425  cu  m/s  (8000  to  15,000  cfs)  range. 

Puget  Sound  tides  are  of  the  mixed  type;  i.e.,  two  highs  and  lows  per  the  24  plus 
hour  tidal  day.  Mean  higher  high  water  (MHHW)  is  about  3.66  m  (+12  feet)  above 
zero  mean  lower  low  water  (MLLW);  extreme  tidal  levels  are  on  the  order  of  +  4.27 
m  (+14  ft)  MLLW  and  -  1 .22  m  (-  4  feet)  MLLW.  Thus,  lockage  elevation 
differences  can  be  as  low  of  about  1.8  m  (6  feet)  to  a  maximum  of  about  7.92  m  (26 
feet).  These  extremes  are  rare,  and  most  lockages  occur  with  a  3.05  m  (10  ft  )  to 
3.66  m  (12  feet)  range  between  Puget  Sound  and  the  fresh  water  lakes. 

As  described  in  paragraph  1 ,  fresh  water  usage  and  in-flow  are  highest  and 
lowest,  respectively,  during  the  summer  months.  On  inbound  lockages,  fresh  or 
brackish  lake  water  is  gravity  fed  into  the  lock  chamber.  Because  the  lock  is  filled 
with  Puget  Sound  salt  water  prior  to  filling,  its  salinity  is  on  the  order  of  30  parts 
per  thousand  (ppt)  by  weight.  After  the  Puget  Sound  gate  closure  and  filling  of  the 
lock  chamber  with  essentially  brackish  water,  salinity  values  in  the  chamber  are 
about  20  ppt  for  filling  conditions  that  occur  at  mid-tide  levels.  For  MLLW  filling 
conditions  salinities  are  about  15  ppt.  For  filling  conditions  at  MHHW  salinities 
are  about  25  ppt.  Brackish  salt  water  inflow  to  the  lock  chamber  is  from  the  lock 
walls  near  the  concrete  lock  floor,  and  mixing  within  the  lock  chamber  is  essentially 
100%  during  the  filling  operations.  During  opening  of  the  upstream  mitre  gates  a 
surface  inflow  of  water  into  the  lock  chamber  is  readily  apparent;  this  is  a  result  of  a 
few  inches  of  higher  lake  vs.  lock  chamber  head  difference  and  the  intrusion  of 
denser  saline  water  into  the  fresh  water  lake  while  being  replaced  by  fresh  water  lake 
surface  flow  into  the  lock  chamber. 

For  outbound  passages  to  Puget  Sound  the  lock  chamber  waters  are  semi- 
stratified  as  a  result  of  the  upper  lock  sill  being  3.96  m  (13  feet)  above  the  lock  floor 
and  by  the  salt  water  barrier,  which  is  normally  maintained  in  its  upright  position, 
and  is  10  m  (33  feet)  above  the  chamber  floor.  After  lowering  lock  chamber  waters 
to  the  Puget  Sound  tide  level,  the  lower  mitre  gate  is  opened,  vessels  exit  the 
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chamber,  and  then  the  remaining  lock  chamber  water  is  essentially  intruded  and 
replaced  by  denser  Puget  Sound  saline  waters. 


Salinity  Control  Measures 

During  the  critical  summer  low  fresh  water  inflow  period  to  Lake  Washington  a 
number  of  measures  have  been/are  employed  to  minimize  or  prevent  salt  intrusion 
into  Lake  Washington.  Operating  criteria  is  for  a  maximum  of  one  ppt  salinity  at 
the  9  m  (30-foot)  deep  Montlake  Cut  channel,  located  about  four  miles  above  the 
locks. 

Salt  water  intrusion  control  measures  include  the  following: 

a.  Relatively  high  fresh  water  discharge  out  of  Lake  Washington. 

b.  Use  of  the  small  lock  as  much  as  possible,  and  delaying  large  lockage 
use  until  sufficient  numbers  of  vessels  occupy  the  lock  chamber. 

c.  The  salt  water  scavenging  drain  just  upstream  of  the  locks. 

d.  Maintaining  the  salt  water  barrier  structure  in  its  upright  condition, 
and 

e.  Mini-flushing  of  the  lock  chamber  on  eastbound  lockages. 

These  measures  are  described  below. 


Measure  1 


During  winter,  high  flow  discharge  holds  back  intrusion  of  salt  water  and  also 
flushes  saline  waters  that  have  accumulated  in  deeper  water  sinks  above  the  locks. 
During  low  summer  flow,  salt  water  intrudes  up  the  channel  and  most  then  sinks  to 
the  deeper  waters  of  Lake  Union.  Once  the  deeper  waters  of  Lake  Union  are 
replaced  by  saline  waters  the  intrusion  continues  further  up  the  canal.  Critical  to 
project  operation  is  to  prevent  intrusion  of  salinities  greater  than  one  ppt  at  the 
Montlake  Cut  channel  and  into  the  very  deep  Lake  Washington  where  salt  would 
sink  and  essentially  accumulate  without  bottom  flushing  action.  While  lock 
operations  generally  prevent  intrusion  to  the  Montlake  Cut  channel  during  most 
years,  this  condition  has  occurred  during  severe  droughts.  Critical  to  halting  the 
intrusion  is  the  discharge  of  fresh  water  flow  with  sufficient  velocity  head  to  force 
retreat  of  the  saline  wedge.  During  these  critical  conditions,  operations  are  closely 
coordinated  with  the  Seattle  Water  Department.  One  method  to  mitigate  the 
intrusion  is  to  increase  discharge  from  the  Water  Department’s  reservoir,  which  can 
result  in  restrictions  on  municipal  use  in  the  area.  Another  method  is  to  increase 
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discharge  by  lowering  Lake  Washington’s  surface  elevation  to  provide  necessary 
outflow  through  the  Montlake  Cut  channel.  A  third  method  is  to  curtail  use  of 
lockages  until  there  are  sufficient  vessels  to  essentially  fill  the  lock  chamber.  This 
operation  saves  use  of  fresh  water  infilling  of  the  locks  and  also  reduces  the  volume 
of  salt  intrusion  above  the  locks.  These  alternatives  are  all  carefully  and  thoroughly 
coordinated  with  the  Seattle  Water  Department,  environmental  agencies,  and  the 
general  public  prior  to  any  implementation. 


Measure  2 

Use  of  the  small  lock  is  routinely  performed  as  often  as  possible  because  of  its 
relatively  small  water  use  volume  compared  to  that  of  the  large  lock.  The  small  lock 
is  9.14  m  (30-foot)  wide,  45.72  (150-foot )  long  with  a  floor  elevation  of  4.9  m  (16 
feet)  MLLW.  It  requires  about  1/30  th  of  the  volume  of  water  required  for  filling 
the  large  locks  chamber.  The  upper  sill  elevation  of  the  small  lock  is  at  +  1.22  m  (+ 
4  feet)  MLLW.  Use  of  this  lock  is  ideal  for  recreation  and  most  fishing  vessels:  salt 
water  intrusion  is  basically  insignificant  compared  to  the  large  locks  and 
filling/emptying  times  are  only  about  five  minutes  compared  to  about  15-20  minutes 
for  the  large  locks.  The  small  locks  transit  about  60  percent  of  some  75,000  vessels 
per  year  and  use  of  the  small  lock  is  significant  in  that  use  of  the  large  lock  is 
reduced;  and  the  salt  volume  intmsion,  fresh  water  filling  volume,  and  lockage  times 
are  greatly  reduced.  For  purposes  of  this  report,  salt  water  intrusion  and  fresh  water 
usage  for  the  small  lock  is  so  relatively  small  compared  to  the  large  locks  that  the 
small  locks  factors  are  not  considered  in  further  discussions.  What  must  be  kept  in 
mind,  however,  is  that  without  the  small  lock  structure,  the  large  lock  would  often  be 
required  to  transit  only  a  few  small  vessels  at  a  time.  For  purposes  of  this  report, 
salt  water  intrusion  and  fresh  water  usage  of  the  small  locks  is  not  considered 
further. 


Measure  3 

Incorporated  in  the  original  construction  of  the  locks  was  a  salt  water 
scavenging  drain.  A  610  m  (2,000-foot)  long  and  61  m  (200-foot)  wide  sump 
channel  with  a  3.72  sq  m  (40  square  foot)  intake  culvert  was  originally  constructed. 
Invert  of  the  culvert  is  at  -9. 14  m  (-30  feet)  MLLW  and  water  is  then  discharged 
into  Puget  Sound  waters  via  a  213  m  (700-foot)  long,  2.8  sq  m  (30  square  feet) 
culvert.  To  reduce  intake  velocities,  which  resulted  in  draw-in  of  mid  and  upper 
surface  less  saline  waters,  the  U.S.  Army,  Committee  on  Tidal  Hydraulics  (CTH) 
recommended  enlarging  the  culvert  intake  area  and  adding  a  cover  structure  over  the 
intake  to  increase  the  discharge  of  higher  salt  water  content  water  in  lieu  of  that 
partially  diluted  from  the  intake  of  mid  and  upper  surface  waters.  The  scavenging 
culvert  intake  was  increased  to  1 7. 1  sq  m  ( 1 84  square  feet)  in  the  mid- 1 960’s  and 
its  intake  lowered  a  few  feet;  however,  the  lid  cover  has  not  been  constructed. 
Increasing  the  culvert  intake  area  has  reduced  its  velocity  from  about  seven  to  two 
feet  per  second;  this  reduction  significantly  increasing  the  outflow  percentage  of  salt 
water  discharge. 
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As  a  fish  and  water  conservatory  measure,  a  pipe  diversion  from  the  salt  water 
drain  was  installed  and  diverted  into  the  fish  ladder  for  upstream  migrating  salmon 
attraction.  The  salt  water  drain  is  generally  open  24-hours  per  day  throughout  the 
year  and  is  the  largest  use  of  water  from  above  the  locks,  generally  surpassing  lock 
usage  by  factor  of  about  two-to-three,  about  2  to  3.7  cu  m/s  (70  -130  cfs)  for 
lockages  during  winter  and  summer,  to  about  8.5  to  9.3  cu  m/s  (300  -  330  cfs)  for 
the  salt  water  drain  during  winter  and  summer,  respectively.  Fish  water  attraction 
water  diverted  from  the  salt-water  drain  system  is  only  about  10  percent  of  the  salt 
water  scavenging  drain  system. 


Measure  4 

The  salt-water  barrier  for  the  large  lock  was  originally  constructed  in  1966.  The 
barrier  consists  of  a  46  metric  tons  (42-ton)  hollow,  welded  steel  structure,  24.4  m 
(80-foot)  wide  by  6.1  m  (20-foot)  high,  and  is  6.7  m  (2  2-foot)  thick.  The  structure 
is  divided  into  16  buoyancy  compartments  by  steel  plating  arranged  in  two 
horizontal  rows  across  the  width  of  barrier  and  vertically  into  3  m  (10-foot) 
spacings  by  steel  plates.  The  barrier  is  hinged  to  the  lock  sill  floor  by  five  bearing 
pedestals.  The  barrier  is  raised  /  lowered  by  injecting  air  into  the  chambers  and 
forcing  water  out;  or  vice  versa,  allowing  flooding  of  the  chambers  to  lower  the 
barrier.  Barrier  details  are  shown  on  Plate  4.  The  barrier  is  generally  operated  in  the 
upright  position  as  deeper-draft  vessels  or  tows  with  chain  are  now  relatively 
infrequent.  With  its  elevation  of  1 .2  m  (4  ft)  MLLW,  lowering  of  the  barrier  is 
seldom  required  as  vessel  drafts  rarely  exceed  4.5  m  (15  feet). 

When  the  upper  mitre  gates  are  opened  for  inbound  passages  the  barrier  acts  as 
a  physical  obstacle  to  intrusion  of  the  lower  depth  lock  chamber  saline  waters. 

Upon  opening  of  the  mitre  gate,  salt  water  within  the  chamber  is  generally  on  the 
order  of  20  ppt,  and  is  well-mixed  both  vertically  and  horizontally.  What  the  barrier 
does  is  essentially  minimize  the  salt  water  intrusion/fresh  water  inflow  layer  into  the 
lock  chamber  to  an  5.5  m  (18-foot)  layer  compared  to  1 1  m  (36  feet)  without  the 
barrier.  Thus,  with  the  barrier  in  an  upright  condition,  the  two-layer  interchange  of 
fresh  and  salt  water  is  decreased  by  one-half,  and  overall  should  result  in  a  50 
percent  volume  decrease  of  salt  water  intrusion  to  the  lakes.  Upon  leaving  the  lock 
chamber,  the  denser  salt  water  sinks  into  the  above-locks  deep  scavenging  channel, 
most  of  which  is  intercepted  at  the  drain  inlet  and  discharged  to  Puget  Sound  waters. 
The  denser  salt  water  that  does  escape  the  scavenging  drain  does  then  progress 
upstream  as  a  low  bottom  wedge  of  saline  waters. 

Salinity  measurements  are  taken  at  a  number  of  locations  along  the  fresh  water 
canal.  The  first  is  at  the  upstream  nose  pier  of  the  large  lock.  May  1999  surface  to 
bottom  salinities  are  shown  in  Figure  3.  Sensor  A  represents  readings  taken  0.6  m  (2 
feet)  above  bed  and  sensor  D  represents  readings  taken  0.6  m  (2  feet)  below  water 
surface.  Comparing  daily  logs  of  tidal  heights,  barrier  in  its  up  or  down  position  at 
times  of  lockage,  the  salinity  peaks  along  the  bottom  “spike”  when  the  barrier  is  in  the 
down  position,  and  to  some  extent  when  lockages  occur  at  high  tide  when  filling  of  the 
lock  chamber  requires  a  lesser  volume  of  fresh  water.  Immediately  upstream  of  the 
lock,  low  bottom  salinties  are  on  the  order  of  12-20  ppt  with  the  barrier  down.  When 
in  its  upright  position  these  salinities  are  only  on  the  order  of  5-10  ppt. 
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The  barrier  at  the  LWSC  project  has  shown  to  be  an  effective  deterrent  to  salt 
water  intrusion  of  the  upstream  fresh  water  regime.  It  has  been  operational  for  about 
30  years  and  is  an  integral  operational  method  to  reduce  salt  water  intrusion. 


Measure  5 

Mini-flushing  of  the  large  lock  chamber  was  begun  and  tested  about  1990  to 
supplement  other  measures  in  preventing  excessive  salt  water  propogation  during 
drought  years.  The  mini-flushing  is  only  employed  on  inbound  transits,  and  two 
methods  of  mini-flushing  were  tested,  both  with  the  upstream  barrier  in  its  raised 
position. 

The  first  method  is  the  “pre-lockage  mini-flushing”  and  has  been  used  for  a 
number  of  years.  In  this  method,  inbound  vessels  are  brought  into  the  lock  chamber 
and  the  Puget  Sound  mitre  gates  closed.  Salt  water  is  flushed  from  inside  the 
chamber  by  opening  both  the  above  lock  intake  and  below  lock  emptying  culverts. 
During  the  first  five  minutes,  little  salinity  decrease  was  noted  in  the  lock  chamber 
as  water  drawn  into  the  lock  was  of  relatively  high  saline  values  as  a  result  of 
previous  lockages;  and  only  after  additional  flushing,  was  a  significant  lowering  of 
lock  salinity  noted  in  the  lock  chamber.  Upon  closing  the  Puget  Sound  discharge 
gate,  raising  the  lock  chamber  waters  to  lake  levels,  and  opening  the  upbound  mitre 
gates,  only  small  reductions  of  salinity  were  observed  in  the  upstream  30-foot 
channels. 

The  second  method  tested  is  the  “post-lockage  mini-flushing”.  This  consists  of 
closing  the  Puget  Sound  mitre  gate,  raising  the  water  level  of  the  lock  chamber,  and 
then  beginning  the  “post-lockage  mini-flushing”  operation.  Once  the  chamber  is 
filled,  the  upstream  intake  valves  are  closed,  and  the  lake  mitre  gates  as  well  as  the 
Puget  Sound  culvert  discharge  valves  are  opened.  These  latter  steps  flush  salt  water 
from  the  lower  levels  of  the  lock  chamber  as  fresher  water  enters  the  chamber,  and 
significant  salt  water  decreases  in  the  chamber  resulted  during  five  minutes  of 
flushing.  Salinities  in  the  upstream  channels  were  also  significantly  reduced  by  this 
operations.  Only  two-to-five  minutes  of  flushing  with  the  Puget  Sound  valves  open 
are  required  to  significantly  reduce  salt  intrusion  above  the  locks.  Salinity  intrusion 
increased  during  curtailment  of  the  operation  in  early  August  of  1992  (a  drought 
year),  but  decreased  significantly  when  this  operational  practice  was  again 
implemented  two  weeks  later. 

The  “post-lockage  mini-flushing”  was  clearly  determined  to  be  efficient  in 
preventing  salt  intrusion  and  requiring  less  fresh  water  useage  than  the  previously 
described  “pre-lockage  mini-flushing.” 
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Field  Data 


In  their  paper  entitled  “Occurrence  and  control  of  salt  intrusion  through 
navigation  locks”,  Long  et  al,  (1957)  have  reported  field  data  on  salinity 
measurements  conducted  at  the  Chittenden  Locks  in  Lake  Washington  Ship  Canal. 
They  have  reported  on  the  removal  of  salt  water  from  the  fresh  water,  which  was 
achieved  by  flushing  through  a  sluice  at  the  lock.  The  following  conclusions  were 
drawn  from  the  field  measurements: 

a.  Flushing  through  sluice  was  remarkably  effective.  The  sluice  effluent 
consisted  of  approximately  two  parts  water  and  one  part  salt.  It  was 
necessary  to  use  fresh  water  equivalent  to  41  percent  of  the  lock  volume 
measured  between  the  upper  sill  and  pool  water  levels. 

b.  There  was  a  rational  relationship  between  the  amount  of  salt  intrusion 
through  lock  gate  and  the  time  the  gate  is  opened. 

c.  Salt  intrusion  is  a  function  of  dilution  and  degree  of  mixing  attained  in  the 
lock,  and  relative  elevations  of  floor,  sill  and  pool  surface.  These  factors 
reduce  the  salt  intrusion  potential  through  the  upper  gates  to  75  percent  of 
the  potential  at  the  lower  gates. 


Summary  and  Conclusions 

Puget  Sound  salt  water  intrusion  into  Lake  Washington  is  a  serious  problem 
during  low-flow  summer  fresh  water  discharges.  If  the  wedge  of  salt  water  does 
enter  Lake  Washington,  the  denser  water  will  sink  into  its  deep  water,  stagnating 
bottom  waters  and  depleting  its  dissolved  oxygen  content.  Determining  the 
discharge  uses  of  fresh  water  for  lockage  use,  fresh  water  requirements  for  the 
metropolitan  area  of  Seattle,  and  conservancy  of  andromenous  fish  are  the  important 
issues  related  to  prevention  of  salt  intrusion  in  Lake  Washington. 

The  critical,  low  fresh  water  inflow  period  is  closely  monitored,  both  with 
salinity  value  measurements  upstream  of  the  locks  and  available  fresh  water 
discharge  projections;  and,  in  perhaps  the  unique  case  of  the  Pacific  Northwest,  the 
status  of  salmon  production. 

The  intrusion  of  salt  water  into  the  fresh  water  regime  is  controlled  through  a 
number  of  operational  and  structural  methods.  Foremost  are  use  of  the  salt  water 
scavenging  drain,  maintaining  the  salt  water  barrier  in  the  upright  condition,  and  use 
of  the  small  lock  to  prevent  wedge  intrusion  into  Lake  Washington.  During  critical 
years,  the  potential  for  intrusion  becomes  an  intensified  coordination  process  of 
water/lock  useage  with  municipal  water  supply  authorities,  environmental  agencies 
and  the  general  public. 
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Figure  C2.  General  project  features 
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Figure  C3.  Field  data  collected  in  May  99  at  the  upstream  nose  pier  of  large  Chittenden  Lock 
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Appendix  D 

The  TABS-MD  System 


TABS-MD  is  a  collection  of  generalized  computer  programs  and  utility 
codes  integrated  into  a  numerical  modeling  system.  TABS-MD  is  capable  of 
one-,  two-,  and/or  three-dimensional  computations;  however,  only  the  one-  and 
two-dimensional  vertically  averaged  capability  will  be  discussed  in  this 
summary.  The  system  is  used  for  studying  hydrodynamics,  sedimentation,  and 
transport  problems  in  rivers,  reservoirs,  bays,  and  estuaries.  A  schematic 
representation  of  the  system  is  shown  in  Figure  Dl.  It  can  be  used  either  as  a 
stand-alone  solution  technique  or  as  a  step  in  the  hybrid  modeling  approach.  The 
basic  concept  is  to  calculate  water-surface  elevations,  current  patterns,  sediment 
erosion,  transport  and  deposition,  the  resulting  bed  surface  elevations,  and  the 
feedback  to  hydraulics.  Existing  and  proposed  geometry  can  be  analyzed  to 
determine  the  impact  on  sedimentation  of  project  designs  and  to  determine  the 
impact  of  project  designs  on  salinity  and  on  the  stream  system.  The  system  is 
described  in  detail  by  Thomas  and  McAnally  (1985). 

The  three  basic  2D  depth-averaged  components  of  the  system  are  as  follows: 

a.  “A  Two-Dimensional  Model  for  Free  Surface  Flows,”  RMA2. 

b.  “Sediment  Transport  in  Unsteady  2-Dimensional  Flows,  Horizontal 
Plane,”  STUDH. 

c.  “Two-Dimensional  Finite  Element  Program  for  Water  Quality,”  RMA4. 

RMA2  is  a  finite  element  solution  of  the  Reynolds  form  of  the  Navier-Stokes 
equations  for  turbulent  flows.  Friction  is  calculated  with  Manning’s  equation  and 
eddy  viscosity  coefficients  are  used  to  define  the  turbulent  exchanges.  A  velocity 
form  of  the  basic  equation  is  used  with  side  boundaries  treated  as  either  slip  or 
static.  The  model  has  a  marsh  porosity  option  as  well  as  the  ability  to 
automatically  perform  wetting  and  drying.  Boundary  conditions  may  be  water- 
surface  elevations,  velocities,  discharges,  or  tidal  radiation. 


Appendix  D  The  TABS-MD  System 


Dl 


Figure  D1.  TABS-MD  schematic 


The  sedimentation  model,  STUDH,  solves  the  convection-diffusion  equation 
with  bed  source-sink  terms.  These  terms  are  structured  for  either  sand  or 
cohesive  sediments.  The  Ackers-White  (1973)  procedure  is  used  to  calculate  a 
sediment  transport  potential  for  the  sands  from  which  the  actual  transport  is 
calculated  based  on  availablity.  Clay  erosion  is  based  on  work  by  Partheniades 
(1962)  and  Ariathurai  and  the  deposition  of  clay  utilized  Krone’s  equationns 
(Ariathurai,  MacArthur,  and  Krone  1977).  Deposited  material  forms  layers  and 
bookkeeping  allows  up  to  10  layers  at  each  node  for  maintaining  separate 
material  types,  deposit  thickness,  and  age.  The  code  uses  the  same  mesh  as 
RMA2. 

Salinity  calculations,  RMA4,  are  made  with  a  form  of  the  convective- 
diffusion  equation  which  has  general  source-sink  terms.  Up  to  six  conservative 
substances  or  substances  requiring  a  decay  term  can  be  routed.  The  code  uses  the 
same  mesh  as  RMA2.  The  model  accommodates  a  mixing  zone  outside  of  the 
model  boundaries  for  estimation  of  retrainment. 

Each  of  these  generalized  computer  codes  can  be  used  as  a  stand-alone 
program,  but  to  facilitate  the  preparation  of  input  data  and  to  aid  in  analyzing 
results,  a  family  of  utility  programs  was  developed  for  the  following  purposes: 

a.  Digitizing 

b.  Mesh  generation 

c.  Spatial  data  management 

d.  Graphical  output 

e.  Output  analysis 

f  File  management 

g.  Interfaces 

h.  Job  control  language 
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Finite  Element  Modeling 


The  TABS-MD  numerical  models  used  in  this  effort  employ  the  finite 
element  method  to  solve  the  governing  equations.  To  help  those  who  are 
unfamiliar  with  the  method  to  better  understand  the  system,  a  brief  description  of 
the  method  is  given  here. 

The  finite  element  method  approximates  a  solution  to  governing  equations  by 
dividing  the  area  of  interest  into  smaller  subareas,  which  are  called  elements. 

The  dependent  variables  (e.g.,  water-surface  elevations  and  sediment 
concentrations)  are  approximated  over  each  element  by  continuous  functions 
which  interpolate  based  on  unknown  point  (node)  values  of  the  variables.  An 
error,  defined  as  the  deviation  of  the  governing  equations  using  the  approximate 
solution  from  the  equation  using  the  correct  solution,  is  minimized.  Then,  when 
boundary  conditions  are  imposed,  a  set  of  solvable  simultaneous  equations  is 
created.  The  solution  is  continuous  over  the  area  of  interest. 

In  one-dimensional  problems,  elements  are  line  segments.  In  two- 
dimensional  problems,  the  elements  are  polygons,  usually  either  triangles  or 
quadrilaterals.  Nodes  are  located  on  the  edges  of  elements  and  occasionally 
inside  the  elements.  The  interpolating  functions  may  be  linear  or  higher  order 
polynomials.  Figure  D2  illustrates  a  quadrilateral  element  with  eight  nodes  and  a 
linear  solution  surface  where  F  is  the  interpolating  function. 

Most  water  resource  applications  of  the  finite  element  method  use  the 
Galerkin  method  of  weighted  residuals  to  minimize  error.  In  this  method  the 
residual,  the  local  error  in  the  equations  use  of  the  approximate  and  solution,  is 
weighted  by  a  function  that  is  identical  to  the  interpolating  function  and  then 
minimized.  Minimization  results  in  a  set  of  simultaneous  equations  in  terms  of 
nodal  values  of  the  dependent  variable  (e.g.  water-surface  elevations  or  sediment 
concentration).  The  time  portion  of  time-dependent  problems  can  be  solved  by 
the  finite  element  method,  but  it  is  generally  more  efficient  to  express  derivatives 
with  respect  to  time  in  finite  difference  form. 


The  Hydrodynamic  Model,  RMA2 


Applications 

This  program  is  designed  for  far-field  problems  in  which  vertical 
accelerations  are  negligible  and  the  velocity  vectors  at  a  node  generally  point  in 
the  same  directions  over  the  entire  depth  of  the  water  column  at  any  instant  of 
time.  It  expects  a  vertically  homogeneous  fluid  with  a  free  surface.  The  model 
will  define  the  response  to  a  specified  horizontally  inhomogeneous  fluid.  Both 
steady  and  unsteady  state  problems  can  be  analyzed.  A  surface  wind  stress  can 
be  imposed  and  the  effects  of  the  earth’s  rotation  can  be  included. 
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Figure  D2.  Two-dimensional  finite  element  mesh 


The  program  has  been  applied  to  calculate  water  levels  and  flow  distribution 
around  islands;  flow  at  bridges  having  one  or  more  relief  openings,  in  contracting 
and  expanding  reaches,  into  and  out  of  off-channel  hydropower  plants,  at  river 
junctions,  and  into  and  out  of  pumping  plant  channels;  circulation  and  transport 
in  waterbodies  with  wetlands;  and  general  water  levels  and  flow  patterns  in 
rivers,  reservoirs,  and  estuaries. 


Limitations 

This  program  is  not  designed  for  near- field  problems  where  flowstructure 
interactions  (such  as  vortices,  vibrations,  or  vertical  accelerations)  are  of  interest. 
Areas  of  vertically  stratified  flow  are  beyond  this  program’s  capability  unless  it  is 
used  in  a  hybrid  modeling  approach.  It  is  two-dimensional  in  the  horizontal 
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plane,  and  zones  where  the  bottom  current  is  in  a  different  direction  from  the 
surface  current  must  be  analyzed  with  considerable  subjective  judgment.  It  is  a 
ffee-surface  calculation  for  subcritical  flow  problems. 


Governing  equations 

The  generalized  computer  program  RMA2  solves  the  depth-integrated 
equations  of  fluid  mass  and  momentum  conservation  in  two  horizontal  directions. 
The  form  of  the  solved  equations  is 
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where 

h 

depth 

u,v 

— 

velocities 

x,y7t 

— 

Cartesian  coordinates  and  time 

p 

= 

density  of  fluid 

e 

= 

eddy  viscosity  coefficient,  for  xx  =  normal  direction  on  x-axis 

g 

surface;  yy=  normal  direction  on  y-axis  surface;  xy  and  yx  = 
shear  direction  on  each  surface 
acceleration  due  to  gravity 

a 

= 

elevation  of  bottom 

n 

— 

Manning’s  n  value 

I486 

conversion  from  SI  (metric)  to  non-SI  units 

C 

empirical  wind  shear  coefficient 

va 

= 

wind  speed 

$ 

= 

wind  direction 

O) 

= 

rate  of  earth’s  angular  rotation 

4> 

= 

local  latitude 
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Equations  Dl,  D2,  and  D3  are  solved  by  the  finite  element  method  using 
Galerkin  weighted  residuals.  The  elements  may  be  one-dimensional  lines  or  two- 
dimensional  quadrilaterals  or  triangles  and  may  have  curved  (parabolic)  sides. 
The  shape  functions  are  quadratic  for  velocity  and  linear  for  depth.  Integration  in 
space  is  performed  by  Gaussian  integration.  Derivatives  in  time  are  replaced  by 
a  nonlinear  finite  difference  approximation.  Variables  are  assumed  to  vary  over 
each  time  interval  in  the  form 

fit)  =  /( 0)  +  at  +  btc  to<  t<  t  (D4) 

which  is  differentiated  with  respect  to  time,  and  cast  in  finite  difference  form. 
Letters  a,  b,  and  c  are  constants.  It  has  been  found  by  experiment  that  the  best 
value  for  c  is  1.5  (Norton  and  King  1 977). 

The  solution  is  fully  implicit  and  the  set  of  simulataneous  equations  is  solved 
by  Newton-Raphson  non  lineaarr  iteration.  The  computer  code  executes  the 
solution  by  means  of  a  front-type  solver  that  assembles  a  portion  of  the  matrix 
and  solves  it  before  assembling  the  next  portion  of  the  matrix.  Hie  front  solver’s 
efficiency  is  largely  independent  of  bandwidth  and  thus  does  not  require  as  much 
care  in  formation  of  the  computational  mesh  as  do  earlier  traditional  solvers. 

The  code  RMA2  is  based  on  the  earlier  versions  (Norton  and  King  1977)  but 
differs  in  several  ways.  It  is  formulated  in  terms  of  velocity  (v)  instead  of  unit 
discharge  (vh),  which  improves  some  aspects  of  the  code’s  behavior;  it  permits 
drying  and  wetting  of  areas  within  the  grid;  it  permits  specification  of  turbulent 
coefficients  in  directions  other  than  along  the  x-  and  z-axes;  it  accommodates  the 
specifications  of  hydraulic  control  structures  in  the  network;  it  permits  wetlands 
to  be  simulated  as  either  totally  wet/diy  or  as  gradually  changing  wetting  and  it 
permits  input  in  either  English  or  system  international  units.  For  a  more 
complete  description,  see  Appendix  F  of  Thomas  and  McAnally  (1985). 

The  Sediment  Transport  Model,  STUDH 


Applications 

STUDH  can  be  applied  to  clay  and/or  sand  bed  sediments  where  flow 
velocities  can  be  considered  two-dimensional  (i.e.,  the  speed  and  direction  can  be 
satisfactorily  represented  as  a  depth-averaged  velocity).  It  is  useful  for  both 
deposition  and  erosion  studies  and,  to  a  limited  extent,  for  stream  width  studies. 
The  program  treats  two  categories  of  sediment:  noncohesive,  which  is  referred  to 
as  sand  here,  and  cohesive,  which  is  referred  to  as  clay. 


Limitations 


Both  clay  and  sand  may  be  analyzed,  but  the  model  considers  a  single, 
effective  grain  size  for  each  and  treats  each  separately.  Fall  velocity  must  be 
prescribed  along  with  the  water-surface  elevations,  x-velocity,  y-velocity, 
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diffusion  coefficients,  bed  density,  critical  shear  stresses  for  erosion,  erosion  rate 
constants,  and  critical  shear  stress  for  deposition. 

I 

The  program  does  not  compute  water-surface  elevations  or  velocities; 
therefore  these  data  must  be  provided.  For  complicated  geometries,  the 
numerical  model  for  hydrodynamic  computations,  RMA2,  is  used.  However, 
STUDH  can  only  accept  a  two-dimensional  network. 


Governing  equations 

The  generalized  computer  program  STUDH  solves  the  depth-integrated 
convection-dispersion  equation  in  two  horizontal  dimensions  for  a  single 
sediment  constituent.  For  a  more  complete  description,  see  Appendix  G  of 
Thomas  and  McAnally  (1985).  The  form  of  the  solved  equation  is 


where 


dC  dC  dC 

—  +  u —  +  v — 
8t  dx  dy 
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C 

u 

D* 

Dy 

a 2 


concentration  of  sediment 

depth-integrated  velocity  in  x-direction 

dispersion  coefficient  in  x-direction 

dispersion  coefficient  in  y-direction 

coefficient  of  concentration-dependent  source/sink  term 

coefficient  of  source/sink  term 


The  source/sink  terms  in  Equation  D5  are  computed  in  routines  that  treat  the 
interaction  of  the  flow  and  the  bed.  Separate  sections  of  the  code  handle 
computations  for  clay  bed  and  sand  bed  problems. 


Sand  transport 

The  source/sink  terms  are  evaluated  by  first  computing  a  potential  sand 
transport  capacity  for  the  specified  flow  conditions,  comparing  that  capacity  with 
the  amount  of  sand  actually  being  transported,  and  then  eroding  from  or 
depositing  to  the  bed  at  a  rate  that  would  approach  the  equilibrium  value  after 
sufficient  elapsed  time. 

The  potential  sand  transport  capacity  in  the  model  is  computed  by  the 
method  of  Ackers  and  White  (1973),  which  uses  a  transport  power  (work  rate) 
approach.  It  has  been  shown  to  provide  superior  results  for  transport  under 
steady-flow  conditions  (White,  Milli,  nad  Crabbe  1975)  and  for  combined  waves 
and  currents  (Swart  1976).  Flume  tests  at  the  U.S.  Army  Engineer  Research  and 
Development  Center,  Waterways  Experiment  Station  have  shown  that  the 
concept  is  valid  for  transport  by  estuarine  currents. 

The  total  load  transport  function  of  Ackers  and  White  is  based  upon  a 
dimensionless  grain  size 
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(D6) 


where 

D  =  sediment  particle  diameter 
5  =  specific  gravity  of  the  sediment 
v  =  kinematic  viscosity  of  the  fluid 

and  a  sediment  mobility  parameter 


Fgr  =  |  - - - -  | 

L  PgD(J-l)  J 

where 

r  =  total  boundary  shear  stress  =  pgRS 
where 


(D7) 


R  =  hydraulic  radius 
S  =  slope  of  water  surface 

n  =  a  coefficient  expressing  the  relative  importance  of  bed-load  and 
suspended-load  transport,  given  in  Equation  D9 
NOTE: 

ft  =  1  for  fine  sediments 

ft  =  0  for  coarse  sediments 

r  =  boundary  surface  shear  stress 

The  surface  shear  stress  is  that  part  of  the  total  shear  stress  which  is  due  to  the 
rough  surface  of  the  bed  only,  i.e.,  not  including  that  part  due  to  bed  forms  and 
geometry.  It  therefore  corresponds  to  that  shear  stress  that  the  flow  would  exert 
on  a  plane  bed. 

The  total  sediment  transport  is  (in  Kg/m3)  expressed  as  an  effective 
concentration 


where  U is  the  average  flow  speed,  and  for  l<Dgr<60 

ft’  =  1.00  -  0.56  log  Dgr 


(D8) 


(D9) 
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0.23 


A  = 


sj  Dgr 


0.14 


(DIO) 


log  C  =  2.86  log  Dgr  —  ( log  Dgr )  —  3.53  (Dll) 


For  Dgr<60 


.  »•«  +  1.34 

(D1 2) 

Dgr 

n  ’  =  0.00 

(D13) 

A  =  0.17 

(D14) 

C  =  0.025 

(D15) 

m  =  1.5 

(D16) 

Note  the  Ca  has  units  consistent  with  Gp  (kg/m3  for  STUDH). 

Equations  D6-D1 6  result  in  a  potential  sediment  concentration!  Gp.  This 
value  is  the  depth-averaged  concentration  of  sediment  that  will  occur  if  an 
equilibrium  transport  rate  is  reached  with  a  nonlimited  supply  of  sediment.  The 
rate  of  sediment  deposition  (or  erosion)  is  then  computed  as 


R  = 


(D17) 


where 

C  =  present  sediment  concentration 
tc  =  time  constant 

For  deposition,  the  time  constant  is 


tc  =  larger  of 


At 

or 

CjH 

Vs 


l 


(D18) 


and  for  erosion  it  is 
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to  =  larger  of  \  (D19) 

I  U 

{ 

where 

At  =  computational  time-step 

Cd  =  response  time  coefficient  for  deposition 

Vs  =  sediment  settling  velocity 

Ce  =  respoonse  time  coefficient  for  erosion 

The  sand  bed  has  a  specified  initial  thickness  which  limits  the  amount  of  erosion 
to  that  thickness. 


Cohesive  sediments  transport 

Cohesive  sediments  (usually  clays  and  some  silts)  are  considered  to  be 
depositional  if  the  bed  shear  stress  exerted  by  the  flow  is  less  than  a  critical  value 
xd.  When  that  value  occurs,  the  deposition  rate  is  given  by  Krone’s  (1962) 
equation 


(D20) 
(D21) 

where 

S  =  source  term 

K  =  fall  velocity  of  a  sediment  particle 
h  =  flow  depth 

C  —  sediment  concentration  in  water  column 
r  =  bed  shear  stress 
Td  =  critical  shear  stress  for  deposition 
Cc  =  critical  concentration  =  300  mg/C 
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h 

2  Vs 


1  - 


[  hcT 


C5/3 


Xd  ) 

1  - 


for  C<Cc 


Xd 
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If  the  bed  shear  stress  is  greater  than  the  critical  value  for  particle  erosion  ze, 
material  is  removed  from  the  bed.  The  source  term  is  then  computed  by 
Ariathurai’s  (Ariathurai,  MacArthur,  and  Krone  1977)  adaptation  of 
Partheniades’  (1962)  findings: 
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P 
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where  P  is  the  erosion  rate  constant,  unless  the  shear  stress  is  also  greater  than 
the  critical  value  for  mass  erosion.  When  this  value  is  exceeded,  mass  failure  of 
a  sediment  layer  occurs  and 


TlPl 

hAt 
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where 


Tl  =  thickness  of  the  failed  layer 
pL  =  density  of  the  failed  layer 
A,  =  time  interval  over  which  failure  occurs 
ts  =  bulk  shear  strength  of  the  layer 

The  cohesive  sediment  bed  consists  of  1  to  10  layers,  each  with  a  distinct 
density  and  erosion  resistance.  The  layers  consolidate  with  overburden  and  time. 


Bed  shear  stress 

Bed  shear  stresses  are  calculated  from  the  flow  speed  according  to  one  of 
four  optional  equations:  the  smooth-wall  log  velocity  profile  or  Manning 
equation  for  flows  alone;  and  a  smooth  bed  or  rippled  bed  equation  for  combined 
currents  and  wind  waves.  Shear  stresses  are  calculated  using  the  shear  velocity 
concept  where 
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where 

vb  =  bed  shear  stress 
w*  =  shear  velocity 

and  the  shear  velocity  is  calculated  by  one  of  four  methods: 

a.  Smooth-wall  log  velocity  profiles 

—  =  5.75  log  f  3.23—1  (D25) 

u*  \  v  J 

which  is  applicable  to  the  lower  15  percent  of  the  boundary  layer  when 
u*  h 

- >30 

v 

where  0  is  the  mean  flow  velocity  (resultant  of  u  and  v  components) 

b.  The  Manning  shear  stress  equation 
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where  CME  is  a  coefficient  of  1  for  SI  (metric)  units  and  1.486  for  non-SI 
units  of  measurement. 

c.  A  jonsson-type  equation  for  surface  shear  stress  (plane  beds)  caused  by 
waves  and  currents 
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where 


fw-  shear  stress  coefficient  for  waves 
Mom  maximum  orbital  velocity  of  waves 
fc  =  shear  stress  coefficient  for  currents 

d.  A  Bijker-type  equation  for  total  shear  stress  caused  by  waves  and  current 

IT  -2  1 

“* =  iifcU  +  4fwU™  (D28) 


Solution  method 

Equation  D5  is  solved  by  the  finite  element  method  using  Galerkin  weighted 
residuals.  Like  RMA2,  which  uses  the  same  general  solution  technique,  elements 
are  quadrilateral  and  may  have  parabolic  sides.  Shape  functions  are  quadratic. 
Integration  in  space  is  Gaussian.  Time-stepping  is  performed  by  a  Crank- 
Nicholson  approach  with  a  weighting  factor  (0)  of  0.66.  A  front-type  solver 
similar  to  that  in  RMA2  is  used  to  solve  the  simultaneous  equations. 


The  Water  Quality  Transport  Model,  RMA4 


Applications 

The  water  quality  model,  RMA4,  is  designed  to  simulate  the  depth-average 
advection-diffusion  process  in  most  water  bodies  with  a  free  surface.  The  model 
is  used  for  investigating  the  physical  processes  of  migration  and  mixing  of  a 
soluble  substance  in  reservoirs,  rivers,  bays,  estuarines  and  coastal  zones.  The 
model  is  useful  for  evaluation  of  the  basic  processes  or  for  defining  the 
effectiveness  of  remedial  measures.  For  complex  geometries  the  model  utilizes 
the  depth-averaged  hydrodynamics  form  RMA2. 
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The  water  quality  model  has  been  applied  to  define  the  horizontal  salinity 
distribution;  to  trace  temperature  effects  from  power  plants;  to  calculate 
residence  times  of  harbors  or  basins;  to  optimize  the  placement  of  outfalls;  to 
identify  potential  critical  areas  for  oil  spills  or  other  pollutants  spread;  to  evaluate 
turbidity  plume  extent;  and  to  monitor  other  water  quality  criterion  within  game 
and  fish  habitats. 

Limitations 

The  formulation  of  RMA4  is  limited  to  one-dimensional  (cross-sectionally 
averaged)  and  two-dimensional  (depth-averaged)  situations  in  which  the 
concentration  is  fairly  well-mixed  in  the  vertical.  It  will  not  provide  accurate 
concentrations  for  stratified  situations  in  which  the  constituent  concentration 
influences  the  density  of  the  fluid.  In  addition,  the  accuracy  of  the  transport 
model  is  dependent  on  the  accuracy  of  the  hydrodynamics  (e.q.  as  supplied  from 
RMA2). 


Governing  equations 

The  CEERD  version  of  RMA4  is  a  revised  version  of  RMA4  as  developed 
by  King  (1989).  The  generalized  computer  program  solves  the  depth-integrated 
equations  of  the  transport  and  mixing  process.  The  form  of  the  equations  solved 
is: 
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where 

h 

c 

t 

14,  V, 

Ac  Dy 

k 

a 


water  depth 

constituent  concentration 
time 

velocity  components 
turbulent  mixing  coefficients 
first  order  decay 
source/sink  of  constituent 


Note  that  the  basic  governing  equation  for  RMA4  is  the  same  as  for  the  sediment 
transport  model,  STUDH.  The  differences  between  the  two  models  lies  in  the 
source/sink  terms. 


Equation  D29  is  solved  by  the  finite  element  method  using  Galerkin 
weighted  residuals.  As  with  the  hydrodynamic  model,  RMA2,  the  transport 
model  RMA4  handles  one-dimensional  segments  or  two-dimensional 
quadrilaterals  or  triangles  with  the  option  for  curved  sides.  Spatial  integration  of 
the  equation  is  performed  by  Gaussian  techniques  and  the  temporal  variations  are 
handled  by  nonlinear  finite  differences,  consistent  with  the  method  described  in 
paragraph  15  above.  The  frontal  solution  method  is  also  used  in  RMA4,  as  with 
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the  other  programs  in  the  TABS-MD  system,  to  provide  an  efficient  solution 
algorithm. 

The  boundary  conditions  for  RMA4  are  specified  in  several  optional  ways. 
The  boundary  concentration  may  be  specified  absolutely  at  a  certain  level 
regardless  of  the  flow  direction;  the  concentration  can  be  specified  to  be  applied 
only  when  the  water  is  leaving  the  model;  or  a  mixing  zone  may  be  specified  just 
beyond  the  model  boundary  to  provide  the  possibility  of  reentertainment  of 
constituent  into  the  model  that  may  have  crossed  the  boundary  earlier.  For  a 
more  detailed  description  of  the  constituent  transport  model,  RMA4,  see  King 
and  Rachiele,  1989. 

Within  the  one-dimensional  formulation  of  the  model,  there  is  a  provision  for 
defining  the  constituent  concentration  mixing  and  transport  at  control  structures 
as  they  may  have  been  specified  in  RMA2.  These  allow  for  either  a  flow  through 
condition,  as  for  example  for  a  wier  type  flow,  or  for  a  mixing  chamber  type  of 
flux,  which  would  be  appropriate  for  a  navigation  lock. 
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An  analysis  of  results  showed  that:  a.  Freshwater  consumption  can  be  reduced  by  up  to  about  60  percent  by 
recirculating  lift  water  through  holding  ponds,  albeit  at  the  cost  of  increasing  salinity  in  the  ocean  approach 
channels  and  ultimately  in  the  uppermost  locks.  If  this  option  is  employed,  the  potential  for  increased 
salinity  intrusion  and  necessary  remedial  measures  should  be  addressed,  b.  A  sump  in  the  channel 
upstream  of  the  locks  can  capture  salt  water  for  withdrawal.  It  is  most  effective  for  higher  salinity  water  in 
the  locks;  however,  higher  salinity  in  the  locks  risks  greater  salinity  contamination  of  the  lake  if  any  lock 
water  escapes  the  sump.  A  sump  will  require  some  sort  of  confined  channel  upstream  of  the  lock  to  train 
the  densimetric  flow  out  of  the  lock.  c.  A  salinity  higher  than  0.3  ppt  in  the  uppermost  locks,  including  a 
Syncrolift  lock,  risks  unsatisfactory  buildup  of  salt  water  in  Lake  Gatun  unless  some  means  of  capturing 
most  of  the  salt  water  is  incorporated  in  the  facility,  d.  A  gated  chamber  on  the  lake  side  of  the  uppermost 
lock  offers  the  option  of  draining  salty  water  from  the  lock  and  replacing  it  with  fresh  water  before  opening 
for  ship  passage.  This  technique  provides  positive  salinity  control  at  the  standard  or  Syncrolift  locks  at  the 
cost  of  a  delay  in  the  vessels’  exit  of  the  lock. 


The  following  recommendations  are  made:  a.  Further  consideration  may  be  given  to  recirculating  lift  water 
through  holding  ponds  as  a  means  of  freshwater  conservation,  b.  Further  consideration  may  be  given  to  a 
gated  chamber  at  the  uppermost  lock  entrance  on  each  side  to  reduce  salinity  intrusion  into  Gatun  Lake  if 
either  the  Syncrolift  lock  or  the  holding  pond/recirculation  alternatives  are  considered  for  design,  c.  Tools 
developed  may  be  used  in  the  design  phase  of  any  canal  improvement  plans.  They  should  be  verified  to  a 
comprehensive  field  data  set  of  water  levels,  velocities,  salinities,  and  temperatures  before  use  in  a  design 
effort,  d.  A  numerical  model  of  the  sea  approaches  may  be  developed  to  calculate  the  salinity  there  under 
any  alternative  that  might  alter  the  salinity,  e.  Field  data  collected  in  support  of  further  modeling  should  be 
synoptic  to  the  maximum  degree  possible,  should  include  acoustic  Doppler  current  profiling,  and  should 
include  salinity  meters  capable  of  measurements  to  hundredths  of  parts  per  thousand  concentration. 


